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SUM. MAR  Y 


1  ht*  purpose  ol  this  work  is  to  perform  a  prool'-of-pr inciplo  test  of 
a  unique  eoncepl  for  auxiliary  electric  power  g£^he ration  based  on 
thenri'so  'of  a  piezoelectric  venerator.  The  demonstration  device  con¬ 
sists  01  a  modifier!  St  irling-cycle  engine  which  delivers  pneumatic 
pressure  pulses  to  a  coupler  consisting  of  an  oscillating  mercury 
column  near  ambient  temperature.  The  coupler  replaces  the  fly¬ 
wheel  a  nr!  mechanical  coupling  of  a  conventional  engine -generator 
and  drives  a  hydraulic  load  which  simulates  a  piezoelectric  generator. 

The  work  is  divided  into  two  phases.  Phase  I  includes  piezoelectric 
generator  tests  and  a  system  design  based  on  computer  simulation 
anc:  requisite  mechanical  design  analysis  to  develop  complete  system 
layout  drawings.  Preliminary  piezoelectric  generator  tests  per¬ 
formed  by  Physics  International  were  inconclusive  in  certain  areas. 
Subsequent  testing  at  DWDL  using  an  electrohydraulic  materials 
loading  system  established  the  basic  feasibility  of  the  piezoelect ric 
generator  concept.  Measurements  at  60  Hz  included  power  densities 
to  18  w/in.  efficiency  of  SO";,  to  7r>%  and  continuous  running  time  to 
20  hours. 

The  final  system  design  which  includes  substantial  margins  in 
several  areas,  is  predicted  to  deliver  6  18  v.itts  ( mecha nical)  to  the 
simulated  load  with  an  overall  effic  iency  of  24'r».  \\  hen  coupled 

with  a  piezoelectric  generator,  410  watts  (electrical)  at  J9"u  over¬ 
all  system  efficiency  is  expected.  Safety  and  operational  problems 
with  the  mercury  coupler,  while  not  inhibiting  a  proof-of-principle 
demonstration  of  the  concept,  indicate  the*  need  for  change  in.  anv  proto 
type  system. 


The  originally  planned  Phase  II  efforts  consisted  of  detail  design,  fab¬ 
rication,  and  testing  of  the  engine-coupler-simulated  load  system,  but 
later  was  redirected  toward  life  testing  of  piezoelectric  generators. 
Piezoelectric  generator  feasibility  was  demonstrated  at  ambient 
temperatures  from  -8S°  to  Operating  times  of  more  than 

2000  hours  nt  °000  p-si.  peak  stress  levels  resui  ed  in  about  ^.  perfor¬ 
mance  degradation. 
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FOREWORD 


This  report,  was  prepared  by  the  Donald  W.  Douglas  Laboratories 
Richland,  Washington,  a  .Subdivision  of  the  McDonnell  Douglas 
Astronautics  Company,  under  U.  S.  Department  of  the  Army  con¬ 
tract  DA  AK02-7  1 -C -0405  to  perform  a  proof-of -princ iple  demon¬ 
stration  of  a  500-w  Stirling  engine  pie/.oelectr ic  (STEP/.  1  generator. 
The  work  was  monitored  by  J.  Harrison  Dai.ic-1,  Jr.,  Project: 
Engineer,  U.  S.  Army  Mobility  Equipment  Research  and  Develop¬ 
ment  Center,  Fort  Belvoir,  Virginia. 

This  Final  Technical  Report  covers  the  contract  period  from  >0  June 
1971  through  30  November  1971  for  Phase  I  and  from  26  January  1973 
through  2b  July  1973  for  Phase  II.:  Phase  I  consisted  of  developing 
layout  drawings  for  the  system  and  demonstrating  practicability  of  the 
piezoelectric  generator.  Phase  II  was  originally  planned  to  include 
detail  design,  fabrication,  and  testing  of  a  system  utilizing  a  hydraulic 
load  which  simulates  the  piezoelectric  generator,  but  later  was  redi¬ 
rected  to  provide  2.000-hour  piezoelectric  generator  life  tests  under 
varied  environmental  conditions. 

The  program  was  managed  by  M.  A.  White  under  the  direction  of 
Dr.  W.  R.  Martini.  Piezoelectric  generator  demonstrations  at 
Physics  International  Company  during  Phase  1  were  performed  under 
the  direction  of  P.  Smiley. 
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1.1  P  KOC.K  A  M  O  1  *.l  KCTIV  KS 

The  Ions’- range  object »ve  of  this  proc/am  is  to  develop  a  long-life, 
high-efficiency.  maintenance-free  Pame- fired  engine-generator  as  a 
silent,  portable,  electric  power  source.  The  effort  reported  here 
encompasses  layout  design  of  a  pr.»of-of-principle  svstem  as  well  as 
performance  and  1  i  f<*  testing  of  pivzoelec!  ric  generators.  C'tirrentlv. 
in  both  civilian  and  military  nuxi'ia  rv  power  appi  icat  inns,  gasoline 
or  diesel  engines  operating  rotating  elect  romagnot  ic  genern'ors  are 
almost  universally  used  for  portable  auxiliary  electric  power.  These 
units  are  reasonable  in  cost  but  are  noisv.  relatively  unreliable,  and 
require  regular  knowledgeable  maintenance  of  the  engine.  Particularly 
in  the  Army,  there  is  a  need  for  something  better  at  a  reasonable  price 
Many  applications  in  which  reliability,  low  noise  emission,  and  dura¬ 
bility  under  rugged  use  conditions  arc  important  require  continuous 
electric  power  in  quantities  larger  than  batteries  can  supply  over  long 
periods. 

The  STEP/  concept  has  the  potential  for  satisfying  these  require¬ 
ments,  but  it  involves  certain  innovative  and  unconventional 
approaches.  The  initial  .objective  of  this  program  has  accordingly 
been  established  as  a  proof-ot.-prineiple  reduction  to  practice  of 
the1  key  components.  Included  are  verification  off])  pie/.oe\cct  ric 
generator  feasibility,  (2>  resonant  operation  of  the  coupler  to  ob¬ 
tain  the  required  phase  lag  without  mechanical  linkages,  and 
f  1)  adequate  heat  exchange  capabilities  for  n  simple  displacer- 
regenerator  engine  with  no  extended  heat  transfer  surfaces  for  use 
in  the  nmltihundrcd-waf  t(e)  range. 

In  the  Phase  I  effort,  an  experimental  evaluation  of  the  piezoelectric 
generator  was  conducted  a  nr!  a  layout  design  of  a  proof-of-principle 
system  was  completer!.'  The  system  consists  of  a  modified  Stirling 
engine,  a  -simulated  piezoelectric. generator,  and  a  mercury -filled 
hydraulic  coupler  to  match  engine  output  w  ith  generator  input. 

Basic  feasibility  of  the  piezoelectric  generator  concept  was  established 
in  Phase  1,  but  the  question  of  long- term  degradation  led  to  a  redirected 
Phase  II  effort.  lest  fixtures  suitable  for  life  testing  piezoelrr  t  ru 
generators  was  developed  and  applied  to  2000-hour  tests  of  four  stacks 
under  various  environmental  conditions.  1  lie  life-test  fixtures  consist 
of  a  rigid  loud  frame  which  encloses  an  .SO-disc  piezoelectric  slack  and 
associated  force  and  displacement  measuring  instrumentation.  Seventy 
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of  the  discs  are  operated  as  a  driver  stack  from  a  high  voltage  ac  power 
supply.  The  resultant  stress  fluctuation  produces  ac  power  in  the 
10-disc  generator  stack;  power  density  and  efficiency  are  measured. 
Long-term  performance  degradation  was  not  severe. 

77  2  PHASE  I  PIE '/OE  LEG  TRIC.  GENERATOR  I^STING 

Piezoelectric  generators  are  novel,  partly  because  practical  materials 
have  only  recently  come  into  use.  In  addition,  the  generator  requires 
pulsating  mechanical  power  at  high  force  and  low  displacement  which 
is  not  a  typical  output  characteristic  for  prime  movers.  A  ‘ 1 7 .  T’i.  effi¬ 
cient  piezoelectric  generator  driven  by  a  similar  piezoelectric  motor 
has  recently  been  reported  (Reference  It. 


This  generator,  however,  operates  at  very  high  frequencies  and  low 
stresses  compared  with  the  STEP/,  ope  rating  regime.  Acknowledged 
experts  in  the  piezoelectric  materials  field  differ  widely  in  their  eval¬ 
uation  of  practical  operating  stress  levels  and,  therefore,  powe  r  density. 


Because  of  the  relative  uncertainty  of  piezoelectric  generator  feasi¬ 
bility  compared  with  that  of  Stirling  engines,  the  -Phase -I  effort 
included  component  testing  of  piezoelectric  stacks.  Testing  of  two 
considerably  different  piezoelectric  materials  was  originally  per¬ 
formed  under  subcontract  to  Physics  International,  (results  reported 
in  Appendix  A),  but  limitations  of  the  stack  stressing  equipment 
prevented  extended  operation  as  well  as  efficiency  and  temperature- 
rise  measurements.  These  data  have  now  been  obtained  at  DUDL 
and  are  presented,  along  with  supporting  analysis,  in  Appendices  13  and  C. 

Piezoelectric  generator  operation  is  shown  in  Figure  1-1.  The  disks 
ave  made  from  a  lead  ziromate -titannte  ceramic  and  metallized  on 
both  flat  surfaces  for  electrical  contact.  Warm  disks  are  polarized 
(analogous  to  magnetizing  a  permanent  magnet!  by  applying  a  high 
voltage  to  align  the  dipoles.  These  arc  stacked  with  alternating 
polarity  (electrically  in  parallel!  with  one  side  connected  to  the  load, 
the,  other  to  ground.  Compression  of  the  stack  produces  charges 
delivered  as  load  current  in  one  direction,  followed  by  a  reversed 
'  polarity  and  current  when  compression  is  relaxed.  Theory  and 
experiment  show  that  typically  only  about  half  the  applied  strain 
energy  tjrodtiees  a  charge.  I  he  balance  is  stored  elastically,  and 
in  real  system,  is  utilized  to  effect  the  compression  stroke  on 
the  engine*.  A  portion  of  the  generated  power  is  dissipated  in  stack 
losses.  This  must  be  minimized  to  maintain  an  acceptable  efficiency 
and  maintain  temperature  rise  in  the  stack  within  toloiable  limits. 
Experiments  confirm  the  feasibility  ol  these  objectives. 


preliminary  experimental  results  obtained  at  Physics  International 
are  summarized  in  Figure  1-2  with  measurements  normalized  to 
specifie  output  power.  Early  stack  stress  fluctuations  to  5500  psi 
u.-re  obtained  using  a  piezoelectric  driver.  Extrapolation  of  these 
data  is  indicated  as  the  original  projection  in  Figure  1-2.  The 
tests  rondurted  as  part  of  the  Phase  I  ellort  used  a  rotary  valve 
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and  high-pressure  hydraulic  pump  to  apply  stress  fluctuations.  Two 
'types  of  pie/.oceramic  materials  were  tested.  PZT -5H  is  a  piezo- 
electrically  soft  material  with  a  high  power  density  at  a  given  stress 
level,  but  with  relatively  low  stress  limitations  and-;  high  internal 
loTTTses.  '  As  indicated  in  Figure  1-2,  P/.T-5H  exhmited  degradation 
with  short-term  tests  near  5000  psi  compressive  stress.  The  other 
material  tested  was  LTX-1  (equivalent  to  PXT-4),  a  pie/.oelectri- 
cally  hard  material  with  basically  opposite  properties  to  those  of 
P/T-5H.  This  stack  also  exhibited  an  apparent  leveling  off  in 
power  density  around  9000  psi,  but  this  was  attributed  to  high 
pressure  binding  in  the  drive -piston  seal. 

Some  of  the  pertinent  data  obtained  at  DWDL  utilizing  an  electro- 
hydraulic  mate  rials  -loading  system  are  presented  in  Figure  1-3. 
These  data  were  obtained  using  sinusoidal  stressing  at  60  Hz  with  a 
minimum  compression  ol  1000  psi,  and  / 000  to  13,  000  psi  maximum 
compression. 

Experimental  power  density  follows  the  predicted  curve  closely  to 
8000 -psi  stress  fluctuation,  with  decreased  performance  at  higher 
levels.  Following  these  measurements,  a  20-hour  continuous  run 
at  12,  000  psi  stress  fluctuation  exhibited  continuous  and  irreversible 
degradation  at  a  decreasing  rate.  A  similar  continuous  run  at  6000 
psi  stress  fluctuation  three  weeks  later  showed  only  an  initial, 
re ve rs Ible  de gr adat ion  while  the  hydraulic  driver  approached  thermal 
equilibrium.  Inductance  was  then  used  in  the  load  circuit  to  com¬ 
pensate  for  capacitance  inherent  in  the  generator.  Power  density 
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Figure  1-1.  Piezoelectric  Generator  General  Arrangement 
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OUTPUT  ELECTRIC  POWER  (W  /IN  OF  PZ  MATERIAL) 
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Figure  1-2.  Measured  PZ  Generator  Output 
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Figure  1-3.  Piezoelectric  Generator  Test  Results 

increased  by  aTactor  of  2.  5  as  compared  with  a  predicted  value  of 
with  .  A piezoelectric  generator  efficiency  of  67  %  was  measured, 
with  optimum  inductance  in  the  circuit  at  6000  psi  and  a  power  den¬ 
sity  ot  15  w/cu  in.  r 

These  data  show  that  an  efficient  piezoelectric  generator  with 
reasonable  power  density  is  possible.  Excellent  analytical  corre¬ 
lations  with  experiment  give  confidence  in  designing  generators  for 
specific  applications.  Further  experiments  will  eventually  be 
required  to  determine  the  effect  of  disk  geometry  and  surface  finish 
on  long-term  operating  properties  and  allowable  stress  levels. 

‘  I.  •  '  _ 

1.3  LAYOUT  DESIGN  OF  PROOF-OF-PBINCIPLE  SYSTEM 

The  Proof ~of- principle  system  is  designed  to  show  that  the  basic 
concept  shown  in  Figure  1-4  is  feasible;  the  system  utilizes  a  simu¬ 
lated  piezoelectric  generator.  The  Stirling  engine  produces  pres- 
sure  pulses  from  heat.  After  a  phase  shift  by  Inertia  in  the  Coupler, 
these  pressure  pulses  can  be  applied  to  a  piezoelectric  generator  to 
3  electricltV-  Gas  in  the  Stirling  engine  acts  as  a  spring  and 

a  egative  damper,  one  that  produces  instead  of  consuming  power. 
The  coupler  acts  as  a  mass,  and  the  piezoelectric  generator  acts 


5 


i  STIRLING  COUPLER  PIEZOELECTRIC 

•  ENGINE  (INERTIAL)  GENERATOR 


V 


DP.I'l  .Vt  V.b  Nt  •  I’tlrtf  It  PlSTilM 


Figure  1-4.  Basic  Concepts  of  the  STEPZ  Generator 


as  another  spring  combined  with  a  conventional  damper.  Figure  1-4 
shows  how  forces  from  the  engine  and  generator  act  on  the  coupler 
which  serves  as  an  inertial  energy  storage  device  analogous  to  the 
flywheel  in  a  conventional  engine..  Figure  1-5  shows  the  principal 
features  of  the  proof-of-principle  system  to  demonstrate  the  basic 
concept.  . 

The  Stirling  engine  is  similar  to  other  silent,  long-lived  power 
sources  developed  and  tested  at  DWDL.  The  engine  consist-:  of  a 
lightweight  displacer  which  oscillates  inside  the  engine  cylinder 
with  a  five-mil  jradial  clearance  between  cylinder  wall  and  displacer. 
This  separation  is  maintained  by  a  flexure  at  the  hot  end  and  a  bear¬ 
ing  (part  of  the  displacer  drive)  at  the  cold  end.  This  bearing  may 
ultimately  be  replaced  by  a  flexure.  The  gap  acts  as  a  gas  heater, 
regenerator,  and  gas  cooler.  The  use  of  this  gap,  instead  of  many 
small  tubes  manifolded  into  regenerator  assemblies  as  in  conventional 
Stirling  engines,  greatly  simplifies  the  engine  anti  significantly  reduces 
co s t  for  small  Stirling  engines.  Oscillation  of  the  displacer  heats 
and  cools  the  confined  helium  and  creates  the  pressure  pulses.  An 
electric  motor  for  startup  and  fine  frequency  control  moves  the  dis¬ 
placer  through  a  crank  mechanism.  These  components  will  probably 
be  replaced  by  an  electromagnetic  drive  when  prototype  development 
is  achieved.  Once  operating  speed  is  reached,  the  gas  pressure 
difference  applied  to  the  displacer  drive  piston  varies  over  a  cycle 
in  such  a  way  as  to  apply  power  to  the  displacer  drive.  If  friction 
is  not  too  great,  the  engine  will  run  itself.  The  crank  mechanism 
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F;v  >re  1-5.  Proof-of-Principle  Test  Setup 
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is  useful  to  control  motion  of  the  displacer  to  fuH  stroke  without 
danger  ol  hitting  the  ends.  A  similar  smaller  engine  operated  at 
20  Hz  and  =>  watts  with  a  crank  and  flywheel  ior  over  5000  hours  at 
DWDL. 

‘ !  '  '  ,  .  •  .  '  -  ■  ■  V  .  ' 

The  coupler  considered  for  the  present  study  isjt mercury- filled 
— 4+rbe-with  diffuse  rs  at  each  end.  A  corrugated  cflaph  ragm  separates 
the.  mercury  from  the  engine  gas  at  one  end  and  from  the  hydraulic 
fluid  in  the  gone  ra  tor  simulator  at  the  other.  The  necked-down 
configuration  produces  the  effect  of  a  much  larger  mass  than  is 
actually  present  (similar  to  a  flywheel  attached  to  a  speed-increase’* 
go  a  r  > .  therefore  the  name  "fluid  flywheel.’  For  the  reference 
design,  the  coupleir  is  96%  efficient.  One  inconvenience  in  this  test 
is  the  reaction  forces  of  the  coupler  which  require  rigid  anchoring 
for  the  test.  In  a  practical  system,  imbalanced  forces  are  elimin¬ 
ated  by  two  opposing  couplers  or  inertial  members.  A  liquid 
coupler  can  also  be  coiled  to  reduce  system  dimensions.  Another 
-maior  change  anticipated  for  a  practical  generator  is  in  the  dis¬ 
placer  drive  system.  The  present  crank  and  flywheel  mechanism  and 
sliding  seal  will  be  replaced  by  a  linear  electromagnetic  drive 
mechanism  with  full  flexural  support.  A  dynamic  jcoil;  motor  has 
been  concbptualiy  designed  with  27%  efficiency  (Appendix  E)  at  the 
ope  -watt  level  lor  a  50  \v(e '  STEP/,  generator.  This  approach  has 
j  more  basic  simplicity,  which  is  a  fundamental  objective  of  the 
1  STEP/  effort,  than  does  the  proof-of- principle  design. 

The  generator  simulator  is  a  liquid  pump.  Pump  compliance  and 
pressures  in  the  accumulator  and  reservoir  are  chosen  to  simulate 
actu,»u  of  the  piezoelectric  generator.  Power  output  of  the  simulator 
can  be  checked  by  measuring  the  steady  .flow: and  pressure  drop  at 
the  throttle  valve. 

Figure  1-6  shows  the  proof-of-principle  device  based  on  current 
designs.  The  coupler  is  mounted  vertically  on  a  firmly  anchored 
I-beam.  The  engine  operates  horizontally  and  is  attached  to  the  top 
of  the  beam.  The  simulated  piezoelectric  generator  is  on  the  side. 
The  lest  system  is  about  5b  in.  high.  24  in.  wide,  and  24  in.  deep. 


1. 4  DESIGN  GOAL  CONCEPT 

The  proof-of-principle  device  is  not  intended  to  be  a  prototype  of  a 
practical  silent  power  source.  Once  the  proof-of-principle  is 
demonstrated,  a  number  of  changes  and  additions  must  be  marie  to 
develop  a  practical  prototype.  These  are  expected  to  be: 

1.  A  piezoelectric  generator  instead  of  the  simulation. 

2.  A  counterbalanced  coupler  to  eliminate  reaction  forces. 

5.  A  simplified  displacer  drive  mechanism. 
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Figure  1-6.  Assembly  Drawing  of  Proof-of-Principle  Device 
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4.  Automatic  control  to  give  desired  output  voltage  over 

the  full  power  range.  ; 

5.  Burner  heat  source.  j* 

6.  Radiator  heat  exchanger. 

1.5  PIEZOELECTRIC  GENERATOR  LIFE  TESTS 

Life  tests  were  conducted  to  establish  the  effect  of  extended  operation 
(2000-hour  goal)  on  several  piezoelectric  stacks  with  different  surface 
finishes  at  various  environmental  temperatures.  Four  life  test 
assemblies  were  fabricated  similar  to  an  IRAD  prototype  system  which 
was  designed  and  built  in  the  interim  between  Phase  I  and  Phase  II. 

Each  test  fixture  consists  of  a  stack  of  eighty  0.  5-in.  diameter  by 
0.  040-in.  thick  piezoelectric  discs,  restrained  within  a  rigid  load 
frame,  together  with  force  and  displacement  measuring  instrumenta¬ 
tion.  Electric  excitation  at  60  Hz  is  applied  to  70  discs  of  each  stack, 
which  generates  the  driving  force  to  excite  the  10-disc  generator  stack. 
Specified  test  conditions  for  the  four  stacks  are  summarized  in 
Table  1-1.  |  .  ■  /  ; 

In  addition,  the  IRAD  prototype  test  fixture  was  reassembled  in  a  new 
configuration  and  life  tested  concurrently  with  the  contract  stacks  at 
room  temperature.  The  IRAD  stack  is  unique  in  that  it  utilizes  a 
60-disc  driver  stack  and  the  generator  stack  is  only  0.  25  in.  diameter. 
This  enables  the  generator  to  experience  the  full  desired  stress  swing 


Table  1-1 

SPECIFIED  LIFE  TEST  CONDITIONS 


j  ;  . 

5  ■ 

Maximum 

Minimum 

Vendor  Specified 

!  j  Ambient  Temp 

Str es  s 

Stres  s 

Disc  Surface 

|  |  Stack  No.  (°  F) 

(psi) 

(psi) 

Finish** 

j'  1  <20* 

i 

9000 

As  low  as 
pos  sible 

±0.  0005  in.  flatness 
(standard  double  lap) 

2 

<20- 

9000 

±0.  00005  in.  flatness 
(optical  lap) 

3 

Approx  70 

9000 

tO.  0005  in. 

'  4 

Appr-.x  150 

9000 

±0.  0005  in. 

'••■After  1000  hours  of  life  testing  at  cold  point  of  standard  freezer, 
ambient  test  ten-  >erature  will  be  reduced  to  <-25°F  for  100  hours 
of  operation.  A'.  er  this  100  hours,  life  test'  will  continue  at  <20°  F 
ambient.  . 

■:*As  specified  in  fiulton  Industries,  Inc.  Quotation  No.  13K2,  12/15/71. 
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of  1000  to  9000  psi,  in  contrast  to  the  contract  stacks  where  the 
generator  stress  swing  is  typically  around  5000  to  9000  psi  as  a  result 
of  excessive  stack  compliance. 


All  stacks  have  survived  ,the  life  test  with  total  times  of  2300  to  3300 
ho.urs.and  are  continuing  to  accumulate  ho-urs,  v*£lh  the  exception  of 
stack  number  3  which  suffered  irreparable  damage  as  a  result  of 
accidental  overstressing  during  a  preload  operation.  Several  isolated 
discs  in  'driver  sections  suffered  electrical  breakdown  during  •'he  life 
test,  but  in  each  case,  the  appropriate  10-disc  segment  was  removed 
from  the  circuit  and  operation  continued.  No  breakdown  or  severe 
degradation  occurred  with  any  generator  stack. 


Quality  control  with  commercially  supplied  piezoelectric  discs  is  a 
practical  problem.  Standard  lapped  discs  were  microscopically 
reasonably  smooth,  hut  polarization  after  lapping  caused  warping  with 
0.  002-  to  0.004-in.  deviation  from  flatness.  The  optically  lapped 
discs  were  much  better,  but  only  about  as  smooth  as  the  vendor 
specified  finish  for,  standard  discs,  and  parallelism  of  the  flat  surfaces 
was  poor.  The  warped  discs  flatten  with  a  relatively  small  preload, 
but:  compliance  of  the  standard  discs  was  substantially  higher  than  for 
the  optically  lapped  discs.  Those  stacks  received  from  the  supplier  in 
a  permanently  bonded  »tate  exhibited  even  higher  compliance  and  were 
not  "used  for  life  testing.  High  compliance  is  a  severe  problem  for  the 
driver  stack s  asN  used  for  life  testing,  but  may  actually  be  beneficial 
in  a  pure  generator  application.  To  maintain  operating  voltages  at 
manageable  leive|!s,  many  relatively  thin  discs  are  used.  The  resultant 
large  number  of  interfaces  is  the  principal  source  of  high  compliance, 
with  the  result  that  much  of  the  strain  energy  produced  is  not  applied 
to  the  generator  stack. 


One  result  of  nearlv  all  stack  testing,  or  even  preloading  to  a  sub¬ 
stantial  axial  load,  is"  that  many  of  the  discs  crack  into  two  or  more 
:  pieces.  Based  pn  MTS  .tests'  using  discs  before  and  after  cracking 
;  occurred,-  no  apparent  performance  degradation  results  from  this 
disturbing  phenomenon.  Some  of  the  life  test  stacks  were  found  to  have 
cracked  discs  after  assembly. 


Through  7  September  1973,  the  1RA1)  stack  had  accumulated  3338  hours 
of  excitation  for  a  total  of  0.72  1  x  10  cycles.  Of  the  3  118  hours 
accumulated  after  the  final  load  was  established,  the  reduction  in  nor¬ 
malized  output  power  was  less  than  5".>,  discounting  occasional  statis¬ 
tical  or  other  variations.  The  normalization  applied  to  the  power 
output,  based  on  theoretical  expectations,  was  to  divide  the  power 
density  by  the  square  of  the  applied  stress  swing.  Efficiency  data 
showed  considerable  scatter,  especially  about  900  to  1200  hours  after 
life  tests  began,  but  the  efficiency  trend  was  basically  constant. 
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n,av  :  v  Prohibitive  £££  "co^p  VeU^^t’"* 

1 80- force  low-displacement  characteristic  of  piezoelectric  generators 
P^sents  significant  problems  as  well.  If  successful  application  of 
Stirling  engine  piezoelectric  generators  is  made,  it  mav  initiallv  hP 

very  small  power  levels  where  the  piezoelectric  generator  is  unique  in 
maintaining  a  high  efficiency.  ^  s  unique  m 


'  ;  -  Section  2 

i _  INVESTIGATION  AND  DISCiJSS^fN 

This  section  discusses  the  analytical  methods  applied  to  all  the 
components  ot  the  system,  how  these  methods  were  used  to  arrive  at 
a  reference  design  for  the  proof-of-principle  test  device,  a  piezo¬ 
electric  generator  design  to  meet  the  same  requirements,  and 
mechanical  design  considerations  incorporated  into  the  proof-of- 
principle  test  device. 

2.1  ANALYTICAL  METHODS 

Concurrent  with  this  work,  a  computer  simulation  has  been  developed 
for  a  Stirling  engine,  coupler,  and  simulated  load  system  in  an 
independent  research  and  development  program.  Based  on 
experience  with  this  simulation,  a  reference  design  was  chosen  and 
the  mechanical  design  layouts  were  completed. 

...  •  ’  1  ’  •  " 

2.  1.1  Engine  Analysis 

Figure  2-1  shows  a  simplified  schematic  of  the  engine.  The  cylinder 
wall  consists  of  discrete  heater,  cooler,  and  regenerator  sections. 
Inner  surfaces  of  the  heater  and  cooler  are  held  fixed  at  the  heat 
source  and  sink  temperatures,  respectively.  The  regenerator  section 
has  an  approximately  linear  temperature  profile  between  the  heater 
and  cooler  temperatures.  As  the  displacer  reciprocates  in  the 
cylinder,  engine  gas  is  shuttled  between  the  hot  and  cold  ends  through 
the  thin  annula r .clearance.  The  power  piston  follows  displacer 
motion  with  a  phase  lag  near  90s.  The  expansion  stroke  occurs 
while  the  displacer  moves  from  mid-stroke  to  bottom  dead  center  and 
back  to  mid-stroke  (while  most  of  the  engine  gas  is  in  the  hot  region). 
Conversely,  power  piston  compression  occurs  while  the  displacer 
moves  between  mid-stroke  and  top  dead  center  (gas  mostly  in  the  cold 
region).  Because  hot  expansion  work  output  is  greater  than  cold 
compression  work  requirement,  net  power  is  delivered  by  the  engine. 

i  Various  approaches  to  analyzing  modified  Stirling -cycle  engines  have 
been  developed  at  DWDL  over  the  past  five  years.  The  most  exten¬ 
sively  used  is  the  isothermal  analysis  in  which  gas  volumes  in  the 
engine  are  assumed  to  have  a  steady  temperature  throughout  one  cycle 
This  is  an  integrated  approach  in  which  all  computations  are  time 
averaged  ove  r  one  cycle.  All  loss  mechanisms  a  re  assumed  to  act 
independently  so  that  the  principle  of  superposition  can  be  applied. 
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Figure  2-1.  Simplified  Engine  Schematic 
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Ain-jlysis  can  he  reduced  to  the,  solution  of  a  set  of  algebraic  equations. 
1  hi;s  approach  has  the  advantage  of  short  computation  times  and  is 
well  suited  to  parametric  survey  and  optimization  analyses.  It  has 
been  used  with  good  success  in  accurately  analyzing  engines  lor  the 
National  Ilea rt*  and  Lung  Institute. 


The_jLiih.ee  basic  type  of  program  is  called  nordsothe^hna  ■  because  gas 
volumes  in  the  engine  have  a  changing  temperature  throughout  one 
cycle.  The  computer  program  utilizes  a  time-step  solution Which 
models  the  system  with  simultaneous  interactions  and  inherently 
includes  nonisot he r mal  heating  and  cooling  effects  caused  by  pressure 
changes.  The  accuracy  of  the  results  is  limited  only  bv  the  accuracy 
ol  the  analytical  model  and  the  time  .step  size,  which  can  ho  arbitrarily 
small  in  principle.  Computer  running  times  are  typically  two  or  three 
orders  of  magnitude  longer  for  the  latte,  r  program  because  a  full  set 
of  calculations  must  he  -made  for  each  time  step  in  a  cycle  and  the 
solution  must  continue  through  se vc  ra  1  "cycles  to  eliminate  transients. 


Computer  programs  representing  both  types  of  analytical  models 
have  been  used  in  the  design  of  the  STEP/,  proof -of- principle  system. 
The  isothermal  program  requires  about  seven  seconds  per  case  and 
is  set  up  for  running  parametric  surveys  with  minimum  effort.  It 
has  proved  to  he.  satisfactory  for  optimizing  all  parameters  except 
compression  ratio.  If  the  latter  is  restricted  to  values,  known  by 
experience  with  the  nonisothormal  simulation  to  he  reasonable,  the 
isot  he  rma  r  prog  rain  j  ca  n  bo '  used  with  confidence  to  design  new 
engine  configurations,  and  in  practice,  is  the  only  feasible  means 
for  so  doing.  The  nonjisot he  mini  simulation  includes  dynamics  of 
the  entire  system  and  provides  t  he  best  a'  /liable  accuracy  for 
predicting  operation;  of  a  given  systc*  •  Comparison's  are 

made  in  Section  .i.  2.  I  between  the  isot  tie  final- and  nonisothormal 
performance  prediction's. 


The  isothermal  program  was  developed  to  operational-status  and 
used  with, extensive  pa  rametric' surveys  to  establish  the  STEP/, 
engine  geometry  prior  to  initiation  of  this  program.  The  simu¬ 
lation' was  then  utilized  during  phase  1  to  verify  the  earlier 
engine  design  and  to  re  tine  the  operating  parameters  and  portor- 
mance  specifications  where  necessary. 


i.  I.  1.  I  Isothermal  Analysis 

In  this  analysis,  g.is  temperatures  in  a  given  region  are  assumed 
not  to  vary  in  temperature  over  one  cycle.  Average  hot  and  cold 
gas  temperatures  are  determined  by  an  iterative  process  based 
on  limited  heat  transfer.  power  output  is!  computed  from  the 
S<  hmidt  equation  i  Reference  which  is  an  analytical  solution  of  the 
pro s  u  re  -  vol  no  ie  i  ut  e  g  r.t  1  over  one  eyclo.  Sinusoidal  motion  of  the 
regenerator  and  power  piston  with  a  specified  lixed  phase  angle  is 
assumed,  which  is  very  realistic  tor  conventional  Stirling  engines  with 

lb 
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mechanical  phasing.  In  the  STEP/,  system,  both  power  piston  and 
displacer  move  almost  sinusoidally  but  the  phase  angle  between  the 
two  is  not  fixed. 

Hfcnt  transfer  between  metal  and  gas  in  the  annular  ^clearance  is 
de-fce-rmined  by  assuming  laminar  gas  flow'.  -The  Heat  transfer 
coefficient  depends  on  annular  clearance  and  not  on  gas  velocity. 
The  heat  transfer  calculational  approach  is  given  in  some  detail  in 
Appendix  D. 


Imperfect  heat  transfer  in  the  heat  exchangers  reduces  the  average 
hot  gas  temperature  and  increases  the  average  cold  gas  temperature. 


The  basic  engine  heat  requirement  (second 


from  output  power  by  assuming  the  engine  is  an  ideal  heat  engine 
with  Ca  root  efficiency  computed  on  the  basis  of  average  hot  and 
cold  gas  temperatures.  Another  set  of  heat  requirements  (or  losses), 
only  one  of  which  is  significant  in  the  STEP/,  engine,  applies  when 


the  engine  is  running.  Less -than-perfect  heat  transfer  in  the 


regenerator  region  (main  reheat  loss)  contributes  the  primary  loss 


mechanism  in  the  STEP/  engine.  Main  reheat  loss  results  from 
not  regene  ratively  supplying  all  the  heat  necessary  to  increase  the 
gas  temperature  from  the  cold  to  hot  values  as  gas  flow’s  through 
the  regenerator.  Heat  which  is  not  supplied  regene  ratively  must 
be  supplied  by  the  heater.  Calculation  of  main  reheat  loss  is  detailed 
in  Appendix  D.  The  rest  of  the  thermal  losses  apply  whether  the 
engine  is  running  or  net.  These  include  insulation  loss,  cylinder 
and  displacer  w;all  conduction,  and  gas  conduction  and  radiation 


through  the  displacer. 


In  high  power  engines,  mass  flow  around  the  regenerator  is  a 
critical  parameter.  In  the  isothermal  analysis  program,  the  maxi¬ 
mum  and  minimum  mass  of  gas  in  the  hot  and  cold  gas  space  is 
computed  from  a  formula  which  considers  the  sinusoidal  motion 
of.' both  the  displace  r  and  the  power  piston,  fixed  gas  temperatures, 
fixed  total  mass  of  gas,  and  an  instantaneously  uniform  engine  pres¬ 
sure.  Mass  flow’  calculated. in  this  way  is  only  approximate  and 


leads  to  a  significant  error  in  main  reheat.  This  does  not  have  a 
■major  Impact  on  the  optimization  capabilities,  w’hich  represent  the 
prime  usage  of  this  program.  Windage  power  is  based  on  the  fur¬ 
ther  approximate  assumption  that  gas  flows  past  the  displacer 
sinusoidally,  but  windage  calculations  are  nevertheless  very 


accurate. 


Z.  1.  1.2  Non'.isothc  rrnal  Simulation 

The  nnni  sothe  mini  dynamic  simulation  overcomes  most  of  the  basic 
limitations  imposed  by  the  isothermal  model.  -Two  basic  versions  of 
I  he  simulation  are  available  for  specialized  applications.  The  simple 
version,  like  the  isothermal  program,  calculates  performance  of  the 
engine  only.  Power  piston  anrl  displacer  motion  are  specified  as 
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•  sinusoidal  'with  a  fixed  phase  lag.  This  solution  stnbiliz.es  most  rapidly 
because  only  thermal  processes  are  alloefed  by  transients.  A  more 
complete  simulation  includes  an  equation  of  motion  for  the  displace r 
and  coupler  and  a  generator  simulation.  Because  the  displacer 
I _ equat  ion  of  motion  includes  flywheel  torque,  n^fsexte  r  nelly  driven  dis¬ 

placer  can  be  simulated  by  specifying  a  large  moment  of  inertia  for 
the  flywheel.  While  the  nor.isothe  nival  simulation  is  a  substantial 
1  improvement  over  previous  techniques,  it  remains  incomplete  in  two 
general  areas. 

The  major  thermal  inaccuracy  involves  specification  of  linear  temper¬ 
ature  profiles  in  the  cylinder  and  displacer  sidewalls./  Calculations 
showed  that,  on  the  average,  during  the  cycle,  some  portions  of  the 
cylinder  and  the  displacer  walls  in  the  regenerator  were  losing  more 
heat  tha  iv  they  we  re  receiving,  and  no  provision  Was  made  for  supply¬ 
ing  this  heat  from  the  heat  source.  Other  portions  of  the  regenerator 
we  re  gaining  more  heat  I  turn  the  v  were  losing  and  no  provision  was 
' made  for  conducting  this  heat  to  the  heat  sink.  Figure  1-1  shows 
the  average  gas  temperatures  for  the  10  nodes  of  the  displacer  in 
comparison  with  the  assumed  temperatures  of  the  engine  for  the 
reference  design.  In  the  hot  space,  heater,  cooler,  and  cold  space, 

.  it  is  expected  that  the  average  gas  temperatures  will  be  different 
than  the  metal  temperature  because  of  the  net:  heat  flow  through  the 
engine.  However,  in  the  regenerator,  it  is  assumed  that  there  is 
no  net  heat  flow  between  the  gas  and  the  walls.  For  the  displacer 
this  appears  to  be  true.  In  the  case  of  the  cylinder  wall,  for  the 
assumed  temperatures  to  be  true,  additional  heat,  must  be  supplied 
from  the  heater  and  be  conducted  along  the  cylinder  wall  to  supply  • 
nodes  ?.nnd  4.  In  the  same  way,  additional  cooling  must  be 
arranged  for  nodes  6  and  7.  Under  such  conditions,  an  accurate 
heat  balance  was  not  possible.  However,  the  effect  on  overall 
performance  is  . very  minor.  In  actual  operation,  experience  with 
heart  engines  shows  that  an  S-shaped  wall  temperature  profile  will 
result  with  minimal  effect  on  engine  thermal,  requirements.  The 
baste,  heat  requirement  for  the  engine  is  determined  by  ,/pdV  in 
the  hot  space  plus  a  calculation  of  the  main  reheat  loss.  These  flu 
not  depend  on  the  computed  heat  transfer  from  and  to  the  nodes. 

The  problem  could,  be  remedied  by  adding  nodes  in  the  metal  walls 
of  the  regenerator,  each  with  an  appropriate  heal  capacity  and 
thermal  conductance  to  adjacent  -nodes,  but  the  simulation  then  would 
require  reprogramming  for  a  larger  computer  and  become  costly  to  run. 

The  other  thermal  inaccuracy  concerns  the  uncertainty  of  heat  transfer 
rate  at  the  end  walls  in  the  hot  and  cold  regions.  A  conservative 
assumption  is  now  being  used  in  both  hot  and  cold  spaces.  It  is  assumed 
that  no  bulk  .mixing  occurs,  and  heal  transfer  occurs  through  stagnant 
gas.  Only  with  experimental  data  can  this  inaccuracy -be  assessed. 

A  unique  method  of  computer  simulation  has  been  developed  which 
allows  a  much  larger  time  step  than  is  stable  for  a  conventional 
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.simulation  scheme,  t  his  improvement  lias  allowed  limited  use  of 
.lonisothcnnal  parametric  survey  techniques  for  this  report. 

Some  inore  subtle  differences  between  the  isothermal  and  nonisothe rmal 
simulations  require;  discussion.  In  the  isothermal  simulation,  second 
law  heat  input  is  determined  by  applying  the_Ca  rn^eqaai  ton  to  the  time- 
aveTnged  ga s  temperatures.  This  is  only  an  approximation  at  best, 
because  the  nature  ,of  the  temperature  variations  over  one  cycle  have  an 
impact  on  heat  transfer.  For  ('sample,  in  the  STEP/,  reference  design 
case  using  the  nonisothe  rma  l  simulation,  hot  gas  temperature  at  one 
point  reaches  GO '  above  the  heat  source  temperature.  This  results 
in  some  heat  flow  from  the  gas  to  the  heater  which -must  be  made  up 
during  other  parts  of  the  cycle.  The  net  result  is  that  most  heat  is 
transferred  tb  the  gas  when  it  is  at  be  low-ave  rage  temperatures. 

The  Effective  ‘temperature  for  determining  second  law  heat  input  from 
the  Carnot  relationship  is,  therefore,,  below  the  time-averaged  gas 
teinpe ratu to ,  and,  for  this  reason,  the  isothermal  program  gives  a 
low  value  for  second  law  heat  input.  However,  the  nonisothe  rmal 
simulation  can  exactly determine  the  value  by  numerically  intog  rating 
the  pressure  with  respect  to  volume  in  the  hot  region  over  one  cycle. 
Work  output  is  similarly  calculated  by  integrating  cold  region  pres¬ 
sure  With  respect  to  total  gas  volume.  Main  reheat  is  determined 
by  the  same  basic  equations  as  in  the  isothermal  program,  except 
that  it  can  be  computed  more  accurately  because  it  is  based  on 
instantaneous  mass  flow.  This  results  in  a  significant  increase  in 
calculated  main  reheat,  as  discussed  in  Appendix  IV 

.  '  . 

In  the  more  complete  engine  simulation,  displacer  piston  motion  is 
determined  by  a  force  balance  on  the  displacer  which  includes  such 
effects  as  drive  chamber  pressure,  flywheel  torque,  I  fiction,  windage, 
gravity,  and  spring  forces.  Flywheel  torque  is  determined  by  solving 
a  torque  balance'  equation  on  the  tlywheel. 


The  equation  for  displacer  dynamic  force  balance  is  (Figure  d-Sh 
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angle  between  drive  rod  and.  connecting  rod 
frictional  drag  on  drive  rod 

reciprocating  mass  of  displacer  «. 

.  —  ■ 

linear  acceleration  of  displacer 


Tito  connect  ing  rod  force,'  F^.,  is  determined  by  a  torque  balance'  on 
the  flywheel 

^C1  rank  sin  (0  -  ?>)  -  Frictional  Torque-*  External  Torque  10  (2) 


where 

r  ,  crank  arm  radius 
crank 

0  angular  displacement  of  connecting  rod  from  top  dead  center 

I  -  flywheel  moment  of  inertia 

•  *  \  , 

;  b  angular  acceleration iof  the  flywheel 

/  i  :  ‘  . 

Equation  1  is  solved  for  displacer  acceleration,  Xj^g,  at  time  t,  and 
this  value  is  used  to  determine  the  velocity,  X^^.  and  displacement 
from  midpoint,1  at  time  t  4  At  by 

xdis (t  "  At\  xdis  +  xnis  (t)  A1 

Xl)IS(  t  f  At)  “  X1)I.S  ,t-  ‘i  XI)LS  At  +  X1)IS  ^  At  (■ 


The  external  torque  term  in  Equation  2  represents  the  torque  applied 
by  an  electric  motor.  If  this  is  y.erp,  there  is  a  net  force  balance?  on 
the  displacer  at  any  time  which  moves  the  displacer  in  some  manner. 
It  is  desirable  to  have  sinusoidal  motion  with  the  displacer  leading 
the  coupler  by  approximately  40“. 


2.  1.  2  Couple  r  A  nalys  is 

Figure  2-4  shows  the  arrangement  of  the1  fluid  coupler.  The  basic 
element  is  the  fluid  inertia  column  contained  within  the  coupler  tube. 
This 'inertia  column  permits' opt imum  engine  and  load  performance 
by  eliminating  the  requirement  of  instantaneous  engine  and  load  pres¬ 
sure  equality.  This  is  marie  possible  by  the  large  pressure ‘differ¬ 
ence  required  to  accelerate  the  liquid  column  at  the  velocities  and 
displacements  of  inti*  rest. 
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i  he  fluid  in  contained  at  each  end  of  the  coupler  by  relatively  large 
diameter  corrugated  metal  diaphragms.  Diffusers  are  provided  at 
each  end  to  minimi/.e  coupler  tube  entrance  and  exit  flow  losses. 

1  he  corrugated  diaphragms  are  thin  and  cannot^Jupport  a  significant 
-pee  s' sure  difference.  With  reference  to  Figure  2-4,  the  equation  of 
motion  of  the  coupler  is 

pe  \c  “  pl  Ac  "  pc  Lc  Ac  *c  _  1/2  cc  Ac  Pc.  ^c;2  "  0 

where 

Pj.  engine  pressure 

A^.  coupler  tube  cross-sectional  area 
Pj  load  pressure 

p^.  coupler  fluid  density 

Lv.  •  coupler  tube  length 

V.-*  •  » 

X  ;  -  fluid  acceleration  in  coupler  tube 

\  •  . 

'  X  fluid  velocity  in  coupler  tube 

Vi*  ' 

C  .  flow  loss  coefficient 

The  flow  loss  coefficient  is  given  by  1 

CC  "  CD  7  CT  (6) 

who  re  1 

diffuser  coefficient  . 

C.j.  coupler  tube  friction  coefficient 


The  diffuse  r  coefficient  Cj,  includes  friction  and  entrance  and  exit 
losses  between  the  coupler  throat  and  the  corrugated  diaphragms.  It 
is  a  complex  tutu  lion  of  throat  angle,  throat  length,  and  Reynold’s 
number.  As  used  here*,,  it  includes  flow  losses  in  the  diffuser  proper 
plus  ,i  full  head  loss  assuming  sudden  expansion  at  the  end  of  the 
diffuser.  These  two  loss  components  are  plotted  in  Figure  <i  -  a  as 
'.unctions  ol  the  area  ratio  at  the  two  ends  of  the  diffuser.  The  two 
diffuser  loss  curves  are  experimental  correlations  from  Reference  i 
tor  loss  coefficients  at  the  specified  upstream  Reynold’s  numbers 
with  an  optimum  total  diffuser  angle  of  7.  5°.  These  essentially 
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R  «  1.5  x  10' 
DIFFUSER  LOSS 


bracket  the  range  ol  loss  Coefficient  values.  The  th i i*cl  curve  represents 
the  total  head  loss  for  infinite  expansion  from  the  large  end  of  the 
diffuser.  The  STEP/  coupler  design  has  an  area  ratio  of  about  6 
which  gives  a  near  optimum  total  loss  coefficient  uffO.  12  to  0.  16.  A 
conservative  value  of  0.2  is  used  in  the  simulation.  The  coupler  tube 
friction  coefficient  is  expressed  in  general  terms  by  the  Blasius 
formula  (Reference  4)  for  turbulent  flow  through  smooth  pipes  and 
represents  all  friction  losses  in  the  necked  down  region  of  the  coupler. 

It  is  expressed  as 


0.  i  164 
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whe  re 


Re  -  Reynolds  number 


I.„  coupler  tube  length 

d  coupler,  tube  .diameter 

C is  typically  0.  2  to  0.4  which  is  on  the  order  of  C^. 


2.  1.  1  Piezoelectric  Generator  Analysis 

2.  1.  1.  I  Material 

The  piezoelectric  generator  consists  essentially  of  a  stack,  of  piezo¬ 
electric  ceramic  crystals  which  are  subjected  to  compressive  stress 
by  the  engine  and  coupler.  While  several  na  tu  rally  been  r  ring  crystals 
exhibit  piezoelectric  properties,  in  the  last,  two  decades  work  has 
centered  on  developin.:.  piezoelectric  materials  with  improved  electro¬ 
mechanical  constants  and  better  temperature  and  time  (aging)  stability. 

Progress  has  <  entered  on  ferroelectric  crystals  with  a  perovskite 
structure.  This  ’structure  is  a  cubic  c  lose -packed  arrangement  of 
oxide  ions  and  the  larger  rations,  with  the' smaller  cations  occupying 
octahedral  interstices  in  an  ordered  pattern.  Useful  pe  rovsk  ite 
pe  izoelei  t  r  ic  cum.pos.it  ions  of  lead  /.irconate  titanate  tP/.T)  have  boon 
patented  in  the  United  Stales  by  the  Clevite  Corporation.  I  .end 
zirconate  titanate  compositions  have  170  to  400  C  curie  points  .{the 
temperature  at  which  materials,  permanently  lose  their  piezoelectric 
propertiesi  which  permit  high  operating  le  mpe  rat  ures.  This  eharac- 
leristii  .  together  with  lavorable  piezoelectric  elect  romeehanicnl 
(onstanls,  enaliles  power  densities  of  a  I  least  10  w/in.  ^  to  he  achieved. 

polye  ry-sla  Mine  ceramic  P/,T  material  has  an  advantage  of  more 
s{, able  polarization  than  single  crystals  and  can  be  molded  into 


desired  shapes  before  polarizing.  Table  2-1  shows  the  pertinent 
properties  of  some  important  piezoceramic  materials. 

Appendix  A  describes  the  work  performed  by  Physics  International 
working  with  P/.T-5H,  a  material  capable  of  producing  a  high  output 
power. for  a  given  stress,  and  LTZ-]  {similar  to  P/.T-4),  a  material 
whrcrh  produces  a  lower  output  power  for  the  giveifFstress,  but  which 
can  withstand  high  compressive  stresses  without  degrading  its  piezo¬ 
electric  properties. 

Appendix  C  describes  the  analytical  techniques  utilized  at  DWDL  for 
predicting  pe rfortna nee  and  correlating  theory  with  experimental 
data.  Included  are  correlations  of  the  P/.T-5H  and  LT/.-l  data  with 
theory. 

,2.  1.  3.2  Power  Output,  Losses,  and  Efficiency 

When  the  piezoelectric  crystal  is  compressed,  energy  is  stored  both 
in  mechanical  and  electrical  form,  typically  about  equal  in  magnitude. 
Electric  energy  is  stored  as  the  charge  on  an  effective  capacitance; 
the  mechanical  energy  is  stored  as  the  potential  energy  of  a  com¬ 
pressed  spring.  The  piezoelectric  crystal  may  be  connected  to  a  load 
continuously  intermittently.  In  the  latter  case,  the  material  func¬ 
tions  as  a  voltage  generator  with  a  voltage  proportional  to  the  applied 
force.  This  energy  can  then  be  discharged  periodically  into  a  load. 

For  energy  conversion,  the  fraction  of  energy  stored  as  electrical 
energy  is  all  that  is  available  to  the  load;  however,  most  of  the  stored 
mechanical  energy  is  not  lost;  it.  is  returned  to  the  engine  when  the 
compressing  force  is  relieved.  This  fraction  of  input  mechanical 
energy,  then,  does  not  enter  into  an  efficiency  calculation  for  time- 
averaged  power. 

Both  electrical  and  mechanical  losses  do  occur  anc!  must  be  accommo¬ 
dated  in  any  efficiency  calculation.  Electrical  losses  result  from  Joule 
heating  and  electrical  hysteresis  in  the  stack.  Mechanical. losses  are 
both  internal  and  interfacial  and  typically  total  a  few  percent  of  the  out¬ 
put,  with  best  performance  occurring  at  resonance.  Internal  losses 
result  from  domain  motion  within  the  ceramic  material  and  interfacial 
losses  from  relative  motion  of  the  ceramic  with  respect  to  electrodes 
and  epoxy  at  the  piezoelectric  disc  interfaces  in  the  stack.  The  elec¬ 
trical  losses  are  dependent  on  the  electric  field  on  the  disc  and  on 
stress  and  load  resistance.  Mechanical  losses  depend  on  stress  and 
piezoelectric  stack  geometry. 

Stipe  r  tn  i  posed  on  these  de  pe  inle  nr  ies  is. a  temperature  dependence. 
Tempe  rat  are  affects  all  piezoelectric  properties.  Furthermore, 
mechanical  stresses  occur  with  temperature  changes  at  the  disc  inter¬ 
faces  bee  a  ns*-  of  <  1  i  f Te  re  nt  in  1  expansion.  Increases  in  stack  tempera  - 
tore  can  result  from  internally  generated  iiea :  associated  with  these 


Loss  mechanisms  or  from  increases  in  ambient  temperature.  Prelim¬ 
inary measurements  suggest  some  reversible  degradation  in  power 
output  with  an  increase  in  temperature  above  ambient. 

For  large,  high-power  piezoelectric  stacks,  operating  temperature 
jmist  be  maintained  low  enough  that  the  maxi operating  stress 
does  not  introduce  permanent  degradation.  The  higher  the  stack 
operating  temperature,  the  lower  is  this  allowable  maximum  oper¬ 
ating  stress.  For  PZT-4  material,  the  manufacturer's  catalog 
t  Reference  5) 'indicates  this  stress  as  12,000  psi  at  25°CPbut  only 
t)000  psi  at  100'' C .  For  long-term  stable  operation,  the  allowable 
stress  may  be  only  50°o  to  60"n  of  these  values  according  to  some 
manufacturer  representatives  or  may  be  much  higher  according  to 
others.  Only  life  testing  of  stacks  under  various  stress  and  fabrica¬ 
tion  conditions  can  unequivocally  answer  the  questions  of  maximum 
allowable  stress  for  given  requirements. 

[n  the  discussion  which  follows,  the  nomenclature,  involved  in  power 
output  and  power  loss  calculation  is  introduced.  {This  enables 
discussion  of  performance  expectations,  parametric  dependencies, 
•.Assumptions;  uncertainties,  and  relative  merits  of  continuous  oper¬ 
ation  vo  rsus  pulsed  ope  ration.  '  .v  -'. 

The  electric  field  gradient  generated  in  a  piezocrystal  is  proper-  •_  ' 

tiottnl  to  the  applied  stress 

i  K  g  f.F/A)  '.1 

where  the.  constant  of  proportionality  ,  g,  is  called  tin*  piezoelectric 
Voltage  coefficient,  E  is  the  elect  ric  field,  V  the  generated  voltage, 
i  the  ceramic  thickness,  and  F/A  the  jsiressi  or  force  applied  per  unit 
area.  In  this  equation,  the  use  of  peak  force  results  in  peak  voltage, 
and  rms  force  results  in  rms  voltage:. 

To  determine  the  output  power  generated  by  a  piezoelectric  crystal,  the- 
dependence  of  lorce  with  lime  and  the  equivalent  electric  circuit  must  be 
established.  At  the  relatively  low  frequencies  t<500  Hz)  being 
co  ns  if  le  red  lor  compressing  the  piezoelectric  material,  the  material 
is  considered  electrically  as  a  voltage  generator  in  series  with  a 
resistor  and  capacitor  {  Figure  2-61.  In  Appendix  C,  a  more  complete 
equivalent  circuit  involving  series  and  shunt  resistances  is  utilized  as 
a  basis  for  detailed  analysis. 

2.  1.  ?>.  S  Continuous  Operation 

For  coin  i  minus  (  non -pul  sed )  operation,  the- optimum  load  resistance  for 
a  non  reactive -load  .equals  the  internal  impedance  of  the  stack  plus  the 
impedance  of  the  leads.  Kecnuse  t he  »  apac ilive  reactance  is  so  much 
greater  than  the  stack  and  lead  resistance  at  the  frequencies  being 
considered,  these  can  be  disregarded.  Then,  maximum  power  trans¬ 
ferred  is  given  by 


For  Sinusoidal  compression,  this  reduces  to 


where 

X  ,  l./wC'c); 

with  -  the  angular  frequency  and  C  the  pie/.oelectric  (ceramic  capac¬ 
itance  t  A  /t 

I'  \  r  ;  :■ 

.  ,  .  \ 

where  £r  is  the  relative  dielectric  constant  of  the  material  and  t  is  the 
dielectric  constant  of  free  space  (8.85  :-:  10“^  farads/meter).. 


Appendix  A  provides  a  stepwise  derivation  leading  to 


"  r  iv'1  -  »*:  VP / 


who  re 


<r  is  the  peak  stress  -  2F  /A 

Vp y  is  1 1 1«*  volume  of  pie/.oelectric. 'material  --  At 
d  is  the  pie/.oelectric  charge  coefficient  -  c  g 
For  this  case,  the  equivalent  circuit  is  shown  in  Figure  2-7. 


Localise  ho  net  energy,  is  dissipated  by  the  capacitor,  essentially  all 
power  is  applied  to  the  load,  and  overall  efficiency'  is  near  100'!1>. 
Losses  arc  calculated  using  the  electrical  dissipation  factor  (D  or 
tanoi  which  equals  the  ratio  of  effective  series  resistance  divided 
by  Xc.  Mechanical  losses  within  the  crystals  arc  generally  small 
compared  to  these  electrical  losses.  Mechanical  losses  associated 
with  stacking  several  c  ryslals  together  and  with  coupling  the  stack  to 
the  compressing  mechanism  (e.g.,  piston)  may  be  appreciable  in 
practice.  Electrical  losses  which  do  occur  result  from  lead  resistance 
and  internal  resistance  of  the  pie/.oelec rric  stack,  which  were  assume: 
in  the  first  paragraph  of  this  section  to  be  much  smaller  than  Xc  and 
consequently  much  smaller  than  Rj  also.  .Thus  the  Joule  healing  is 
much  less  in  these  (resistances  than  in  which  results  in  the  high 
efficiency.'  While  these  losses  must  eventually  be  considered,  addi¬ 
tional  work  is  required  before  a  meaningful  estimate  can  be  made.  In 


tho  literature,  mechanical  losses  are  considered  as  those  effects  that 
Mtc  rcase  .compliance  ot  the  stack.  However,  because  stack  spring- 
back  is  used,  mechanical  losses  are  only  those  frictional  or  hysteresis 
effects  that  result  in  heat.  * 


Inductance  can  be  added  to  the  load  circuit  to  reduce  total  reactance. 

W  hen  this  is  clone,  the  optimum  load  resistance  decreases,  and  maxi¬ 
mum  power  available  Irom  a  given  stack  increases  as  the  inductive 
reactance  approaches  the  capacitive  reactance.  At  optimum,  the 
inductive  reactance  ^equals  the  capacitive  reactance,  and  stack  and  lead 
resistance  are  the  only  remaining  components  of  impedance  in  the  cir¬ 
cuit  with  the  load  resistor.  For  the  same  load  resistor  (R  =  X  )  as 
gave  .maximum  output  power  without,  inductance,  adding  inductance 
in  this  formulation  permits  four  times  as  much  power  because  the 
inductance  has  effectively  cancelled  half  the  circuit  impedance, 
thereby  doubling  the  current.  Power  to  the  load  is  I2R  and  is  the  re- 
lore  increased  four  times.  The  efficiency  of  the  circuit  is  equal  to 
the  power  delivered  to  the  load  divided  by  the  power  to  the  load  plus 
power  dissipated  internally.  When  the  internal  resistance  and  its 
field  dependence,  inductor  resistance,  and  frequency  of  operation 
are  included  as  in  Appendix  C,  this  improvement  in  output:  power 
from  adding  inductance  becomes  a  factor  of  two  instead  of  four. 

Maximum  power  transfer  can  be  obtained  by  reducing  the  load  resis¬ 
tor  until  it  equals  the  sum  of  the  stack,  inductor,  and  lead  resistances7. 
The  maximunji  power  transfer  is  over  four  times  greater  than  what  was 
possible  without  the  inductance.  However,  efficiency  has  decreased 
from  nearly  100%  to,  at  best,  50%,  if  mechanical  losses  are  disre¬ 
garded.  Approximately  half  of  the  power  is  now  dissipated  in  the  stack 
and  half  is  transferred  to  the  load. 

While  this  condition  permits  the  smallest  stack  for  a  given  output 
power,  overall  efficiency  and  stack  overheating  problems  dictate  oper¬ 
ating  at  a  load  resistance  well  above  the  level  for  maximuni  power 
transfer.  As  the  various  power  loss  mechanisms  within  the  stack  are 
better  understood,  the  design  point  can  be  optimized.  Experimental 
measurements  to  date  show  that  60%  of  this  maximum  power  transfer 
can  be  obtained  with  electrical  losses. in  the  stack  reduced  by  more 
than  half.  Mathematical  and  experimental  details  of  these  results  are 
presented  in  Appendices  B  and  C. 

The  force-displacement  diagram  for  the  continuous  connection  (noil- 
pulsed  operation)  is  approximately  as  shown  in  Figure  2-8. 


2.  1.  3.4  Pulsed  Operation 

Because  energy  stored  in  a  capacitor  is  E  =  1/2  CV  stack 

,  '  '  or 

energy  becomes 

'  '  30  ... 

A 

•j  f 

% 


\\h<Tc  (.  ^  is  the  capacitance  of  the  st;irk, 


this  energy  is  available  each  time  the  pie/.oelebt  r  ic  crystal  is  com¬ 
pressed  or  relaxed,  or  twice  each  cycle.  Average  power  is  .then 


The  to  rmulat  ion  ot  this  equation  and  the  mathematics  of  pulsed  oper¬ 
ation  with  only  resist  ive  or  resist  ive  a'nd  iiuluctive  loads  are 
presented  in  Appendix  B. 


A  ve  ra  ye  po  we  r  <!  u  r  i  ng  coin  i  mums  ape  rat  ion  (without  inductance)  can  be 
compared  with  pulsed  operation  by 


) 


This  factor  ot  n->  (is  reduced  by  various -interactions  of  the  system 
operating  in  the  pulsed  mode.  For  a  given  stack  configuration,  there 
is  an  optimum  load  resistor  to  which  power  is  ..transferred.  Energy 
discharges  from  the  stack  witha  time  constant  proportional  to  the 
capacitance  of  the  stack  and  resistance  of  the  load.  By  the  time  the 
force  is  again  applied,  not  all  of  the  energy  has  been  discharged:  in 
fact,  only  a  small  fraction  may  be  discharger!  he l ore  the  stack  is 
nga  in'  rceha  rgi ng  nit  t ht*  next  cycle. 

Allot. let*  limitation  on  pulsed  'operation  is  the  extent  to  which  the  piston 
supplying  stress  to  the  stack  can  maintain  this  stress  when  the  stack 
is  being  '.iisv’ha  rged.  This  will  be  determined  quantitatively 'by  the 
interaction,  of  the  electrical  discharge  time  constant  and  the  mechan¬ 
ical  displacement  lime  constant  of  the  piston.  Analysis  must  be  done 
mi  this  time  constant  interaction  and  the*  use  of  an  inductor  ill  the  load 
circuit  to  dele  rtnine  pracl  ica  1  operating  parameters.  Such  system 
analysis  will  provide  belter  understanding  of  the  relative  merits  of 
pulsed  and  continuous  operation. 

Qualitatively  i  d  i  s  re  ga  r  !  i  ig  time  constants).,  tin*  force  displacement 
diagram  for  pulsed operation  mav  appear  as  shown  in  Figure  J-S. 

The  stack  is  compressed  from  A  to  B  on  open  circuit,  discharged 
without  relaxing  the  force  (B  to  Cl,  stress  is  release*!  on  uprn  *  ir- 
cuit  iC  to  l)i.  and  the  slack  is  again  discharged  with  opposite  polar¬ 
ity  iB  t o  A  > ,  ' 


•J2 


2.  1.  3.  i  Parametric  Dependencies 

Both  continuous  and  pulsed  power  output  depend  on  several,  classes  of 
parameters.  System  variables  include  maximum  stress  level,  stress 
wave-- form,  and  frequency.  Stack  variables  rncludf*%e  ramie  diameter 
and  thickness  and  piezoelectric  material  parameters.  Pertinent  load 
system  characteristics  are  maximum  allowable  voltage,  load  resist¬ 
ance,  and  load  inductance. 

In  general,  lor  a  given  stress  (force /area),  output  power  increases 
with  frequency  and  with  the  volume  of  piezoelectric  material. 

Output  is,  therefore,  proportional  to  the  height  of  the  stack  and 
the  squa  re'  of  stack  diameter.  Because  the  output  voltage  is  pro¬ 
portional  to  ceramic  thickness,  the  same  output  power  can  lie 
delivered  at  a  lower  voltage'  by  .using  riiany  thinner  ceramic  discs 
instead  of  fewer  thicker  discs.  This  assumes  negligible  stacking 
losses  at  the  crystal  interfaces.  To  the  extent  that  these  losses 
occur,  output  power  and  efficiency  will  be  decreased. 

Associated  with  piezoelectric  ceramic. crystals  of  decreasing  thickness 
is  the  cost  of  fabricating  additional  ceramic  discs  for  a  given  output 
power.  A  0.  020-in.  thick  disc  costs  about  25%  more  than  a  0.  040-in. 
thick  disc.  Disregarding  potential  losses  at  the  additional  disc  inter¬ 
faces^  twice  as  many  discs  are  required  to  supply  a  given  power; 
therefore,  stack  cost  is  about  2.  5  times  as  much.  Reducing  the  disc- 
thickness  to  0.003  in.  (to  supply  30  volts  rms  without  a  transformer), 
increases  the  cost  by  3  to  7  times  above  the  0.  020-in.  thick  disc. 
Further,  7  times  as  many  discs  are  required  for  a  total  cost  of  20 
to  50  times  as  much  as  the  stack  using.  0.  020-in.  thick  discs. 

Stacking  is  much  more  difficult  and  manufacturing  reliability  decreases 
as  disc  thickness  is  reduced  below  0.020-in,  However,  the  thicker 
the  disc,  the  higher  the  voltage  generated,  and  the  higher  the  heat 
produced within  a  single  disc.  These  factors  limit  maximum  accept¬ 
able  disc  thickness.  If  stack  potential  is  more  than  500  volts  rms,  addi¬ 
tional  insulation  must  be  added  to  tin'  stepdown  transformer  winding 
to 'prevent  breakdown.  This  results  in  no  significant  penalty  until 
potentials  of  a  few  thousand  volts  are  reached.  .The  0.040-in.  thick 
disc  dcs  ign  provides  less  ihan400  volts  rms  at  the  l  ransformcr  pr  ima  ry . 

A  -more;  significant  constraint  is  heat  generated  per  disc  and  the 
ability  of  the  system  to  remove  this  heat  without  the  piezoelectric 
ceramic  overheating  and  degrading.  The?,  brass  interdisc  electrical 
connectors  also  remove  heat  radially  from  the  stack.  .Because  brass 
has  a  thermal  conductivity  about  50  times-  that  of  the  piezoelectric 
ceramic  material,  the  brass  thickness  can  be  adjusted  so  that  most 
of  the  thermal  gradient  in  the  piezoelectric  disc  occurs  over  the 
short  distance  between  the  brass  and  the  niiclplane  of  the  piezo¬ 
electric  disc  ,  thereby  making  the  temperature  increase  proportional 
to  the  piezoelectric  disc  thickness. 
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lo  t  noil  it  ate  removal  ol  heat  irom  the  stack,  the  diameter  of  the 
connectors  can  be  larger  than  the  piezoelectric  ceramic  discs, 
brass  connectors  then  serve  as  fins  to  dissipate  heat.  Further 
analysis  and  design  studies  are  needed  to  determine  a  desirable 


brass 

The 

brass 


connector  configuration.  These  studies  must  colder  stack  geometry, 
pT^er 'and  efficiency,  ambient  conditions,  and  assembly  constraints. 


Other  hies  ij>h  uncertainties  can  only  be  resolved  by  testing  the  piezo¬ 
electric  stacks  under  system  operating  conditions  for  extended 
periods.  These  uncertainties  are  primarily  in  requirements  for 
stack  preparation  and  limitations  on  stack  operating  conditions.  The 
uncertainties  in  stack  preparation  include  piezoelectric  disc  surface 
finish  requirements,  epoxy  modulus,  and  joint  thickness.  Uncertain¬ 
ties  in  stack  operating  conditions  include  interaction  of  maximum 
operating  stress,  frequency,  and  temperature,  as  these  variables 
promote  and  retard  degradation. 


2.  1.3.6  Simulated  Generator 

Two  generators  were  modeled  and  used  with  the  dynamic  simulation 
program.  One  represented  continuous  power  extraction,  the  second 
pulped  power  extraction  from  the  piezoelectric  generator.  Initial 
simulation  involved  continuous  power  extraction  because  this  was  the 
operating  mode  utilized  in  the  testing  program.  Early  test  results 
we  re  disappointing  from  a  power  density  standpoint,  and  alternate 
approaches  were  evaluated.  The  most  promising  concept  was  pulsed 
powe  r  extraction,  which  could  be  simulated  by  a  hydraulic  pump. 
Evaluation  of  this  simulated  load  concept  showed  several  advantages 
over  that  for  continuous  operation,  and  it  was  subsequently 
adopted  for  the  reference  design.  Recent  studies  and  test  results 
discussed-  in  Appendices  13  and  C  suggest  that  continuous  power 
extraction  may  w<Jll  be  the  best  approach  for  the  piezoelectric  gener¬ 
ator.  Experience  with  the  system  simulation  program  has  shown 
that  the  type  of  load  has  no  significant  impact  on  system  performance 
at  the  desired  operating  point.  For  this  reason  and  to  avoid  duplica¬ 
tion  of  effort,  the  hydraulic  load  has  been  maintained  in  the  reference 
design'.'  The  paragraphs  which  follow  describe  the  two  simulated  load 
approaches. 

In  the  complete  non isothe  r nial  simulation  of  the  system,  n  force-' 
displacement,  sc  ha  racte  ristic  must  be  included  for  the  piezoelectric 
generator.  Displacement  of  the  piston  which  compresses  the  stack 
is  keyed  to  displacement  of  the  power  piston  in  the  test  setup  because 
the  coupler  has  very  little  compliance.  Hjpwever,  force  applied  to  the 
'.generator  piston  by  compressing  the  piezoelectric  stacks  depends  on 
the  electrical  Interconnection  of  the  stacks.  A  continuous -ope  ration 
generator  simulation  test  setup  which  utilizes  hydraulic  flow  dissi¬ 
pation  and  pneumatic  bounce  chamber  is  shown  in  Figure  2-9.  Pres¬ 
sure  recovery  is  obtained  at  the  gene  r.yto.t3  end  of  the  coupler  by  a 
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d  it  fuse r.  1  his  is  the  point  at  which  work  would  be  applied  to  the 
piezoelectric  stack  if  a  real  generator  were  used.  A  metal  diaphraem 
can  bo  inserted  at  this  point  if  different  fluid  properties  are  required 
m  the  coupler  and  load  regions  for  the  real  or  simulated  load.  In  the 
simulated  load,  the  fluid  necks  flown  to  a  variable ^Ffow  restriction 
which  is  characterized  by  both  linear  and  quadratic  terms  in  flow 
velocity  to  represent  any  desired  combination  of  laminar  and  turbulent 
llow.  After  passing  through  the  load  restriction,  the  load  region 
-fluid  operates  against  a  pneumatic  bounce  chamber  acting  as  a  spring. 
It  there  were  no  flow  restriction,  pressure  volume  changes 
(Figure  2- 10)  would  move  up  and  down  the  dashed  line,  with  no  work 
performed.  T  he  overall  volume  of  the  bounce  chamber  determines 
the  stiffness  and  nonlinearity  of  this  pneumatic  spring  characteristic. 

\\  ith  a  flow  restriction^  the  generator  side  coupler  pressure  is 
higher  as  liquid  moves  to  the  right  and  compresses  the  bounce 
chamber.  Conversely,  this  pressure  is  lower  when  liquid  flow 
reverses.  1  he  resulting  P - V  diagram  is  shown  by  the  solid  curve  in 
Figure  2-  10.  Work  at  the  flow  restriction  is  converted  to  heat  and 
dissipated  by  the  cooling  coils. 


The  pulsed -ope ration  generator  simulation  test  setup  selected  for  the 
proof-of-priiieiple  system  is  shown  in  Figure  2-11.  This  is  a  very 
close  analog  to  a  piezoelectric  generator  optimized  for  pulsed  power 
extraction  (Appendix  B).  The  coupler  works  into  the  load  diaphragm 
from  the  diffuser  on  the  generator  si.de  of  the  coupler.  The  load 
diaphragm  separates  mercury  in  the  coupler  from  hydraulic  fluid  in 
the  compliance  chamber,  but  docs  not  sustain  any  pressure  difference. 
With  the  diaphragm  in  its  extreme  left  position  and  the  compliance 
chamber  at  reservoir  pressure  (FigutV  2-12),  as  the  mercury 
column  and  diaphragm  move  to  the"' right,  hydraulic  fluid  in  the  com¬ 
pliance  chamber  is  pressurized  to  Pjj  (A  to  15).  Sizing  of  the  compli¬ 
ance  chamber  and  compressibility  ol  the  hydraulic  fluid 'deter  mine 
the  slope  of  A -15.  At  Point  B,  the  high  pressure  check  valve  opens. 

As  the  load  diaphragm  continues  to  the  right,  hydraulic  fluid  is 
pumped  into  the  high  pressure  accumulator  at  constant  pressure. 

When  the  diaphragm  reverses  direction,  pressure  drops  (C-I)),  and 
the  low  pressure  check  valve  opens.  Fluid  is  returned  to  the  compli¬ 
ance  chamber  at  reservoir  pressure,,  p  ,  ns  the  diaphragm  continues 
its  leftward  motion.  At  Point  A,  the  . cycle  repents. 

The  .accumulators  allow  hydraulic  fluid'  to  bleed  through  the  throttle 
valve  at  a  constant  rate,  even  though  fluid  is  delivered  in  pulses.  A 
throttle  valve  simulates  work  (Figure  2-1  1)  by  dissipating  energy 
stored  by  the  high  pressure  hydraulic  fluid.  Heat  generated  by  flow- 
loss  in.  l-he  throttle  valve  is  dissipated  by  cooling  coils.  This  system 
has  'he  advantage  of  providing  a  ready  measurement  of  simulated 
power  delivered  by  measuring  flow  through  the  simulated  load  and 
pressure  drop  across  it.  It  can  also  be  readily  converted  to  produce 
real  work  1  > v  replacing  the  throttle  valve  with  a  hydraulic  motor. 


vj 


35 


Figure  2-10.  P  V  Characteristics  of  Continuous  Power  Load  Simutator 


Performance  calculations  with  the  two  simulation  models  show  vir¬ 
tually  identical  system  performance  at  optimum,  as  a  result  of  the 
large  coupler  inertia.  The  shape  of  the  frequency  response  curve 
at  off-optimum  conditions  differs  in  the  two  cases,  but  this  does  not 
affect  the  reference  design  system  performance^ 

_ - .  - '  —  -  ■ 

2.2  PERFORMANCE  ANALYSIS 

Using  the  analytical  methods  described  in  Section  2.  1,  performance 
of  the  proof-of-principle  device  and  a  suitable  piezoelectric  gener¬ 
ator  was  computed. 

2.2.1  P roofrof-Pr inciple  Test  Device 

i  Although  the  isothermal  analysis  was  used  to  select  the  first  refer¬ 
ence  design,  real  understanding  of  the  problems  came  with  the  use 
j  of  the  nonisothermal  simulation..  Effects  of  nonlinear  springs  and 
i  the  effect  of;  displacer-power  piston  interaction  were  observed  early. 

The  analysis  which  led  to  selection  of  the  reference  design  is  presented 
in  this  section.  The  design  was  selected  using  a  very  high  moment  of 
inertia  for  the  displacer  flywheel.  This  simulates  a  driven  displacer 
which  makes  it  much  easier  to  achieve  stable  and  consistent  data  for 
parametric  studies.  It  was  subsequently  shown,  as  part  of  the  startup 
and  transient  load  change  analysis,  that  a  pneumatically  sell -driven 
free-running  displacer  is  possible  with  a  realistic  flywheel  inertia. 

It  was  also  found  that  startup  can  be  accomplished  smoothly  by  crank¬ 
ing  the  displacer  with  an  electric  motor..  , 

2.2.  1.  1  Effect  of  Nonlinear  Springs 

Some  unexpected  and  apparently  nnpmnlbus  results  emerged  during 
development  c»f  the  dynamic  simulation.  For  example,  coupler 
oscillation  sometimes  occurred  at  a  frequency  different  from  that  of 
the  displacer,,  typically  in  integral  ratios,  such  as  five  power  piston 
strokes  for  four  displacer  strokes.  Bi-Stable  modes  were  also 
observed,  e.  g.  ,  two  cases  with  identical  parameters  and  slightly 
different  starting  transients  stabilized  with  vastly  different  operating 
conditions';  These  phenomena  were  rat  first  attributed  to  the  usual 
anomalies  which  appear  in  a  complex  new  computer  program.  Eval¬ 
uation  of  simplified  nonlinear  spring  mass  analogies  to  the  STEP/ 
system  showed,  however,  that  these  were  real  effects  rather  than 
computer  idiosyncrasies... 

Thd  STEP/,  system  is  highly  resonant  and  has  relatively  little- damping 
This  results  from  the  large  amount  of  inei  'ial  energy  storage  required 
by  the  coupler,  as  compared  with  the  work  output  of  the  system.  As  a 
result,  system  performance  is  quite  sensitive  to  operating  conditions 


»im;I  ties  ign  parameters.  1  ho  coupler  is  analogous,  to  a  driven  spring 
mass  system,  with  the  oscillating  liquid  column  as  the  mass  attached 
to  the  hydraulic  loaf!  springon  one  end  and  the  pneumatic  engine  spring 
on  the  other.  The  problem  is  complicated  by  the  fact^that  both  springs 
and  tyftgtuo- drive  force  are  nonlinear.  These  nohlinelfrities  are  the 
direct  cause  of  the  bi-stable  operating  inodes  (Reference  6).  Simi¬ 
larly,  the  differing  displacer  and  coupler  frequencies  can  be  caused 
by  either  the  nonlinearities  or  transient  effects.  These  effects  were 
greatly  reduced  in  the  final  reference  design  as  a  result  of  reduced 
compression  ratio. 

2.2.  1.2  Effect  of  Displacer  and  Power  Piston  Interaction 

The  sensitivities  and  dynamic  interactions  between  the  displacer  and 
power  piston  caused  enough  problem  in  the  simulation  runs  that  an 
early  decision  was  made  to  add  an  auxiliary  electric  motor  drive  for 
the  displacer.  This  approach  essentially  removes  one  degree  of  free¬ 
dom  from  the  system  and  greatly  simplifies  control.  The  final 
reference  design  exhibited  much  less  operating  sensitivity  than  the 
initial  design.  This  resulted  primarily  from  reducing  the  compres¬ 
sion  ratio  of  the  system,  which  effectively  increases  the  ratio  of  work 
removed  from  the  System  to  energy  stored  in  the  coupler.  The  system, 
nevertheless,  remains  highly  resonant.  « 

2.  2.  1.3  Compression  Ratio  Effects 

The  classic  solution  by  Schmidt  as  refined  by  (ioransori,  et  nl. 
(Reference  2)  is  based  on  the  assumption  of  isothermal  conditions  in 
the  hot  and  cold  gas  regions.  In  the- isothermal  program,  this  assump¬ 
tion  is  carried  through  into  thermal  input  computations.  The  primary 
factors  which  contribute  to  making  this  assumption  unrealistic  for  the 
STEP/,  engine  are  pressure  changes  over  one  cycle  and  operating 
frequency.  The  combination  of  large  pressure  changes  over  one  cycle 
with  high  cyclic  rates  results  in  rapid  pressure  changes  with  time. 

This  tends  to  cause  pressure  changes  to  be  more  nearly  adiabatic  thru 
isothermal,  with  large  resultant  temperature  changes  or  nonisotherr. i 
effects  over  one  cycle.  The  relative-  problem  between  engines  for  the 
artificial  heart  and  the  STEP/,  application  is  that  the  operating  pres¬ 
sure  ratio  is  increased  from  1.2  to  1.9  and  the  frequency  from  20  'o 
60  Hz.  The  problem  is  further  augmented  by  an  increased  heat  trans¬ 
fer  requirement  per  unit  heat  exchanger  area  as  discussed  in 
Appendix  D. 

This  problem  was  recognized  early,  but  the  isothermal  program  was 
the  oriy  means  available  to  establish  an  initial  reference  design  with 
which  to  begin  the  mechanical  design  effort.  Early  results  with  the 
nonisothc mini  program  led  to  the  conclusion  that  the  compression 
ratio  selected  for  the  initial  reference  design  .was  too  high.  The 


stroke*!  volume  tradeoff  studies,  with  the  nonisothermal  simulation, 
as  t  t.si  iibed  in  Section  2.  2.  1.  4,  were  used  to  select  a  near-optimum 
compression  ratio  based  on  solid  analytical  foundation. 

Compression  ratio  is  conventionally,  defined  as^Se  ratio  of  working 
gas  volume  at  the  beginning  of  the  compression  stroke  to  that  at  the 
end.  These  points  are  difficult  to  determine  exactly  in  the  STEP/, 
engine  but  are  very  closely  approximated  by  the  ratio  of  maximum 
to  minimum  working  gas  volume.  Based  on  this  definition,  the 
compression  ratio  of  the  engine  was  reduced  from  3.6  to  1.9  when 
the  nonisothermal  simulation  was  used  to  establish  a  new  reference 
design.  The  resultant  increase  in  displacer  stroke  caused  perform¬ 
ance  red  action  because  of  increased  mass  flow  through  the  regener¬ 
ator  with  attendant  higher  main  reheat  loss. 

2.2.  1.4  Design  Trades  and  Predicted  Performance  of 
Reference  Design 

;  '  f  ■  ."  .  .  • 

Two  reference  designs  were  used  during  this  effort.  The  initial 
reference  design  was  established  prior  to  the  start  of  contract 
work  using  the  isothermal  program.  A  summary  of  computer 
performance,  operating  conditions,  and  geometry  is  provided  in 
[  Table  2-2  for  this  initial  design  (first  column).  Following  completion 
of.  the  nonisothermal  simulation,  it  was  used  to  establish  ?  second 
reference  design  (center  column).  The  last  column  in  Timle  2-2  shows 
i  performance  of  the  second  reference  design,;  calculated  using  the 
isothermal  program;  this  calculation  is  presented  to  contrast  the 
results  obtained  with  the  two  computer  programs. 

The  major  difference  between  the  first  and  second  designs  is  the 
compression  ratio.  The  isothermal  program,  which  was  used  to 
establish  the  initial  design,  predicts  ever  increasing  performance  with 
increasing  compression  ratio.  A  compression  ratio  of  3.  6  was 
arbitrarily  selected  for  the  initial  reference  design  as  representing  a 
practical  upper  limit.  The -compression. 'ratio  in  the  second  reference 
design  was  reduced  to  1.  9  because  of.  the  nonisothermal  heat  transfer 
effects  discussed  in  Section. 2.  1.  L  2.  This  was  accomplished  by 
reducing  coupler  swept  volume  and  increasing  displacer  stroke  length. 

In  arriving  at  the  new  reference  design,  other  geometric  parameters 
were  varied  in  a  limited  parametric  survey  using  the  simplified 
version  of  the  nonisothermal  simulation.  No  significant  changes  were 
indicated.  However,  minor  changes  weire  made  for  design  convenience 
For  example,  engine  diameter  was  reduced  from  1.  30  to  1.44  in.  so 
that  previously  designed  displacer  support  flexures  could  be  used 
directly.  In  the  following  paragraphs,  some  of  the  design  tradeoffs 
and  engine  dynamic  studies  are  discussed. 
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Table  2-2 

REFERENCE  DESIGN  PERFORMANCE  AND  SPECIFICATIONS 


Or  igma  l/lsot  her-mal 

Computed  Performance  1. eve  Is  •  program 

New/Non-lso- 
thcrmal  simul. 

New /isothermal 
Program 

Power  Output  (w) 

Gross  engine -output 

SIS’  . — 

'*?538 

477 

Windage 

IK 

20 

20 

Drive  power 

•  t 

31 

Net  engine  output 

497  ■ 

540 

44V 

Coupler  loss 

Power  delivered  to  simulated  or 

11(1^)* 

22 

22 

real  load 

F.stimatod  mechanical  toss  in  Pa 

48t.  133«)* 

518 

TIT 

stack  driver  (90"„  efficient ) 
Flectrical  and  mechanical  losses 

49  (34)* 

4 

52 

43 

in  P/.  slack 

44  (  101* 

47 

18 

Gross  electrical  p*»\\rr  generated 

393  (..'274)* 

4  10 

,»4.p 

Power  Input  (w) 

Second  law  input 

■  811  ' 

091 

7t.3 

Main  reheat  loss 

•37% 

1108 

hi  3 

Cylinder  wall  conduction  loss 

45 

45 

45 

Displacer  wall  conduction  loss 
Displacer  gas  conduction  and 

1.0 

10 

10 

radiation  loss 

3 

3 

i 

Insulation  loss 

20 

20 

20 

1  otal 

!2o4  . 

2177 

14  74 

Component  and  System  Efficiency  (g«.) 


hngine  onlyj 

34.3  ■ 

24.  8 

30.  5 

Kngine  and  coupler 

38.  4  (26.  7)6 

23.  8 

29.  0 

Complete  system 

ill  1  (21.  7)* 

19.2 

23.  5 

Operating  Ootid  mens 

Power  piston,  phase  lag  (“I 

»() 

00  to  101 

«0 

Ope  ratine  (requeue y  (!!/.) 

oO 

t.O 

t.O 

Maximum  wclo  pressure  (psia) 

2  tO'l;- 

26 12 

227) 

Minimum  cycle  pressure  (psia) 

MV 

1  1  41 

1  J0>. 

Average  cy.  le  pressure  (psia) 

1  <S  V 

1  7  4.8 

1  758 

Maximum  load  pressure  (psia) 

-  • 

2070 

Minimum  load  pleasure  Ipsia) 

•- 

1  502 

'  - 

Over  .ill  pressure  ratio 

4 

1 .  7 

Heat  ii.ource  lempe ratu re  ( 'O 

t.sO 

« i  0 

I.-.0 

Average  hot  gas  tempera! are  (l,0 

t.  37 

5t.l 

1.44 

Heat  >inl:  temperature  ("*(’.) 

-lit  104  *  F 

40 

40 

Average  cold  gas  tempi*  rat  ire  '1  #C  1 

‘ 

70 

.  t  >  8 

Coupler  swept  \tdume  (cti  in.) 

0.  4/4 

0.214 

0.214 

(irnmrirU  para  mete  is 
(  oiipier: 

Throat  diameter  (in.) 

0.  V 

0.  ' 

«1t  s 

K: Jet  »i.ve  length  at  throat  dia  (in.  ) 

i  o. 

!  0.  2 

10.2 

K-ngir.e: 

Cylinder  diameter  (in.) 

1.4 

1.44 

1.44 

Hispl:»e.er vlind«*r  radial  clearance 

0.  Otis 

0.  00  . 

0.  00s 

(in.  1 

hispl.in**'  stroke  length  (in.  ) 

0.  Ot.O 

o;  Ouo 

o.  oqo 

«t  swept  volume  (a  n  in-).. 

0.  10). 

0.  14  7 

0.  14  7 

Hot  end  <  lea  ranee  f  in.  )  * 

n.  op) 

0.010 

0.  010 

Additional  h«*t  .dead  -volume  (cu  in.)  * 

0.  0 

o.  on 

0.014 

Cold  end  «  h'arane**  (in.  ) 

•ft.  >1*1 

0.  00  v 

ft.  ION 

.  Additional  •  old  dead  ^volume  feu  iii.  ) 

0.tll8 

o.tv>; 

ft.  ftV,* 

O*. era M  t  ylinder  length  (in.  ) 

;.  o 

l 

7.1.. 

Heater  length  (in.  V 

l.o 

1 . 4/ 

1 . 42 

‘ !N*ger«er.»tor  length  (in.) 

•1 .  »v 

i.  »'V 

Cooler  length  ( tn.  ) 

*.  \ 

1  .  *>•; 

l  »7 

Drive  piston  diameter  (m.) 

0.  0 

0..V* 

ft.  ft 

Compression  ratio 

y .  t  • 

1 .  •> 

1  •• 

Principal  vale*"  ,i*»!4tunr  *»did  plstmi  cutis idered  lor  original  deMgn. 

f \t rent|»et(ca)  values  N*  mercury  coupler  lot*  direct  temper iwon  with  ether 


•vi  se*  , 
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The  nature  of  the  nonisothermal  design  trades  and  the  relative 
performance  sensitivities  are  shown  in  Figures  2-13  through  2- 16. 
Figure  2-13  shows  performance  (engine  power  output  and  efficiency 
with  no  parasitic  losses)  as  a  function  of  displace r-cylinder  radial 
clearance.  This  parameter  h^s  a  stronger  effect  pcrfperformance 
thah-any  other.  (Optimum  performance  occurs  near'li  mils,  but 
5  mils  was  selected  as  the  design  point  on  the  basis  of  practical 
fabrication  limitations  and  maintenance  of  a  reasonable  windage  value. 


Figure  2-14  shows  the  effect  of  tradeoffs  in  the  heat  exchanger  and 
regenerator  lengths  where  the  total  cylinder  length  is  held  fixed  at 
a  practical  limit  of  7.0  in. ,  which  is  also  a  little  below  optimum. 

The  selected  design  point  sacrifices  some  efficiency  to  permit 
enough  cooler  length  for  close  coupling  of  the  power  piston  diaphragm 
unit  below  the  cold  cylinder  wall  flange. 


The  effect  of  compression  ratio,  in  terms  of  parametric  variations 
between  displacer  stroke  and  piston  pumped  volume,  is  shown  in 
Figure  2-15.  The  original  isothermal  reference  design  is  far  off- 
optimum  with  the  nonisothermal  analysis.  Some  efficiency  com¬ 
promise  was  made  in  the  new  design  to  raise  the  power  level  above 
500  w.  This  reflects  a  power  density  limitation  for  the  simple 
annular  regenerator,  at  a  reduced  compression  ratio,  which  could 
not  be  observed  with  the  isothermal  program.  This  effect  is  empha¬ 
sized  in  Figure  2-16  which  shows  performance  of  the  new  reference 
design  as  a  function  of  ^frequency.  At  20  Hz,  power  has  decreased  to 
200  w,  but  the  efficiency  has  increased  beyond  40%.  Parasitic  losses 
were  not  included  in  these  calculations,  but  their  effect  would  be  to 
decrease  efficiency  at  20  Hz  to  about  36%  with  a  sharp  drop  to  zero 
efficiency  at  zero  frequency. 

In  Figure  2-16,  it  is  assumed  that  the  charge  pressure  remains  fixed. 
This  means  that  coupler  length  must  be  increased  as  frequency  is 
reduced.  The  proof-0f-principle  system  will,  however,  be  capable 
of  demonstrating  increased  efficiency  at  reduced  power  levels  by 
operating  at  an  appropriate  combination  of  reduced  frequency  and 
charge  pressure., 


After  final  selection  of  the  basic  geometry  and  operating  conditions 
with  tiie  simple  version  of  the  nonisothermal  simulation,  the  more 
cortiplete  version  was  used  for  additional  refinements.  With  the 
driven  displacer  and  free  coupler  motion,  fine  adjustments  were 
made  with  load  pressures,  gas  charge,  and  coupler  length  to  achieve 
stable  resonant  operation  at  60  Hz  with  the  engine  diaphragm  nearly 
bottoming  out  bn  the  support  plate  for  minimal  dead  volume. 

Selected  computer  plots  in  Figures  2-17  and  2-18  show  resonant 
operation  of  the  system.  Figure  2-17  shows  hot  and  cold  end  region 
gas  temperatures- a nd  load  and  engine  pressures  as  functions  of  time, 
over  one  cycle.  The  significance  of  nonisothermal  effects  is  apparent 
even  at  the  low  compression  ratio  of  this  system.  Hot  gas  tempo  rntun 


efficiency  (%) 


REGENERATOR  LENGTH  (CM) 


HOT  GAS  temperature  r K) 
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NE  PRESSURE  (DYNE/CM^  X  lO'^ATM) 


varies  irony  60  °C  above  the  heat  source  temperature  to  226  °C 
below  it.  Similarly,  the  cold  gas  drops  iti'C  below  the  heat  sink 
temperature  and  rises  to  99 6 C  above  it.  Load  pressure  exhibits 
relatively:  steep  rise  pnd  fall  characteristics  with  n^&st  of  the  time 
speTTT with  intake  and  discharge.  The  engine  pressure  trace  is 
considerably  broader  on  the  low  pressure  portion  of  the  cycle,  with 
the  arithmetic  average  pressure  of  1876  psia  at  138  psi  above  the 
time  average  pressure  of  1738  psia.  One  pressure-displacement  and 
three  pressure -volume  diagrams  are  shown  in  Figure  2-18.  The 
area  enclosed  by  the  hot  and  cold  region  diagrams  respectively 
rep1-  sents  second  law  heat  addition  to  and  rejection  from  the  system. 
The  area  of  engine  pressure  versus  total  working  gas  volume  diagram, 
gives  the  cyclic  work  output  of  the  engine.  The  area  of  the  load  dia¬ 
gram,  multiplied  by  the  effective  diaphragm  area  of  45.7  cm^,  gives 
the  work  delivered  to  the  load. 

The  displacer  was  driven  at  a  series  of  frequencies  around  resonance 
to  develop  a  frequency  response  curve  for  the  system.  Results  are 
plotted  in  Figure  2-19-  Power  peaks  with  a  fairly  flat  characteristic 
right  at  60  Hz.  Power  drops  off  steeply  below  the  resonance  fre¬ 
quency,  but  comparatively  slowly  with  increasing  frequency.  This 
characteristic  is  the  most  pronounced  effect  of  using  the  pulse 
power  generator  simulation  as  opposed  to  continuous -connection 
piezoelectric  generator  simulation.  With  the  latter,  the  falloff;  from 

resonance  is  sharper  on  the  high  frequency  side  and  shallower  as 
lower  frequencies  are  approached.  This  results  from  the  differing 
characteristics  of  the  nonlinearities.  The  resonant  point  shows  vir¬ 
tually  the  same  Der  forma  nee,  however.  Efficiency  shows  a  much 
sh,!arper  peak  than  power  dens ity,  which  reflects  greater  efficiencies 
possible  at  lower  power  densities.  The  phase  lag  of  the  power  piston 
behind  the  displacer  is  near  zero  at  frequencies  much  below  resonance. 
Phase  lag  increases  sharply  toward  the  resonance  value  of  about  90°, 
then  asymptotically  approaches  180°  with  increasing  frequency/. 

2.2.  1.5  Free  Running  Displacer 

All  computer  runs  discussed  utilized  an  extremely  large  moment  of 
inertia  for  the  displacer  flywheel  to  effectively  drive  the  displacer 
at  a  fixed  frequency  equal  to  its  initial  value.  The  effect  of  a  realistic 
flywheel  inertia  is  consequently  of  interest.  An  inertia  of  450  g-cm^* 
equal  to  that  of  the  system  design,  was  used,  and  the  system  was 
started  from  equilibrium  resonant  conditions.  The  result  is  shown 
in  Figure  2-20,  which  plots  angular  velocity  of  the  flywheel  as  a 
function  of  time.  After  a  slowing  transient  which  lasted  for  about 
5  cycles,  the  system  was  run  for  15  additional  cycles  (only  5  shown) 
aw!  stabilized  nicely.  Because  ol  non-uniform  d  rive  forces,  angular 
velocity  varies  from  about  3  50  to  407  radians  per  second  over  one  cycle. 
This  represents  fully  self  driven  operation  with  the  drive  motor  coast¬ 
ing  awl  all  windage  forces  overcome  with  the  pneumatic  drive  power, 
if  bearing  and  seal  friction  on  the  displacer  drive  shall  are  disregarded 
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Figure  2-19.  New  Reference  Design  Performance  Characteristics  Versus  Frequency 


Figure  2-20.  Transient  and  Stable  Operation  of  Fully  Self  Driven  System 


;2.  2.  1.  6  Startup,  Shutdown,  and  Transient  Load  Change 
j  V'.  Characteristics  ‘ 

Several  studies  were  made 'with  the  nonisothermal  simulation  to 
investigate  startup, ..shutdown,  and  transient  -loac^change  dynamics. 

Startup  and  Shutdown  Characteristics 

Two  types  of  startup  were  investigated.  These  include  (1)  an  increase 
m  displacer  frequency  from  0  to  slightly  above  the  design  value  with 
the  heater  and  cooler  temperatures  held  constant,  and  (2)  an  increase 
in  heater  temperature  from  cooler  temperature  (40°C)  to  the  refer¬ 
ence  design  heater  temperature  with  the  cooler  temperature  and 
displacer  frequency  held  constant. 

The  increasing  frequency  startup  simulates  a  condition  where  heat  has 
been  applied  and  operating  temperature,  has  been  achieved;  at  that  time, 
displacer  motion  is  started.  During  the  "initial  heatup  before  the 
displacer  is  started,  pressure  in  the  engine  increases;  this  moves 
the  engine  power  diaphragm  from  its  initial  position  to  a  new  equilib¬ 
rium  position  ivhich  depends  on  the  displace r  position  (i.  e.  the  distribu¬ 
tion  of  gas  between  the  hot  and  cold  regions).  These  are  the  starting 
positions  for  displacer  and  engine  diaphragm  when  displacer  motion  is 
started.  During  heatup  or  cooldown,  a  condition  is  never  reached  in 
which  either  diaphragm  bottoms  out  against  its  backup  plate  as  a  result 
of  static  pressures. 

In  the  increasing  ifrequgncy  simulation,  a  forcing  function  was  applied  to 
the  displacer  to  increase  its  angular  velocity  from  0  to  400  rad/sec 
(equivalent  to  63.  66  Hz)  in  three  seconds  of  real  time.  Because 
effects  near  each  end  ol  the  speed  range  .are  of  primary  interest,  the 
forcing  function  was 

0  200  [j  -  cos(  I.  04  721)] 

where  0  is  the  displacer  flywheel  angular  velocity  in  rad/sec  and  t  is 
(time  in  sec.  This  function  causes  the  displacer  to  start  with  zero 
frequence  a  nd  acele  rat  ion  and  Increase.' speed  slowly  at  first.  Maxi- 
jmum  acceleration  is  at  an  angular  velocity  of  200  rad/sec,  which  is 
'reached  i.  50  sec  after  start.  The  speed  then  continues  to  increase 
to  400  rad/sec  at  00  sec  while  acceleration  is  decreasing  to  zero. 

The  first  full  cycle  takes  0.56  sec  as  contrasted  with  the  final  ■  cycle  of 
0.0  16  sec.  A  hot  shutdown  cast*  is  simulated  by  this  satin*  plot  with 
a  reversal  of  the  time  scale  or  0  1.66  llz  at  zero  time  decreasing  to 
zero  Hz  at  A.  00  seconds. 

St  a  rf  ing  with  the  rei'e  fence  <lesign  conditions  (except  for  frequency), 
the  computet  simulation  required  15  hours  to  determine  system 
response  to  this  forcing  function.  Power  piston  (diaphragm)  position 
was  plotted  as  ;*  function  of  time  on  a  plot  about  11  ft  long,  to  achieve 
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i-ensojiabie  resolution  of  the  high  frequency  oscillations.  The  length 
ot  this  plot  makes  it  impractical  to  include  in  this  report,  because 
photo  reduction  to  a  practical  size  would  destroy  essential  detail. 

y  1 l?'Toimadt  Witih  the  connectin«  rod  force  (F_,  Equation  ; 
of  Section  2.  1.  J.  i.)  as  the  plotted  variable.  This  showedcthat 

instantaneous  forces  daring  startup  never  exceede^Hhe  design  bear- 
intfloads.  Highest  loads  we  re  experienced  at  frequencies  above 
resonance.  These  plots  are  furnished  separately  for  MERDC  review. 


One  unexpected  but  readily  explained  result  emerged  from  the  first 

run.  A  low- to -mode rate -amplitude  75-Hz  oscillation  of  the  mercury 
column,  about  a  central  point  directly  proportional  to  power  piston 
position,  started  almost  immediately.  This  oscillation  continued  until 
the  displacer  reached  about  50  Hz,  after  which  the  power  piston  syn¬ 
chronized  with  the  displacer  frequency.  In  normal  resonant  opera¬ 
tion,  the  effective  spring  constant  of  the  load  is  proportional  to  the 
slope  of  line  A-B,  Figure  2-12,  during  compression  of  the  fluid  in 
the  hydraulic  chamber.  During  expulsion  of  the  hydraulic  fluid  into 
the  high  pressure  accumulator  from  B  to  C,  the  spring  constant  of  the 
load  is  zero.  The  net  effective  value  of  the  spring  constant  is,  there¬ 
fore,  proportional  to  the  slope  of  a  line  from  A  to  C.  The  system  has 
been  sized  to^produce  60  Hz  operation. with  this  effective  spring  con¬ 
stant.  When  the  frequency  is  appreciably  below  resonance,  however, 
the  power  piston  is  not  moving  far  enough  to  expel  any  fluid,  so  the 
system  oscillates  along  a  portion  of  A-B.  This  effective  spring  con¬ 
stant,  together  with  that  of  the  engine,  is  sudh  that  the  coupler 
oscillates  at  the  observed  value  of  75  Hz.  This  oscillation  occurs 
with  a  minimal  driving  force,  because  the  only  dissipation  mechanism 
is  from  flow  losses  in  the  coupler. 


The  envelope  of  the  power  piston  oscillations  follows  displacer  motion, 
with  no  phase  lag,  to  about  25  Hz,  where  an  irregular  beating  betweeri 
the  driving  frequency  and  the  75-Hz  coupler  frequency  begins  and 
continues  to  47  Hz  where  the  frequencies  synchronize.  Starting  at 
56  Hz,  *he  power  piston  amplitude  increases  rapidly  to  its  resonance 
value  at  60  Hz,  then  dec reases  almost  imperceptibly  at  the  maximum 
frequency  of  63.  66  Hz.  All  observations  arc  consistent  with  pre¬ 
viously  observed  or  calculated  phenomena. 


The  increasing  temperature  case  simulates  a  situation  where  the  dis¬ 
placer  is  brought  to  the  driven  frequency  of  60  Hz,  and  heat  is  then 
applied.  The  simulation  starts  with  the  entire  engine  at  a  tempera¬ 
ture  of  :40°t  and  increases  heater  temperature  linearly  to  650° C  in 
thrde  seconds  of  real  time.  This  simulation  also  requires  15  hours 
of  coriipute r  time  and  results  in  a  very  lengthy  plot  of  power  piston 
(diaphragm)  position  versus  time,  which  is  also  furnished  separately. 


This  type  of  startup  yielded  a  beating  effect  attributable  to  a  starting 
transient.  Amplitude  varies  ns  well  ns  frequency  during  this  starting 
transient.  These  transient  effects  decay  gradually  and  have  almost 
completely  disappeared  after  100  cycles  (1.67  seconds  of  real  time): 
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after  this  point,  amplitude  gradually  increases  uniformly.  At  the 
end  pt  3  seconds,  heater  temperature -is  923°K  and  the- linear  temper¬ 
ature  increase  is  stopped;  the  heater  is  then  held  at"a  constant  tem¬ 
perature.  When  the  heater  temperature  increase  is  removed,  the 
average  position  of  the  power  piston  moves  toward  the  engine;  all 
variables  (ppwe r| piston  position,  pressure,  temperature,  etc.) 
rapidly  stabilize  to  the  reference  design  values.  The  work  output 
and.  efficiency  are  also  equal  to  the  reference  d.esijjef:  values.  This 
curve'also  simulates  loss  of  the  heat  source  with  the  displacer  driven 
at  60  Hz  by  reversing  direction  of  the  time  scale. 

The  .most  significant  result  of  these  startup  runs  and  the  static 
analysis  is  that  no  undesirable  effects,  such  as  a  large  amplitude 
transient,  is  expected  during  startup. 

Load  Change  Characteristics 

The  basic  load  characteristic  for  simulated  pulsed  load  operation  of  a 
piezoelectric  stack  (as  in  Appendix  B,  Case  IV)  is  in  the  form  of  a 
parallelogram  shown  in  Figure  2-12.  The  slope  of  the  sides  is  deter¬ 
mined  by  simulated  stack  compliance  and  minimum  and  maximum 
pressures;  Pmih  and  pmax>  can  Ve  sel  to  a5y  desired  value.  For.  the 
reference  design,  Pmin  =  96  x  10°  dyne/cm2  and  Pmax  =  143  x  10 
dyne/cm  . 

The  effect  of  a  load  change  was  investigated  by  step  changes  in  the  values 
of  p  .  and  P  .,  after  stabilized  reference  design  operation  had  been 
achieved,  suclTVKat  the  average  pressure  remained  the  same.  This  is 
analogous  to  changing  the  load  resistance,  or  the  voltage  level  at  which 
power  is  dumped  to  the  load  in  pulsed  operation.  ,A  summary  of  these 
load  changes  and  results  is  presented  in  Table  2-3. 

In  all  cases,  a  run  was  started  with  the,  reference  load,  the  load  was 
changed  and  then  permitted  to  run  until  stabilization  occurred  or  until 
it  became  apparent  that  stabilization  would  not  occur  for  a  long  period; 
n't  this  point,  the  load  conditions  were  changed  to  the  reference  values 
again  until  stable  operation  resulted.  Cases  1  through  A  were  run  to 
50  cycles;  the  reference  cases  were  terminated  after  10  cycles.  The 
resultant  plots  of  power  piston  position  for  these  cases  were  necessar¬ 
ily  long,  and  are  included  separately  and  appropriately  annotated. 

Cases  1  and  3  had  the  effect  of  decreasing  the  height  and  increasing  the 
width  of  i-he  parallelogram.  Therefore,  the  effective  spring  constant  of 
the  load  (proportional  to  the  slope  of  the  line  from  A  to  C  in  Figure  2-12) 
is  decreased.  These  runs  achieved  a  stable  operating  condition  even 
though  power  output  was  low.  Cases  2  and  4  had  the  elfect  ol  increasing 
the  height  and  decreasing  the  width  of  the  load  parallelogram,  with  a 
consequent  increase  in  effective  spring  constant.  These  cases  never 
stabilized  (pdwer  piston  position  did  not  follow  a  regular  sinusoidal 
motion  of  constant  amplitude  and  frequency);  the  diaphragm  hits  the  end 
stpp  al  low  velocity  in  these  eases,  but  the  overpressure  accumulators 
prevent  da  mage. 
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These  stability  phenomena  are  consistent  with  the  engine  frequency 
response  curve  (Figure  2,  19).  ;  . An  increase  in  sprin^  constant  moves 

0^  C‘irVP  ^  ’^ht,  thereby  placing  the  operating  fre- 

Y  6°  lU  m  a  region  ot  rapidly  changing  response  (the  slope  of 

the  ^wer  and  the  re  fore  the  amplitude  curve  is  ver^teep  in  approach¬ 
ing^  the  pea  k  i  rom  the  lett).  In  such  a  region  of  rapidly  changing  charac¬ 
teristics,  unstable  operation  might  be  expected,  especially  in  a  nonlinear 
system  such  as  this  simulation  represents.  A  decrease  in  spring  con¬ 
stant  places  the  operating  frequency  to  the  right  of  the  peak;  this  region 
has  much  more  gradual  frequency  response  characteristics,  and  the 
system  should  be  considerably  more  stable. 

The  power  output  and  efficiency  associated  with  anv  of  the  load  changes 
is  very  low  and  would  not  be  a  desirable  operating  point,  but  the  most" 
important  result  ot  these  runs  is  that  no  large  amplitude  transients  or 
other  potentially  destructive  effects  occurred  as  the  result  of  change 
from  load  to  no  load  or  partial  load. 

These  load-change  results  emphasize  the  highly  resonant  nature  of  the 
STEP/ .system.  In  the  moderate  load  change  cases,  stabilized  operation 
is  possible  if  the  system  is  tuned  to  the  new  operating  conditions  by 
changing  frequency  or  other  parameters.  .The  most  likely  operating 
mode,  based  on  present  limited  knowledge,  is  to  utilize  a  variable  para¬ 
sitic  load  to  mairita in  fixed  output  voltage  at  thecgenerator  terminals. 

2.  2.  4  Piezoelectric  'Generator  ,|  ■ 

An;  analytical  iormalism!  has  been  developed  which  accurately  predicts 
thd  electrical  output  characteristics  for  a  given  piezoelectric  stack.  1 
Appendix  C  describes  this  formalism  and  analyzes  data  gathered  at 
.  D\V  1)1..  using  this 'approach.  Both  analysis  and  experiment,  show  that 
an  attractive  generator  can  be  designed  using  property  values  and 
limitations  of  the  piezoelectric  material  given  by  the  suppliers. 

Tablr-  2-4  gives  the  results  of  DWD1  calculations  for  a  typical  300  w(e) 
generator  tbgothe  r  with  calculations  and  measured  pe  rformance  for  a 
smaller  stack  operating  under  similar  stress  and  load  conditions. 

These  measurements  verify  the  feasibility  of  pieizoelectric  stacks 
operating  at  a  ijiavc  r  density  of  1  5  watts  /in,  ^  of  piezoelectric  ceramic 
and  a  stack  efficiency  near  70%. 

i-  !■  .  ■ 

Currently,  the  formalism  includes  only  electrical  losses.  The  dif¬ 
ference  between  predicted  and  measured  efficiencies  is  substantial. 

Ibis  fii,  le  rence  is  attributable  to  mechanical  losses  not  accounted 
for,  and  electrical  losses  incompletely  or  inaccurately  predicted. 
Nevertheless,  a  piezoelectric  generator  with  reasonable  power 
density  and  efficiency  can  be  built,  based  on  analytical  and  experi¬ 
mental  results  reported  in  Appendix  C. 
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‘  Table  2-4 

COMPUTED  PIEZOELECTRIC  GENERATOR  CHARACTERISTICS 

Stress  limit:  1000  to  7000  psi 
Frequency:  120  Ilz  '■ 

_ . -  Material:  PXT — 1  or  eqcivaleifF 


Characteristic 

Calculated 
300  w(e) 
Generator 
Module* 

Experi¬ 

mental 

Measure¬ 

ments 

Calculations 
Corresponding 
to  Experimental 
Conditions 

Disc  diameter,  in. 

1.  00 

0.  50 

0.  50 

thickness  (in.  ) 

0.  04 

0.  10 

0.  10 

Number  of  discs  l 

160 

26 

26 

Length  of  stack  (in.  )** 

6.  4 

2.6 

2.  6 

Output  power  (w) 

76 

7.  67 

7.  67 

Output  voltage  (v) 

336 

678 

678 

.  ! ..  3 

Power  Density  (w / in .  ) 

15.  0 

15.  0 

15.  0 

Load  Resistance  (ft) 

1500 

60 

60 

Load  Inductance  (hj 

1.  55 

30 

30 

Efficiency  (%) 

j  '  .  ; .  1  • 

92 

67 

93 

*  4  modules  required 

**  Disregarding  copper  or  brass  interdisc  leads  and  any  epoxy 


2.  3  MECHANICAL  DESIGN 

The  throe  major  components  of  the  design  are  the  engine,  coupler,  and 
load.  This  concept  is  identical  to  hardware  systems  previously  produced 
at  DWDL,1  except  for  the  addition  of  the  fluid  coupler  which  replaces  the 
flywheel  and  mechanical  coupling  of  a  conventional  engine-generator. 
When  the  coupler  operates  in  the  resonant  mode  with  about  90°  phase 
lag,  the  work  diagram  is  nearly  identical  with  one  for  a  conventional 
Stirling  engine.  This  increases  power1  output  by  about  a  factor  of  seven 
compared  with  the  nen-inertial  operating  mode  of  the  artificial  heart 
engine.  Increased  efficiency  also  results  because  parasitic  and  reheat 
losses  change  little  between  the  two  operating  modes.  Compared  with 


the  heart ienginb,  the  STEP/,  engine  operates  with;  about  ]0  times  the 
charge  pressure,  3  times  the  frequency,  and  one-half  the  displacer 
stroke.  Because  these  factors  are  essentially  proportional  to  power 
output,  the  STEP/,  engine  power  level  can  be  scaled l.irom  the  5-watt 
he  a?  I  engine  as  a  watts  x  t  x  10  x  3  x  1/2  =  525  wnTts,  which  is  very 
close  to  the  actual  value. 

r  -  ■  •  V  /  .  ; 

These  differences  in  operating  conditions  have  a  significant  effect  on 
mechanical  design  stresses  and  loads.  Hardware  design  of  the  engine 
and.  simulated  load  is  based  primarily  on  extrapolating  existing  and 
proven  hardware  approaches  to  the  required  operating  conditions  and 
accomodating  system  sensitivities  introduced  by  the  coupler. 

Major  components  of  the  proof-of-principle  generator  are  attached  to 
an  I-beam  (Figure  1—6)  which  is  the  basic  support  structure  for  the 
system.  The  overall  design,  approach  emphasizes  engineering  judgment 
and  economy  in  the  use  of  proven  and  low-risk  hardware  approaches. 

2.  3-  1  Engine  (Figure  2-21) 

The  engine  includes  the  following  components. 

1.  Cylinder  .  .  .f . 

2.  Displace rj  apd  supports 

3.  Displacer  flrive  mechanism  and  housing 
A.  Thermal  insulation  housing 

.  I  u  ;  v  /  . 

T he -j drive  housing  is  the  basic  structure  to  which  major  engine  com¬ 
ponents  attach.  The  engine  cylinder  mounts  to.  the  drive  housing  with  a 
bolted  flange  connection  for  ease  of  construction,  assembly,  and  disas- 
sembly.  The  power  piston  diaphragm  is  attached  to  the  side  of  the 
housing.  !  The  d  rive  shaft  cross  head  is  supported  in  bearings  which  are 
mounted  in  the  housing.  The  insulation  housing  attaches  to  the  upper 
flange  of  the  cooler  section  of  the  engine  cylinder.  The  crankshaft  and 
electric  startup  and  drive  motor  are  assembled  into  the  spindle  housing 
which  is  bolted  to  I  he  main  housing  (Figure  2-22). 

-  ’  •  '  ;  r 

2.  3.  1.  1  Cylinder 

The  engine  cylinder  consists  of  the  heater  and  cooler  section  and  a 
central  section  lo  which  the  hea le  r  a nd‘  coole  r  section  are  connected  with 
flanged  joints.  This  I  h  roe  -piece  construction  allows  the  use  of  high- 
the  rinal-coiuluctanee  male  rials  in  the  healer  and  cooler  ends  of  the 
cylinder  and  low-conductance  material  in  tin1  central  portion  to  mini¬ 
mize  cylinder  sidewall  pa  rasilic  heal  loss. 
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1  his  approach  is  used  in  the  p^oof-of -principle-. tie's ign* to  -avoid  electric 
heatt'r  burnout  problems  and  simplify  the  design  amf  fabrication  of  the 
cold  heat  exchange  rj.  In  a  flame  heated  system,  the  entire  cylinder 
could  ho  fabricated  in  one  piece  from  a  single  high-strength,’  low- 
conductance  alloy.  '  b 

The  central  portion  of  the  engine  cylinder  is  Jneon<^62S,  selected  for 
relaTively  good machinability.  A  production  version  of  the  device  is 
likely  to  use  hone1  41  for  the  cylinder  material  because  of  its  higher 
strength  at  temperature.  The  ID  is  1.437  in.  nominal  with  a  0.]  50-in. 
wall  thickness,  and  operates  lit  a  hoop  stress  of .  14,  000  psi  at  2600  psi 
max  operating  pressure.  This  stress  is  a  factor  of  two  below  the 
Inconel  62  5  value  of  29,000  psi.  for  0.  001%  creep  in  1000  hours.  At 
the  29,000-ps.i  level,  creep  per  year  is  less  than  one  mil  in  length  and 
about  1/10  mil,  in  diameter.  The  cylinder  will  be  proof  tested  at  1.5 
times  operating  pressure  to  ensure . integrity  of  the  structure.  Side 
wall  heat  loss  is  45  watts. 

The  heater  section  is  coupled  to  the  central  portion  of  the  cylinder 
with  a  flanged  joint  using  a  Kaskell  HV3-29  V-type  high-temperature 
self -energising  seal.  This  seal,  which  has  proved  satisfactory  in 
other  DWDL  engines,  is  gold-plated  Rene'  41  operating  with  a 
0,  003-in.  preset.  The  flange  has  eighteen  0.25  in.  diameter  Wasp- 
alloy  high-temperature  bolts  preloaded  to  900  lb.  With  an  effective 
cross-sectional  area  of  0.  0355  sq  in.  for  each  threaded  bolt,  the 
900-lb  preload  yields  an  operating  stress  of  25,400  psi  which  is  well 
below  acceptable  limits.  The  maximum  operating  pressure  in  the 
engine  is  2600  psi,  \vhich  for  the  1.44  in.  diameter  engine,  results 
in  an  axial  force  of  4200  psT  This  is  transmitted  ns  about  235  lb 
force  on  each  of  the  18  bolts.  Because  this  is  well  below  the  preload 
value,  cycling  stress  docs  not  strain  the  bolts  and  no  fatigue  occurs. 

The  10.  OOO-^hr  0.  5%  creep  stress  value  for  Waspalloy  is  60,000  psi, 
which  provides  more  than  a  factor  of  two  margin  above  the  actual 
stress  levels. 

The  heater  block  is  constructed  of  chromium  copper  because  of  its 
high  heat  transfer  characteristics.  It  is  bra/.ed  to  an  Inconel  625  jac¬ 
ket  to  ensure  dimensional  stability  and  carry  the  required  loads. 

Twelve  350-vv  cartridge  heaters  are  used,  resulting  in  a  4200-w  heating 
capability,  well  abov'e  the  expected  nominal  requirement  of  2200-w.  The 
maximum  allowable  heater  element  temperature  is  J600CF,  which 
represents  1400°F  at  the  cartridge  surface.  Therefore,  the  cartridges 
are  silver  soldered  to  the  copper  block  to  minimi/.e  thermal  contact 
resistance.  With  1200°F  at  the  cylinder  surface  of  the' copper  block,  the 
maximum  heater  cartridge  surface  temperature  under  nominal  operating 
conditions  is  1240°F.  Healers  can  be  replaced  by  drilling  out  and 
re brazing,  which  is  standard  operating  nmcp'liiw  at  DWDL  for  this  type 
of  construction. 


Tho  cold-end  heat  exchanger  is  coupled  to  the  central '-cylinder  section 
through  a'  tlange  using  eleven  5/  16  -  in.  diameter  high  strength  holts  pre¬ 
stressed  to  2t>,  000  psi.j  Chromium  copper  is  again  used  for  combined 
high  strength  and  thejrrpal  conductivity  along  with  the  V-tvpe  Haskell 
seal  as  used  in  the:  ilaivge  adjacent  to  the  hot  heat  exchanger.  The  wall 
thickness  is  0.  3  12;  i»y,  yielding  a  maximum  hoop  stress  in  operation  of 
11.  iHQ-ps.i  compared  to  a  room  temperature  yield  s^ess  of  42,000  psi. 

I  he  lowe  r  flange  ot  the  cooler  section  is  bolted  to  the  drive  housing 
and  sealed  with  a  Viton  O-ring. 


Heat  is  removed  tromthe  system  using  water  circulated  through  a 
simlple  co61ingiloop  consisting  of  3/ 16  OD  x  l/B  in.  ID  copper  tubing 
coiled \a round  and  brazed  to  the  periphery  of  the  cold  heat  exchanger. 
The;  pressure  drop  through  this  single-pass  cooling  loop  at  the  design 
fW  of  one  gallon  pe r minute  is  15  psi.  The  sensible  temperature  rise 
of  the  water  passing  through  the  cooler  is  15°F.  The  film  temperature 
drop  is  43  °F  and  the  solid  material  temperature  drop  in  the  cold  heat 
exchanger  is  30°F.  A  maximum  rise  of  88°!'  above  cooling  water 
source  temperature  or  about  150°F  would  therefore  be  expected  at  the 
gas -metal  interface  in  the  cooler. 


2.3.  1.2  Displacer  and  Supports  (Figure  2 - Z 3 ) 

The  displacer  is  a  thin  shell  made  of  Inconel  625  with  a  wall  thickness 
of  0.052  in.  It  is  pressurized  internally  with  argon  to  900  psi  at  room 
temperature  which  increases  to  1900  psi  at:  the  arithmetic  mean  oper¬ 
ating  temperature  of  343 °C.  This  reduces  the  pressure  differential 
expe  rienced  by  the  cylinder  and  minimizes  the  regenerator  wall  thick¬ 
ness  required  for  the  cylinder  buckling  failure  mode.  Hoop  stress 
from  internal  pressure  is  10,  100  psi  maximum  during  operation  of  the 
engine.  The  critical  buckling  pressure  is  4430  psi.  Eight  radiation 
shields  are  mounted  internally. to  reduce  radiation  heal  loss  to  less  than 
3  watts,  fias  conduction  heat  loss  is  onl^  0.  1  watt,  so  the  sidewall  heat 
leak  of  10  watts  makes  a  total  displacer  axial  heat  leak  of  13  watts. 

The  ope  rating  clea  ranee  between  the  displacer  and  cylinder  wall  is 
0.  004  in.  which  is  maintained  by  supporting  the  displacer  at  both 
ends.  It  is  supported  at  the  hot  end  by  flexures  and  is  mounted  to 
the  drive  shaft  at  the  cold  etui  by  a  spherical  joint  which  accomodates 
initial  misalignments  and  eliminates  spurious  side  loads  at  the  flexures. 
This  support,  system  was  •■subjected  to  3.56  x  10^  cycles,  at  operating 
tempo  rutu  re  during 'life  testing  of  a  NI-ILI  engine  (Reference  8)  without 
evidence  of  degradation  under  a  microscopic  examination,  following  a 
crankshaft  bearing  failure. 

The  spherical  joint  at  the  drive  shaft  uses  a  standard  jO.  5-in.  diameter 
r  t  eel  ball;.  The  upper  seat  is  Inconel  625  impregnated  with  molybdenum 
disulfide  dry  lubricant;  the  lower  seat  is  bearing  bronze.  Maximum 
operating  bearing  stress  is  a  (negligible  120  psi.  Precision  in  this 
joint  is  expected  to  result  in  O.OOO^-in.  maximum  lost  motion.  The 
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chspl^Cor  sopaj-atvs  Irom  the  drive  shaft  at" this  joint'  so  that  the  cylinder 
and  displace  r  ca n  be.  removed  from  the  assembly  without  interfering 
with  the  drive  mechanism.  This  enhances  ease  of  checkout  and  mini- 
mt'/.es  turnaround  time  during  the  development  phase  of  the  system. 


I  he  llexures  used  to  support  the  upper  end  of  the  displace r  are  shown 
ill—! igure  2 -24.  The  deflection  mode  of  the  flexiftti  as  the  displacer 
reciprocates  is  .shown  at  A  and  B.  A  pair  of  flexures  is  used  as  shown 
at  C.  ,D  and  E  show  how  lateral  loads  are  reacted  by  the  flexure  pair. 
Each  lie xu re  resists'  lateral  components  in  one  direction  resulting  in 
a  ilexural  System  pi  high  lateral  rigidity  yet  high  axial  compliance. 

Figure  2-25  shows  a;  pair  of  flexures  installed  in  the  hot  end  of  the 
NHLIi Engine  5.  displacer.  These  flexures  were  deflected  from  the 
relaxed  position  0.  080  in.  upward  and  0.  108  in.  downward  during  life 
test  operation  of  the  NHL.I  engine.  The  STEP/.  , flexure  is  identical 
except  for  a  material  thickness  of  0.  004  in.  instead  of  0.005  in. 

(used  in  Engine'  5  because  of  stock  availability).  The  STEP/,  flexure 
is  deflected  !,0.  090  in.  upward  only,  to  utili/.e  existing  tooling  and 
simplify  displacer  hot  end  machining.  Use  of  thinner  stock  reduces 
stress  -  concentrations  in  the  , flexure  to  s a fe  levels. 

.  j  •  ■  'j  ;/  ; 

2.  3.  1 .  ;3  Displacer  Drive  Mechanism  (Figure  2-22) 

The  purpose  of  the  displacer  drive  mechanism  is  to  provide  startup 
and  operation  of  the  displacer  at  3600  rpm  against  the  forces  of 
bearing  and  rod  seal  friction,  inertia,  and  windage.  Seal  friction 
is  about  6  avails;  and  bearing  friction  2  watts,  based  on  manufacturers' 
specifications  for  the  seal  drag  and  bearing  loads  encountered  in 
normal. ope  ration.  The  windage  of  29  watts  was  determined  directly 
from  the  computer  simulation.  The  total  displacer  drive  power  dis¬ 
sipation  is,  therefore  37  watts  under  nominal  operating  conditions. 

Gf  this,  nominally  31  watts  arc  supplied  by  pneumatic  <1  rive  power 
and  6  watts  by  the  drive  motor.  At  ol'I'-optimum  operating  conditions, 
friction  power  varies  directly  with  operating  frequency,  and  windage 
varies  with  frequency  squared.  Drive  power  is  proportional  to  output 
power  arid  is  only  appreciable  near  ^resonance.  The  maximum  load  on 
the  drive  motor,  therefore,  occurs  during  startup,  shortly  before  the 
resonant  frequency  is  reached,  and  amounts  to  about  30  watts.  The 
drive  system  is  shown  in  Figure  2-26.  In  summary,  the  electric  drive 
motor  is  included  in  the  proof-ol'-principle  system  to  (1)  provide  start¬ 
ing  and  trimming  torque,  (2)  act  as  a  regulator  to  maintain  constant 
speed,  and  (3)  permit  investigations  at  off-design  conditions. 


The  basic  linkage  system  is  a  conventional  overhanging  crank  pin  driv¬ 
ing  a  connecting  rod.  The  connecting  rod  reciprocates  the  drive  shaft 
which  kinematically  acts  as  the  crosshead  in  a  conventional  system. 
The  drive  mechanism  is  si/.ed  to  withstand  all  startup  and  steady  slate 
operating  loads  a  sjdete  r  mined  by  the  startup  simulation  described  in 
.‘lec  tion  2.2.  1.6.  Bearings  have  a  90%  probability  of  surviving 
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I0.«j°  hours  of  operation  at  or  below  the  resonant  frequency. 

shows  the  Steady  State  loads  obtained  from  the  computer 
simulation  based  on  a  maximum  operating  pressure  of  2710  psia  a 

f  °n5  slu>i’"and  an  engine  speed  of  3600  rpm. 

7  lu  tlyw  heel  is  si/ed  to  allow  a  maximum  cyclic  angular  velocity 

^nation  corresponding  to  12.4  rpm  resulting  lH§m instantaneous 
unbalance  ol  drive  hnd  energy  absorption  forces. 


1  he  crank  pin  and  main  crankshaft  bearings  are  pre lubricated 
precision  ball  bearings.  The  yvrist  pin  is  a  hardened  0.  164-in. 
diameter  pin  operating  in  an  aluminum  bronze  bushing.  The  support 
beatings  lor  the  drive  shaft  are  of  sintered  oil  impregnated  bearing 
bronze.  The  critical  bearing  is  at  the  crank  pin  and  has  an  estimated 
lite  ot  )«,  500  hours  based  on  data  furnished  by  the  supplier.  Although 
the  loads  tor  this  engine  a  re  significantly  higher,  resulting  in  more 
rugged  members  and  larger  bearings,  the  basic  approach  is  identical 
to  the  NJHLI  Engine  5  which  operated  continuously  for  more  than 
5000  hours  without  .degradation  in  performance. 


The  drive  shaft  sea l  separates  the  engine  cylinder  working  gas,  which 
fluctuates  in  pressure,  iron, i  the  drive  chamber  at  constant  pressure. 
Th|is  seal  mpst  have  low  friction  and  minimum  leakage  across  a  pres¬ 
sure  difference.  o|f  several  hundred  psi.  The  high  pressure  difference 
eliminates  from  icons  ide  ration  the  low  friction  hermetic  seal  provided 
by  a  convoluted  metal  bellows.  A  labyrinth  seal  also  had  excessive 
leakage  at  design  conditions.  However^  two  seals  proposed  by  the 
piston  ring  and  seal  department  of  Hoppers  Company,  Baltimore,  Md  .  , 
were  investigated  and  arc  satisfactory  for  this  application. 

•The  Hoppers  seal  chosen  (Figure  2-28)  is  a  0.250  ID  x  0.390  in.  OD 
split  piston  ring  of  glass-filled  Teflon  with  an  0.062  in.  wall,  split 
|outer  stainless  steel  compression  ring  to  apply  pressure  for  initial 
isealing.  The  split  inner  ring  has  a  lapped  joint  to  jminimizc  leakage. 
jThe  frictional  drag  of  the  jfeal  is  about.  5.  0  lb,  and  leakage  is  calcu¬ 
lated  as.  0.  01:  5,  cu  in.  over  one-half  of  an  engine  cycle,  both  of  which 
are  sat  is  factory  values. 


A  hysteresis  Synchronous  motor  has  been  selected  for  use  with  the  dis¬ 
placer  drive  mechanism.  The  motor  is  inside  the  pressurized  drive 
housing  and  provides  added  inertia  to  promote  stability.  This  motor 
has  the  advantage  of  speed /torque  independence  and  avoids  the  brush/ 
commutator  assembly.  An  appropriate  signal  generate  i  and  power 
amplifier  will  power  this  motor  which  is  capable  of  del  leaving  56  watts, 
or  about  a  factor  of  two  higher  than  the1  expected  transient  startup 
values. 


A  do  motor  was  rejected  for  this  application  because  the  manufacturer 
could,  not  guarantee  performance  and  life  of  the  brush/eommutator 
systems  in  a  high-pressure  helium  atmosphere.  An  induction  motor 
was  rejected  because  of  uncertainties  in  the  electronic  speed  control 
requirement  in  the  resonant  range  of  tjie  STH1V  engine. 
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Figure  2-28.  Drive  Shaft  Seal 

2.  3.1.  4  Thermal  Insulation  Housing 

The  insulation  housing  is  a  4.  62 -in;  dia  aluminum  can  which  attaches 
to  the  upper  flange  of  the  cooler  block.  This  flange  also  serves  as 
a  feedthrough  plate  for  electric  power  and  instrumentation  leads. 

The  insulating  ma teir ia  1  is  Mi ri -  K  20 0  0  ma  c  hi ne d  into  segments  to 
accommodate  irregularities  in.the  iexternal  surfaces  of  the  cylinder 
and  the  power  and  instrumentation  leads.  The  insulation  housing  is 
slipped  over  the  insulation  blocks  after  these  are  set  in  place.  The 
insulation  heat  leak  is  about  20  w. 


2.  3.  2  Coupler  (Figure  2-29) 


The  doupler  function  in  the  proof-of-principle  generator  is  provided 
by  a  mercury  inertia  column  contained  by  a  metal  tube  with  diffusers 
at  each  end.  Corrugated  metal  diaphragms  separate  the  mercury 
from  the  helium  in  the  engine  and  from  the  oil  in  the  generator  system. 
Mercury  fill  and  trim  valves  are  provided  at  each  end  of  the  coupler 
along  with  overpressure  relief  accumulators.  The  coupler  is  designed 
as  a  unit  so  that  it  can  be  removed  for  filling  or  servicing.  Weight 
minimisation  was  not  a  goal  for  this  proof-of-principle  design. 


Mercury  was  selected  as  the  coupler  fluid  primarily  becaus e  of  its 
high  density  and  low  compress ibility.  High  density  minimizes  the 
size  of  the  inertia  column;  low  compressibility  minimizes  lost  motion. 


The  coupler  housing  is  made  of  300-series  austenitic  stainless  steel. 

The  10;  5-in.;  long  coupler  tube  is  standard  0.5-in.  ID  tubing  with  a 
0.  188-in.  wall  thickness.  Hoop  stresses  are  approximately  4000  psi 
at  hydraulic  pressures  of  3000  psi.  A  5.6-in.  long  diffuser  with  a 
7'.5'cie."-rec‘  half-angle  is  provided  at  each  end  of  the  coupler  tube  to 
minim  ike  flow  losses.  End  flanges  wlth-e  fght^V  375  -  in.  diameter 
high  Strength  bolts  prestressed  to  38,  000  psi  provide  attachment  to  the 
;  diaphragm  support  structure  at  each  end.  The  mercury  inventory  in  the 
coupler  is  11.5  cu  in.  "  j 

j.  ^oth  corrugated  diaphragms  are  0.  0015-in.  .  thick  AM  350  stainless  steel. 
)  These  haye  eight  convolutions  with  a  peak-to-peak  height  of  0.  024  in.  and 
]  a  pitch  of  0..  3  5  in.  The  diaphragms  are  capable  of  operating  through  a 
swept  volume  ±0.2  cu  in.  from  the  central  unstressed  position  to  their 
extreme  deflected  position  with  a  maximum  tangential  stress  of  3  1,000 
psi.;  The  design  point  swept  volume  totals  0.  2  1  cu  in.  ,  or  about  one- 
half  the  allowable  value. 

Each  diaphragm  has  a  7075-T6  aluminum  backup  plate  to  limit  its  travel.. 
At  the  engine /coupler  interface,  the  backup  plate  is  on  the  gas  side  of 
the  diaphragm.  Gas  must  be  distributed  over  the  surface  of  the  dia¬ 
phragm  while  minimizing  gas  dead  volume.  /Eight  radial  channels  are 
provided  behind  the  backup  plate,  each  connected  to  eight  0.  015-in. 
diameter  holes  which  vent  gas  to  the  diaphragm.  The  diaphragm  is 
capable  of  supporting  a  1060 -psi  pressure  differential  across  this  hole  , 
area  with  a  maximum  of  20,  000-psi  flexural  stress  when  bottomed 
against. the  backup  plate.  The:  total  dead  gas  volume  between  the  dia¬ 
phragm  and  the  engine  cold  plate  is  0.  04  cu  in.  ;  the  pressure  drop  in 
the  manifolding  is  8  psi.  The  position  of  the  diaphragm  is  sensed  by  1 
a  proximity  transducer  /which  has  a  maximum  linear  range  of  0.  13  in. 

■  compared  to  the  limiting  diaphragm  midpoint  deflection  of  0.  142  in. 
and  the  nominal  design  deflection  of  0.  0/5  in. 

At!the  load /coupler  interface,  the  backup  plate  is  on  the  oil  side  of 
the  diaphragm.  The  vent  holes  and  radial  channels  are  larger  than 
at  the  engine /coupler  Interface  because  entrained  fluid  volume  is 
unimportant,  and  pressure  drops  would  otherwise  be  too  high. 

To  protect  the  diaphragms  from  ove  rpressurization  when  against  the 
backup  plates,  an  overpressure  relief  accumulator  is  provided  at  each 
end  of  the  mercury  column.  These  are  piston-type  accumulators  with  • 
the  pistons  springs  loaded  to  accept  fluid  at  a  preset  value  to  prevent 
.  •  damage  to  the  diaphragm.  Fill,  bleed,  and  trim  valves  are  located 
at  each  end  of  the  coupler  for  charging  and  providing  the  proper 
amount  of  me  rcury. 

'  j  ■  ... 

2.3.3  Generator  Simulator  i 

|  The  function  of  this  component  in  the  proof-of- principle  device  is  pro- 
1  vided  by  a ‘hydraulic  circuit  designed  to  simulate  characteristics  of  a 
piezoelectric  generator.  Primary  elements  of  the  load  system  as  ; 


designed  but  not  fabricated  under  the  present  contract  are  (1)  inlet  and 
.  outlet  power  check  valves,  (2)  a  liquid  compliance  chamber,  (3)  a  low 
pressure  reservoir,  (4)  a  high-pressure  accumulator,  and  (5)  an 
adjustable  main  load  valve  (Figure  1-3).  In  addilfibn,  appropriate 
s"Eutoff  valves  are  provided.  The  load  system  will  be  breadboarded 
using  standard  cornpionents  wherever  possible. 

Power  valves  are  standard  poppet  type  check  valves.  The  pressure 
drop  across  each  valve  at  maximum  flow  is  15  psi;  total  flow  loss  for 
both  valves  if  1%  of  generator  power  output. 

The  accumulator  and  reservoir  are.  standard  10  cu  in.  bladder -type 
accumulators.  Cyclic  pressure  variation  with  a  10  cu  in.  gas  volume 
is  limited  to  2%,  i  .  :  '  ' 

Design  flow  through  the  primary  Load  valve  is  approximately  6  cu  in.  / 
sec,  corresponding  to  a  power  dissipation  of  500  w.  A  0.  25-in,  needle 
valve  with  a  maximum  flow  capacity  of  27  cu  in.  per  sec  has  been 
selected.  The  outlet  Line  from  the  valve  passes  through  a  water-cooled 
heat  exchanger  to  remove  heat  generated  by  the!  power  dissipation  (flow 
loss)  in  the  valve.  All  components  will  be  connected  with  0,5-in.  dia¬ 
meter  lines  to  minimize  flow  losses  elsewhere  in  the  setup. 

2.4  PIEZOELECTRIC  GENERATOR  LIFE  TESTS 

This  section  describes  the  Phase  II  effort  which  consisted  of  life  testing 
four  piezoelectric  generator  stacks  under  varied  environmental  condi¬ 
tions.  Test  results  from  a  fifth  stack,  independently  developed  as  a 
prototype  unit  between  Phase  I  and  Phase  II,  are  also  included.  Widely 
divergent  opinions  of  piezoelectric  materials; experts  and  the  absence 
of  published  test  data  concerning  long-term  stability  of  a  piezoelectric 
generator  operating  at  high  stress  levels  made  life  testing  a  key  element 
in  establishing  feasibility  of  the  STEPZ  concept.  The  Phase  II  effort 
was  accordingly  directed  toward  accumulation  of  such  test  data  for  up 
to  2000  hours.  The  life  test  fixtures  and  associated  hardware  are  dis¬ 
cussed,  followed  by  a~de s cription  of  the  techniques  utilized  in  processing 
the  0.  047) -in.  thick  by  0.  50-in.  diameter  HDT-31  (PZT-4  equivalent) 
piezoelectric  discs  and  interleaving  copper  contacts.  Test  results  in¬ 
clude  both  the  parametric  tests  performed  on  the  MTS  system  and  actual 
life  tests. 

i  I  '  ’  '  '  ‘ 

2.  4.  1  Life  Tost  Approach 

2.  4.  1.  1  Evaluation  of  Alternate  Concepts 

Three  basic,  approaches  for  executing  the  desired  long-term  testing 
were?  independently  evaluated  in  the  interim  between  the  conclusion 
ofi  Phase  I  and  the  initiation  of  Phase  II.  These  included  an  electro- 
hydraulic  system  such  as  used  for  the  Phase  I  generator  tests  by  DWDL, 
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variations  of  a  mechanical  actuation  device,  and  using  a  piezoelectric 
motor  or  driver  stack  to  stress  a  similar,  generator  stack.  A  fourth 
approach  utilizing  a  mechanically  driven  rotary  hydraulic  valve  as  used 
by  Physics  International  and  described  in  Appendix  A  was  not  seriously 
considered.  Drawbacks  with  the  latter  approach  were  undesirable 
pulsating  stress  application  to  the  piezoelectric  stack,  high  friction  in 
the  rotary  valve  with  attendant  heat  generation  and  high  pressure  start¬ 
up  problems,  system  leaks,  and  excessive  but  undetermined  seal 
friction  on  the  stack  driver  piston.  While  none  of  these  appeared  in¬ 
surmountable,  the  required  engineering  design  effort  and  questionable 
likelihood  of  success  led  to  the  decision  not  to  evaluate  this  approach. 


Two  possible  approaches  utilizing  the  electrohydraulic  loading  technique 
used  in  Phase  I  were  seriously  considered.  High  costs  together  with 
basic  limitations  on  the  number  of  stacks  which  could  be  tested  and 
the  attainable  operating  temperatures  led  to  abandonment  of  this  approach 

The  me chanical  actuation  device  which  was  carried  through  layout 
d;esign  utilized  an  oscillating  lever  actuated  by  an  electric  motor  driven 
cam.  Both  ne'edle  bearings  and  tension  flexures  were  considered  for 
the  fulcrum  support  of  the  driving  lever,  and  alternatives  for  compres¬ 
sion  loading  of  the  piezoelectric  stack  were  a  spherical-ended  pushrod 
and  a  compression  flexure.  Factors  which  influenced  the  decision  not. 
to  pursue  this  approach  were  (1)  the  relatively  high  cost,  to  complete 
the  detail  design  and  fabricate  four  test  units  (2)  the  requirement  for 
regular  maintenance  and  compensation  for  wear  associated  with  some 
elements,  and  (3)  the  difficulties  associated  with  ensuring  that  the 
piezoelectric  stack  is  stressed  axially  with  no  lateral  components  of 
applied  force. 


Within  the  constraints  of  the  Phase  II  program,  the  only  practical 
approach  for  long-term  stressing  of  piezoelectric  generator  stacks 
appeared  to  be  with  a  similar  coaxial  piezoelectric  driver  stack  con¬ 
strained  by  a  rigid  load  frame.  This  approach  was  considered  during 
Phase  I,  but  rejected  on  the  basis  that  the  driver  stack  must  have  a 
much  larger  volume  than  the  generator  stack.  This  can,  however,  be 
used  to  advantage  where  a  small  generator  stack  is  acceptable.  The 


driver  stack  experiences  the  same  stress  cycle  as  the  generator,  but 
is  subjected  to  somewhat  higher  electric  fields  and,  therefore,  more 
stringent  operating  conditions.  Conservatively,  such  a  driver-genera¬ 


tor  stack  which  survives  a  given  life  test  has  demonstrated  survival 
of  a  generator  stack  the;  size  of  the  full  stack. 


A  crude  mockup  of  a  piezoelectric  driver-generator  stack  was  estab¬ 
lished  by  stacking  24  piezoelectric  discs  0.  5-in.  diameter  by  0.  1 -in. 
thick  with  a  load  cell  and  proximitor  in  a  1 2 -1  on  hyd rauli c  press.  Up 
to  215  volts  ac  waf  applied  to  12  or  23  discs  in  parallel,  under  various 
preloads,  with  voltage  output  measured  for  the  remaining  discs.  Sub¬ 
stantia!  displacements  but  very  low  forces  arid  generated  voltages  were 
measured.  These  results  were  attributed  to  compliance  in  the  load 
cell  and  frame  members,  which  was  corroborated  by  analysis  of  com¬ 
ponent  compliances.  This  experiment  established  the  necessity  of  a 
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very-rigid  Load  frame  with  no  high-compliance  elements  in  the  system. 
Such  a  load  frame  was  designed  and  fabricated  as  an  IRAD  prototype 
to  verify  performance  prior  to  initiation  of  the  Phase  II  effort. 

Figures  2-30  and  2-31  show  the  prototype  test  fixture,  both  disassembled 
and  assembled.  Components  include  a  Kis tier  Mg£el  901A  low  compli¬ 
ance  load  cell  for  measuring  instantaneous  stack  stress,  a  fixed -pos ition 
reference,  plate,  the  short  generator  stack,  a  holder  for  two  Bently 
Nevada  Model  324L6  proximitor  position  sensors  which  measure  gen¬ 
erator,  stack  displacement  with  respect  to  the  position  reference  plate, 
and  the  piezoelectric  driver  stack  enclosed  by  a  slotted  Pyrex  sleeve. 

The  quartz  load  cell  operates  on  the  piezoelectric  principle,  which  is 
the  only  available  type  with  sufficiently  small  compliance.  The  dis¬ 
advantage  w'ith  using  this  type  of  load  cell  is  that  output  is  not  referenced 
to  an  absolute  base,  and  therefore,  long-term  drift  must  be  accommo¬ 
dated  by  periodically  reducing  the  preload  to  zero  to  provide  a  tempo¬ 
rary  known  reference  point. 

2.4.  1.2  Independent  Compliance  Reduction  Efforts 

High  compliance  that  is  not  dissipative  is  not  a  significant  problem  for 
a  piezoelectric  generator  and  may  even  be  beneficial  in  reducing  the 
required  swept  volume  transformation  between  engine  and  generator. 
When  using  the  stack  as  a  driver,  however,  excess  compliance  results 
in  more  displacement  swung  with  less  stress  fluctuations  developed. 
Substantial  independent  compliance  reduction  efforts  between  Phase  1 
and  Phase  II  met  with  only  limited  success,  and  Phase  II  was  initiated 
on  the  basis  of  accepting  reduced  stress  fluctuation,  but  maintaining 
planned  peak  stresses. 

First  assembly  of  the.  IRAD  prototype  stack  gave  results  which  still 
indicated ;  substantial  compliance,  and  high-voltage  dc  excitation  was 
applied  to  the  stack.  Electronic  dial  indicators  capable  of  measuring 
in  microinches  were  used  to  measure  relative  displacement  at  various 
locations.  These  measurements  showed  that  structural  members 
exhibited  compliances  very  close  to  the  predicted  values,  and  that  excess 
compliance  (or  low  compressive  modulus)  was  restricted  to  the  stack 
itself.. 

A  series  of  stress -strain  measurements  was  made  for  a  variety  of  stack 
assembly  procedures  including  brass  and  copper  electrical  contacts  with 
varying  degrees  of  hardness.  Best  results  were  obtained  using  indium 
solder  which  was  melted  while  maintaining  axial  compression  on  the 
stack.  Practical  problems  included  voids  in  the  interfaces  and  indium 
squeezed  out  on  the  sides.  The  most  reasonable  compromise  was 
judged  to  be  the  use  of  fully  annealed  copper  contacts  with  no  bonding 
agent.  This  technique  produced  moduli’!  of  around  3.  3  x  10^  psi  com¬ 
pared  with  5.  3  x  10?  psi  for  the  indium  bonded  stacks  and  9.4  x  10“ 
psi  for  the  ideal  piezoelectric  material  for  the  short-circuited  condi¬ 
tion,  and  s Lightly  higher  than  this  for  practical  operating  conditions. 

It  increases  to  L  B .  4  x  If/  psi  under  open-circuit  conditions.  A  modulus 
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oi  at  least  7  x  106  psi  under  normal  operation  would  have  been  required 
to  achieve  a  full  1000  to  0000  psi  stress  swing  with  the  55  volts /mil 
excitation  voltage  used  for  life  testing. 


27*4.  1.  3  Test  Fixture  A  ssemblies 

The  Life  test  fixtures  consist  of  a  rigid  load  frame  which  encloses  the 
piezoelectric  stack  with  a.  means  for  applying  uniform  torque -free  axial 
preload  to  the  stack.  Other  elements  mechanically  in  series  with  the 
stack  include  a  quartz  load  cell,  centering  pins,  and  two  proximitors 
which  measure  axial  displacement  across  the  generator  stack.  These 
elements  for  the  IRAD  prototype  load  frame  (Figure  2-3(0  and  included 
in  a  typical  contract  assembly  (Figure  2-32).  The  IRAD  prototype  load 
frame  was  engineered  during  the  interim  between  Phase  1  and  Phase  II, 
and  was  modified  slightly  during  Phase  II  to  incorporate  improvements 
based  on  initial  operating  experience.  The  primary  change  was  an 
increase  in  the  load  frame  opening  for  betLer  assembly  accessibility 
with  only  a  slight  increase  in  compliance.  A  single  large  diameter 
thread  was  also  used  for  the  preload  adjusting  screw  to  simplify  fabri¬ 
cation  and  assembly  and  reduce  galling. 


Assembly  of  the  life  test  fixtures  is  indicated  in  Figure  2-33.  The 
lower  jam  nut  retains  the  locating  plug  which  is  held  concentric  with 
the  load  frame  on  its  outer  diameter  and  locks  the  bottom  load  pad  into 
position  with!  a  pin.  The  load  pad  is  followed  by  a  0.005-in.  thick  mica 
disc  for  electrical  isolation,  a  fully  annealed  E TP  copper  compliance 
pad,  land  the  piezoelectric  driver  stack  consisting  of  seventy  0.  5  in. 
diameter  by  0.040-in.  thick  HDT-31  (P/.T-4  equivalent)  discs  with 
interleaved  copper  contacts  0.  002  in.  thick.  The  assembly  continues 
with  a  compliance  pad,  mica  disc  load  pad,  and  proximitor  mounting 
plate  which  uses  a  locating  pin  to  mate  with  adjacent  load  pads. 

This  is  followed  by  a  similar  stack  buildup  for  the  10-disc  generator 
section.  The  assembly  is  topped  with  the  proximitor  reference  plate, 
load  cell,  spacer,  and  torque  isolation  yoke,  which  are  all  located 
by  ilie  preload  adjusting  screw'.  The  arms  of  the  torque  isolation 
yoke  fit  closely  into  grooves  in  the  load  frame,  which  enables  the  pre¬ 
load  to  be  varied  without  applying  rotational  components  to  the  stack. 
Precision  bore  slotted  Pyrex  tubes  were  used  to  maintain  stack  concen¬ 
tricity  during  assembly,  and  then  broken  away  on  a  previously  scribed 
mark  after  establishing  the  preload.  This  tube  also  locates  up  to  the 
load  pad  on  each  end  of  the  stack,  which  is  in  turn  located  by  the  central 
pin.  A  typical  contract  stack  is  shown  in  Figure  2-34  alter  mechanical 
assembly,  but  prior  to  completing  electrical  and  instrumentation  con¬ 
nections:.  The  Pyrex  alignment  tube  sections  are  still  in  place,  and  the 
edges  can  be  seen  bn  both  sides-  of  the  copper  contact  terminals. 


The  fully  annealed  copper  contacts,  and  especially  the  compliance  pads, 
serve  to  compensate  for  surface  irregularities  and  any  tilling  that  does 
occur  in  the  final  assembly.  These  pads  yield  at  about  9000  psi  com¬ 
pressive  stress,  and  generally  conform  well  to  nonuniform  loading 
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Figure  2-32.  Typical  Contract  Life  Test  Fixtute  Components 
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Figure  2-33.  Load  Frame  Assembly  Details 
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where  the  localized  stresses  can  reach  higher  levels  than  the  initial 
average  preload  of  12,  000  psi  used  to  seal  in  the  stacks.  This  preload 
significantly  reduces  as  stack  elements  relax  and  is  repeated  several 
times  during  initial  assembly  until  yielding  essentially  stops,  as 
indicated  by  load  cell  stress  remaining  high  for  several  minutes.  O'ice' 
this  condition  is  achieved,  the  preload  is  never  fully-: relaxed  until  the 
hig-h-dielect  ric  strength  potting  compound  is  ftarderUTcl  in  place  to 
effectively  maintain  all  relative  locations,  with  or  without  preload. 

The  encapsulating  technique  is  described  more  fully  below  in 
Section  2.  4.  2.  2. 

•  \  •  •  •  '  v  '  : .  '  •  '  i  , 

1  \  '■  '  \  v  !:  ;  I 

2.  4.  1.  4  Driver  Stack  Power  Supply  and  Protection  Circuit 

The  best  means  for  achieving  a  large  voltage  swing  on  the  piezoelectric 
driver  circuit  is  to  allow  the  supply  voltage  to  swing  partly  negative. 

This  reduces  the 'highest  absolute  potential  experienced  by  the  stack. 
Based  on  vendor  recommendations  (Reference  5),  a  negative  swing  of 
10  volts/mil  or  400  volts  is  acceptable  before  depolarization  is  initiated. 
The  55  volts/mil  or  2200  volt  peak-to-peak  power  supply  (Figure  2-35), 
was  therefore  normally  operated  with  a  700  volt  dc  bias  and  an  1100 
volt  ac  amplitude  to  produce  voltage  swing's  from  -400  to  +1800  volts. 

Two  such  circuits  were  used:  one  was  designed  to  handle  three  stacks 
in  parallel,  while  the  other  drives  two  stacks  in  parallel.  Ac  line 
power  is  applied  through  a  variable  autotransformer  to  the  primary 
of  a  30  to  1  step-up  transformer.  Because  the  dc  power  supply  does 
riot  pass  reverse  current,  a  low  impedance  ac  shunt  is  provided  by  the 
high-voltage  1.6  (uf  capacitor. 

The  70-disc  piezoelectric  driver  stacks  are  connected  as  seven  10-disc 
segments,  each  of  which  is  represented  by  a  capacitance  Cg.  Any 
segment  can  be  easily  removed  from  the  circuit  in  case  of  individual 
failure.  Protection  against  high  current  surges  from  complete  voltage 
breakdown  across  a  disc  is  provided  by  resistors  R„.  The  value  of 
R  is  125K  for  the  10-disc  segments  of  the  contract  stacks  and  75K 
for  the  20-disc  segments  of  the  IRAD  driver  stack.  The  variable  re¬ 
sistor  was  included  in  only  one  power  supply  circuit  which  used  a  step- 
change  dc  power  supply  to  limit  voltage  buildup  rates  across  Cs  during 
dc  supply  changes.  The  resistance  was  decreased  to  zero  before  ap¬ 
plying  ac  excitation  so  that  essentially  the  entire  dc  voltage  drop  appears, 
across  Cg.  The  two  power  supply  and  protection  circuit  modules  (not 
including  the  dc  supply)  are  shown  respectively  in  Figures  2-36  and  2-37. 


2.4.  1.  5  Generator  Stack  Electrical  Load  and  Instrumentation  System 


A  piezoelectric  generator  may  be  schematically  represented  by  a  volt¬ 
age  source  in  series  with  an  internal  resistance  and  capacitance.  For 
both  the  parametric  tests  and 'life  tests,  generator  output  was  applied 
to  a  resistance  load  with  a  series  inductance  used  to  compensate  for 
generator  capacitance  by  resonantly  tuning  the  circuit  (Figure  2-38), 
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Generator  output  is  determined  as  V^/R^  where  V j__  is  the  voltage 
measured  across  the  load  resistance  “R \^.  Maxin]«=m  power  to  the  load 
is-aehieved  when  the  load  resistance  R  j-  Ls~equalto  the  generator  stack 
internal  resistance  Rs.  As  with  any  generator,  this  is  an  undesirable 
operating  condition  because  half  the  generated  power  is  dissipated  in 
heating  the  generator,  and  a  maximum  of  50%  efficiency  can  be  achieved 
(if  mechanical  stack  losses  are  negligible).  To  achieve  a  reasonable 
compromise  in 'the  tradeoff  between  power  output  and  efficiency,  load 
resistances  of  120K  (three  times  the  value  for  maximum  power  trans¬ 
fer)  were  used  for  the  contract  life  test  stacks.  The  smalle r-diameter 
IRAD  generator  stack  had  a  maximum  power  load  of  200K,  and  540-K 
resistors  were  used  during  the  life  tests.  ; 

Optimum  inductance  (to  give  maximum  power  transfer  with  the  operat¬ 
ing  load  resistance)  for  the  contract  stacks  ranged  from  1 1  5  to  13  5 
henries  during  the  parametric  tests,  but  was  measured  at  17  5  henries 
for  the  actual  operating  conditions  of  the  life  tests.  An  optimum  value 
of  840  henries  was  experimentally  determined  for  the  smaller  diamete  r 
IRAD  generator  stack. 
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Figure  2-35.  Driver  Stack  Power  Supply  and  Protection  Circuit 
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Figure  2-38.  Typical  Generator  Stack  Electrical  Load 


Instrumentation  to  monitor  stack  performance  included  Tektronix  Model 
P60I3  lOOOx  attenuator  probe  with  a  100  megohm  input  impedance  used 
in  conjunction  with  a  Tektronix  Model  545A  oscilloscope  to  measure 
voltage  applied  to  the  driver  stack  and  output  at  the  terminals  of  the 
generator  stack.  This  not  only  reduced  the;  measured  voltages  to 
manageable  levels,  but  was  the  only  device  available  with  sufficiently 
high  input  impedance  to  avoid  loading  the  circuits  being  monitored. 

The  load  voltage  across  only  was  measured  with  a  Ballantine  Model 
3  55  digital  voltmc-cer. 

Output  from  the  Kistler  Model  901A  quartz  load  ceil  mechanically  in 
series  with  the  stack  was  fed  through  a  Kistler  Model  566  or  504D 
(depending  on  the  stack)  charge  amplifier  into  another  Ballantine  Model 
3  55  digital  voltmeter.  Output  is  directly  proportional  to  axial  force 
and  therefore  also  to  stack  stress. 

The  strain  for  each  generator  stack  was  measured  with  two  Bently- 
Nevada  proximitor  probes  connected  to  Model  3601  proximitors. 
Physical  relationships  between  the  proximitor  probes  and  metal  ref¬ 
erence  surfaces  are  shown  in  Figure  2-33.  The  proximitors  we;e 
calibrated  prior  to  assembly  using  a  Pratt  and  Whitney  Model  G-2100 
Super  micrometer. 
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C.ycltc  work  input  was  determined  by  connecting  the  proximitor  and 
^charge  amplifier  outputs  respectively  to  the  x-  and  y-axes  of  a 
lektronix  Model  RM503  x-y  oscilloscope.  This  produced  a  force- 
d  i  splncement  diagram  for  the  generator  stack  which  was  recorded 
photographically  and  the  enclosed  area  was  later  determined  to  quantify 
, rrrcchan i ca  1  input  to  the  generator  stack.  The  measured  electrical 
output  divided  by  this  mechanical  input  produces  an  efficiency  measure¬ 
ment  for  the  'generator  stack. 

The  life  test  equipment  and  instrumentation  are  shown  in  Figures  2-39 
to  2-41.  Figure  2-29 -includes  the  right  side  of  the  test  facility  including 
the  low  temperature  stacks  inside  the  freezer,  dual  beam  oscilloscope, 
[RAD  test  fixture  on  the  bench  top,  Kepco  Model  1250B  dc  power  supply, 
and  one  power  supply  and  protection  circuit  system.  The  left  side  of 
the  laboratory  (Figure  2-40)  includes  the  Fluke  Model  412B  dc  power 
supply  on  t  op  of  the  I  30  ■  F  oven,  two  digital  voltmeters  on  the  x-y 
oscilloscope,  and;  both  power  supply  and  protection  circuit  systems. 

I  igure  2-41  shows  the  oven  interior  with  the  high-temperature  stack. 
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Figure  2-39.  Life  Test  Laboratory  Showing  Freezer  and  IRAD  Stacks 
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Piezoelectric  Disc  and  Contact  Processing 


Elaborate  cleaning  and  processing  was  used  for  the  piezoelectric  discs 
>  p.f  copper  contacts  lor  the  liie  test  assemblies.:  This  action  was  based 
on  xesp.ts  reported  by  Physics  Internationa  l-(  Refosifnce  7)  for  their 
piezoelectric  motor  driven  artificial  heart  system  and  on  earlier  fail¬ 
ures  experienced  by  DWDL.  The  primary  objective  of  the  procedures 
is  to  eliminate  contaminants  and  sharp  points  which  may  initiate  electric 
discharges  that  produce  carbon  tracks  and  eventually  lead  to  complete 
breakdown.  Another  secondary  effect  of  foreign  particles  between 
discs  and  contacts  or  ot  hard,  non-annealed  copper  contacts  is  to  in- 
ci  ease  system  compliance  which  reduces  the  potential  stress  swing. 


An  example  oi  catastrophic  piezoelect  ric  disc  failure  is  shown  in 
I-  igure  2-42.  Fused  ceramic  and  carbon  deposits  are  visible  where 
the  encapsulant  has  been  removed.  This  failure  occurred  in  a  stack 
with  the  same  size  and  type  of  discs  as  used  for  .the  life  test  assemblies. 
It  was 4  howeve  r,  received  from  the  supplier  as  a  bonded  assembly, 
then  cleaned  and  potted  with  the  silastic  used  as  the  secondary  encap¬ 
sulant  for  the  life  test  stacks.  The  stack  was  operated  under  no  load 
for  about  20  hours  at  peak-to-peak  voltage  swings  up  to  2000  v.  In¬ 
ductance  was  used  in  the  driver  circuit,  which  stored  sufficient  energy 
to  cause  even  the  ceramic  to  melt  in  the  discharge  area.  A  similar 
breakdown  occurred  in  the  first  assembly  of  the  I  RAD  prototype  stack. 


Figure  2-42.  Catastrophic  Failure  of  Bonded  and  Potted  Test  Stack 
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These  components  had  been  thoroughly  cleaned  and  deg-reased,  assembled 
in  a  laminar  flow  bench,  and  the  assembly  operated'  submerged  in  high- 
dielectric  transformer  oil.  .This  unit  operated  only  on  an  intermittent 
basis  v  ith  numerous  minor  audible  discharges  culminating  in  cata¬ 
strophic  failure  |When  voltage  swings  reached  1  500  volts  peak-to-peak. 
Failure  was.  attributed  to  foreign  particles  from  the  oil,  because  some 
such  particles  were  removed  from  adjacent  parts  ofsEhe  stack  prior  to 
fina4-failure.  Improvement  was  obviously  needed  to  achieve  a  2000- 
hour  life  test  with  peak-to-peak  voltages  of  2200  volts  or  55  volts  per 
mil. 

2.  4.  2.  1  Processing  and  Assembly  of  Life  Test  Stacks 

The  processing  procedure  actually  used  for  the  life  test  stacks  was  an 
extensively  modified  version  of  that  developed  by  Physics  International 
(Reference  7).  Many  steps  developed  at  DWDLfor  processing  of  ther¬ 
mionic  diode  components  were  incorporated.  This  procedure  was  first 
applied  to  the  third  assembly  of  the  IRAD  prototype  stack  for  proof 
testing,  although  schedule  committments  required  completion  of  the 
contract  disc  and  contact  processing  prior  to  actual  operation  of  the 
IRAD  stack.  This  procedure  represents  extreme  precaution  to  provide 
maximum  confidence  in  successful  life  testing.  It  is  likely  that  con¬ 
siderably  simpler  procedures  would  prove  satisfactory.  The  sequence 
for  the  piezoelectric  discs  is  summarized  as  follows. 

1.  Discs  are  stacked  and  clamped  with  flat  surfaces  together  for 
surface  protection,  sandblasted  on  the  edges  with  20-micron 
aluminum-oxide  grit,  and  flushed  with  tap  water.  This  removes 
metalizing  left  from  electric  contact  deposition  and  other  con¬ 
taminants  from  the  edges. 

2.  Warm  watet^  and  Alconox  wash  in  ultrasonic  cleaner  with  solu¬ 
tion  being  recirculated  and  filtered  (15  minutes). 

3.  Continuous  flow  flushing  rinse  with  tap  water  in  ultrasonic 
cleaner  (3  0  minutes). 

4.  Distilled  water  rinse  in  ultrasonic  Cleaner  with  water  being 
reei  rqulated  and  filtered  (15  minutes).  Figure  2-43  shows 
some  of  the  discs  at  this  stage  of  processing.  The  black  marks 
on  the  discs  indicate  pola  ;  tty. 

5.  Trichloroethylene  rinse  in  ultrasonic  cleaner  (10  minutes). 

6.  Reagent  alcohol  rinse  in  ultrasonic  cleaner  with  fluid  being 
reci  rculal ed  and  filtered  (15  minutes). 

7.  Freon  rinse  in  ultrasonic  c  leaner  with  fluid  being  recirculated 
and  filtered  (15  minutest. 
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8.  Vapor  degrease  with  trichloroethylene  (1  hour). 

9.  Freon  rinse  in  ultrasonic  cleaner  with  fluid  being  recirculated 
and  filtered  ( 1  5  minutes).  Repeat  with  fresh  fluid  ( I  5  minutes). 

10.  Vacuum  bake  at  150°F  with  discs  shorted  (overnight). 

11.  Recheck  polarity  as  required  by  loss  of  identifying  marks 
during  processing. 

12.  Store  in  cl'ean  Petri  dishes  until  assembly.  Examples  of  this 
for  both  discs  and  completed  copper  contacts  are  shown  in 
Figure  2-44. 

Preparation  of  0.002-in  thick  copper  contacts  is  directed  toward  mini- 
mi/.ing  both  electrical  breakdown  probability  and  stack  compliance. 

The  copper  contacts  are  put  through  the  following  operations. 

■  1..  Freon  rinse  in  ultrasonic  cleaner  with  fluid  being  recirculated 

and  filtered  (15  minutes). 

2,  Discs  are  rolled  flat  on  granite  surface  plate. 

3.  Deburr  in  kitchen  cleanser  and  water  slurry  using  ultrasonic 
cleaner  with  mixture  being  recirculated  (5  hours). 
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Figure  2-44.  Fully  Processed  Piezoelectric  Discs  and  Copper  Contact  Assemblies 


Continuous  flushing  rinse  with  tap  water  in  ultrasonic  cleaner 
(45  minutes). 

Flush  with  distilled  water,  reagent  grade  alcohol,  and  Freon 
(as  required). 

Rinse  with  Freon  in  ultrasonic  cleaner  with  fluid  being  recir¬ 
culated  and  filtered  (Z  hours). 

Rinse  with  reagent  alcohol  (as  required). 

Rinse  with  Freon  in  ultrasonic  cleaner  with  fluid  being  recir¬ 
culated  and  filtered  (15  minutes). 

Rinse with  reagent  alcohol  and  ethyl  alcohol  (as  required). 

Outgas  with  mechanical  vacuum  pump  and  liquid  nitrogen  cold 
tap  ( ove might). 

Anneal  in  10~^  torr  vacuum  at  700*0(5  minutes).  Contact;' 
are  shown  in  the  vacuum  furnace  just  prior  to  annealing  in 
Figure  Z  -45, 

Solder  contact  tabs  together  in  groups  of  five. 


i  •'  • 
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[  13.  Remove  residual  acid  flux  with  boilinjg  distilled  water,  assisted 

’  '  ;h  ‘  by  a  camel  hair  brush,  while  clamped  in  assembly  fixture. 

|  |  r-  14.  Flush  with  distilled  water,  alcohol,  aind  Freon  (as  required), 

f  I  ■  '  •  .  •  ..... 

|  15.  Remove  from  assembly  fixture  and  store  in  Petri  dish  until 

|  ;  final  assembly. 

I  ;  The  piezoelectric  stacks  are  assembled  electrically  in  parallel  by  placing 

:  i  discs  on  edge  with  alternating  polarity  and  interleaving  the  copper  con- 

i  tact  subassemblies  to  form  seven  10-disc  segments.  Initial  buildup 

is  accomplished  in  a  semicylind rical  section  of  precision-bore  Pyrex 
);  ■  tubing.  The  10-disc  generator  stack  or  70-disc  driver  stack  is  then 

transferred  to  a  slotted  Py rex-tube  handling  fixture  which  encloses 
i  about  two-thirds  of  (he  stack  circumference  by  placing  the  Pyrex  sec¬ 

tions  end  lo  end  in  a  steel  vee-block  and  sliding  the  stack  assembly 
from  one  to  the  other  (Figure  2-46).  The  copper  contact  tabs  for  all 
i  seven  sections  of  one  polarity  have  been  soldered  together;  the  sections 

on  the  other  side  are  independently  connected  to  the  power  supply  and 
i  '  ircuit  protection  equipment. 

A  closer  view  of  the  reverse  side  of  the  stack  through  the  glass  sleeve 
is  shown  in  Figure  2-47.  Surface  irregularities  in  the  disc  edges  and 
h  !•<.;).»  misalignments  are  visible.  In  two  or  three-  eases,  copper  eon  - 

Cl 
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Figure  2-46.  Stack  Assembly  in  Pyrix  Sleeve  and  Vee-B)ock 
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2.  4.  2.  2  Dielectric  Environment  for  Life  Test  Stacks 

It  is  imperative  that  a  high  dielectric  strength  environment  be  provided 
around  the  stack,  because  voltage  gradients  are  much  higher,  than  those 
needed  to  ionize  air  and  initiate  an  electrical  discharge.  Oil  is  typi- 
_cally_used  for  this  type  of  application,  and  washed  by  DWDL  for;  early 
work,  including  the  first  assembly  for  operational  testing  of  the  IRAD 
prototype  stack.  Problems  associated  with  using  oil  in  the  present 
application  include. the  following. 

1.  High  dielectric  oils  are  typically  somewhat  hygroscopic,  and 
probably  accumulate  an  unacceptable  amount  of  moisture  over 
a  2000-hour  open  atmosphere  test. 

2.  Tools  used  for  preload  adjustments  during  data  taking,  and 
other  aspects  of  practical  operation,  make  contamination  of 
oil  with  foreign  matter  difficult  to  avoid. 

3.  Most  importantly,  use  of  the  piezoelectric  load  cells  requires 
reducing  stack  preload  to  zero,  for  an  absolute  reference 
point  each  time  data  are  taken.  Stack  slippage  is  possible 
where  it  is  not  externally  held  in  place. 

These  considerations  Led  to  an  attempt  to  find  an  inert  high-dielectric- 
strength  potting  compound  with  sufficient  flexibility  to  accommodate 
stack  vibrations.and  rigid  enough  to  maintain  relative  positions  of  com¬ 
ponents  with  Or  without  the  preload  screw  in  place.  The  basic  encap- 
sulant  selected  was- Dow  Corning  Sylgard  184  with  a  dielectric  strength 
of  500  volts /mil.  Encapsulation  of  a  test  stack,  however,  resulted  in 
a  poor  bond  between  the  stack  and  Sylgard.  Accordingly,  a  5-  to  10- 
mil  thick  layer  of  Dow  Corning  QR-4-3117  circuit  board  coating  with 
a  dielectric  strength  of  2000  volts/mil  and  excellent  adherence  to  copper, 
and  ceramic  was  used  as  the  primary  encapsulant  with  Sylgard  184 
j  applied  for  final  encapsulation.  This  combination  of  encapsulants  tested 
!  with  smooth;  metal  spheres  exhibited  no  breakdown  when  20,  000  volts 
were  applied  with  a  10-mil  separation.  Some  microscopic  bubbles 
were  observed  in  the  primary  encapsulant  applied  to  the  stacks.  Im¬ 
proved  potting  techniques  reduced  the  problem  with  later  stacks,  but 
it  is  possible,  that  some  of  the  observed  breakdowns  in  driver  segments 
might  have  resulted  from  a  series  of  such  bubbles.  Results  of  this 
encapsulating  approach  have  been  basically  good  and  definitely  superior 
to  oil. 


2.  4.  3  Parametric  Tests 

The  primary  purpose  of  the  -parametric  testing  was  to  establish  base¬ 
line  generator  stack  performance  characteristics  and  optimum  load 
values  for  one  standard  and  one  optically  lapped  stack  using  equipment 
with  known  and  controllable  stress  application.  These  objectives  were 


i 

1 
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!  t  n  -  ^  stlB,u,-y  l,uyuml  Uly  flises  and  ar-u  bunt 

Mov, nwrtrrl.  ;  (Mh  rondiiion  (icrumitl  only  in  thu  coni  rn  cl  starks;  For 

nHH.'n.bly  o|  the  prototype  stack,  the  entrance,  end  of  the  elans  tube 
VW,M  )';r  ensler  InMcrtion,  Thin  wan  one  of  several  minor 

upei'nttoiiH  will  eli  were  rmnltflufi  or  e  !1  mi  tittle  H  lo  expedlle  awwomhly  of 
_uie  eohlract  load  f fame  a,  and  which  uene rally  thieved  the  deni  red 
enect.  In  this  instance  however,  the  result  waTan"  unacceptable  con¬ 
dition  where  some  of  the  copper  contacts  had  portions  severed  off  and 
smeared  onto  the  edge  of  the  piezoelectric  discs.  The  extent  of  damage 
is  shown  in  Figure  2-48.  in  the  region  at  one  edge  of  the  Pyrex  sleeve 
where  most  problems  occurred. 

After  considering  several  rework  alternatives,  the  damaged  stacks 
were  assembled  into  the  load  frames,  preloaded  to  12,  000  psi,  align¬ 
ment  s i e c \  e s  r e mo ved ,  and  sandblasted  locally  to  remove  .debris 
Subsequent  microscopic  examination  indicated  removal  of  copper*(and 
substantial  ceramic)  from  the  disc  edges,  but  some  sand  particles 
remained  wedged  between  disc  interfaces;  Such  sharp  edges  are 
sources  of  localized  electric  field  concentration,  but  the  ultimate  di¬ 
electric  environment  was  expected  to  be  capable  of  sustaining  such 
fields  without  breakdown. 


accomplished  on  schedule.  Original  .mentions  were  to  repeat  the  para¬ 
metric  tests  on  surviving  generator  stacks  following  life  testing.  After 
technical  problems  delayed  the  initiation  of  life  tests,  however,  comple¬ 
tion  of  the  final  parametric  tests  would  have  required  pr'ematufe' 
structive  termination  of  the  life  tests  to  meet  the  schedule.  This  ftask;, 
was  accordingly  deleted  by  mutual  as  reement  to^accumulate  morel  time 
Tn  life  testing. 


..4.3.  1  Electrohycraulic  '!  e s ■  :■<  •  r.t 

The  e lectrohyd raulic  materials  loading  system  used  for  the  parametric 
tests  was  leased  from  Battelle  Pacific  Northwest  Laboratories  in 
Richland,  Washington.  Details  of  the  test  equipment  and  associated 
instrumentation  are  described  in  AppendixiC.  The  test  configuration 
used  for  the  parametric  tests  is  shown  in  Figures  2-49,  2 -50,  and  2-51 . 
Figure  2-49  shows  the  load  system  with  electrical  leads  from  the  test 
stack  and  proximitors  leading  to  the  electricall  load  and  instrumentation 
racks.  The  test  stack  in  place  appears  in  Figure  2-51;  electric  cor.nec- 
-ors  from  the  1 0-disc  generator  stack  obscure  the  proximeters  used  for 
displacement  measurements  aero-s  the  stack.  These  are  more  clearly 
shown  in  Figure  2-51  showing  the  opposite  side  of  the  test  stack. 


Figure  2-49  Carsmptnc  Test  Setup 


2.  4.  3.  2  Parametric  Test  Results 

All  parametric  tests  were  run  with  peak  stress  levels  of  9000  psi,  and 
were  preloaded  to  achieve  minimum  stress  levels^of  5000,  3000,  and 
KLOO  psi.  These  varied  conditions  resulted  in  riffpective  stress  swines 
of  4000,  6000,  ant!  S000  psi.  Load  resistances  were  varied  over  a 
wide  range  and  inductance  was  varied  around  the  expected  optimum. 
Highlights  of  the  test  results  are  summarized  in  Table  2-5.  The  data  in 
/Table  2-  5  is  fo  r  ope  ration  at  the  optimum  load  inductance  for  each  stack 
ana  loao  resis  Lances  of  60K  and  80K  as  noted.  This  improves  efficiency 
substantially  above,  the  maximum  power  transfer  condition  at  a  load  re¬ 
sistance  of  40 K.  The  effects  of  varying  load  resistance  and  induc  tance 
on  power  density  and  efficiency  for  the  different  stress  swings  are  shown 
in  Figures  2-  53  through  2-  57.  Efficiency  of  the  ;optically  lapped  stack  was 
very  close  to  50%  in  each  case  where  the  power  output  reached  a  maxi¬ 
mum  value,  or  the  load  resistance  was  equal  to  the  internal  resistance 
of  the  stack.  The  breadth  of  the  maximum  power  peak  and  amount  of 
data  prevent  an  accurate  assessment  of  exactly  where  the  peak  occurs, 
but  the  efficiency  is  clearly  not  far  below  50%,  which  shows  there  are’ 
only  small  mechanical  losses  within  discs  and  at  interfaces.  For  the 
standard  lapped  stack  in  Figure  2-54,  peak  power  occurs  for  a  load 
resistance  of  20K  to  30K.  The  efficiency  in  this  range  varies  from  32% 
to  40%.  Since  the  electrical  efficiency  at  this  operating  pointTs  SO”?,, 
the  reduced  efficiency  is  caused  by  mechanical  losses  ••within  'the  stack. 

Table  2-5 

SUMMARY  OF  PARAMETRIC  TEST  RESULTS  AT  60  HZ  WITH 
OPTIMUM  LOAD  INDUCTANCE 
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Figure  2-52.  Performance  Dependence  on  Load  Resistance  for  Optically  Lapped  Stack 


Figure  2-53.  Performance  Dependence  on  Load  Inductance  for  Optically  Lapped  Stack 
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The  efficiency'  range  of  32%  t<*  40%  means  that  me  chan  real  losses  are 
between  50%  and  1.12%  of  the  electrical  power  actually  delivered  to  the 
load.  It  is  these  substantial  mechanical  interface  tosses  which  cause 
the  maximum  measured  efficiency  of  63%  in  Figure  2-54  to  be  near  the 
asymptotic  limit. 

2.  4.  4  Life  Test  Results 

Life  test  conditions  and  test  results  are  summarized  in  this  section. 

2.4.4.  I  Operating  Conditions 

Electric  excitation  of  the  life  test  fixtures  was  initiated  as  each  assembly 
was  completed.  Varying  amounts  of  time  were  spent  at  relatively  low 
input  voltages  lo  provide  stack  conditioning.  Emphasis  was  on  initiating 
life  testing  ot  alt  units  as  early  as  possible,  and  complete  instrumenta¬ 
tion  was  not  added  to  all  units  until  some  had  accumulated  several 
hundred  .hours,  A  brief  summary  of  accumulated  times  at  various 
conditions  is  reported  for  each  stack  in  Table  2-6. 

Normalized  power  density  and  efficiency  arc  plotted  for  the  operational 
times  of  the  life  test,  as  appropriate  to  each  stack,  in  Figures  2-5S 
through  2-61.  Life  test  hours  do  not  necessarily  agree  withftotal  accu¬ 
mulated  hours  in  Table  2-6,  since  the  life  test  hours  were  only  counted 
from  the  time  when  full  stress  was  applied  to  the  stack  at  the  specified 
environmental  condition.  No  life  test  plot  is  shown  for  stack  No.  3 
because  it  was  physically  damaged  prior  to  installation  of  full  instru¬ 
mentation  for  making  efficiency  measurements.  The  normalization 
factor  used  in  this  data  presentation  was  to  divide  output  power  density 
by  the  square  of  the  developed  stress  swing.  This  accounts  for  varia¬ 
tions  in  applied  stress  swing  as  a  result  of  different  driver  voltages  or 
reducing  the  number  of  driven  stack  segments  when  short  circuits  devel¬ 
oped.  The  quad  ralie  dependence  corresponds  to  the  simple  theory  in 
Appendix  B,  and  correlates  well  with  measurements  over  a  wide  range 
of  conditions.  '> 

both  power  density  and  efficiency  of  the  1RAD  stack  exhibit  erratic  be¬ 
havior  over  the  period  from  about  ‘*00  to  1200  hours  into  the  life  test 
(Figure  2- 38),;  including  some  effic iency  measu remon ts  above  100%.  The 
actual  variations  in  powe  r  density  ove  r  this  period  are  only  on  the  order  of 
7%,  although  these  a  re  somewhat  magnified  by  the  use  of  an  expanded  scale. 
The  efficiency  variations  are  much  more  severe,  however,  and  suggest 
a  temporary  error  in  instrumentation  calibration  or  utilization.  The 
circles  and  triangles  represent  efficiency  as  measured  by  the  two  prox- 
imilors  on  both  sides  of  the  stack.  Except  for  the -erratic  data,  the 
correlation  between  the  two  proximitors  was  generally  good,  indicating 
very  little  tipping  of  the  proximilor  mounting  plate  or  reference  plate. 

The  long-terni  trend  in  normalized  power  density  for  tin*  TRAD  stack 
show's  less  than  a  3".  reduction.  The  efficiency  trend  was  essentially  flat. 
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■!-■  .  Table  2-6 

SUMMARY  OF  LIFE  TEST,  RESULTS  THROUGH  7  SEPTEMBER  1973 
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Figure  2-58.  I  RAD  Stack  Efficiency  and  Normalized  Power  Dependence  on  Elapsed  Time 


Figure  2-59.  Stack  No.  1  Efficiency  and  Normalized  Power  Dependence  on  Elapsed  Time 


Figure  2-60.  Stack  No.  2  Efficiency  and  Normalized  Power  Dependence  on  Elapsed  Time 


Freezer  Stacks  No.  1  and  2  (Figures  2-59  and  2 -6 1 ,  respectively) 
again  showed  quite  stable  power  output,  with  scatter  in  the  efficiency 
data.  The  increased  power  density  at  about  1000  hours  for  Stack  No.  1 
and  800  hours  for  Stack  No.  2  are  real  effects  associated  with  de¬ 
creasing  the  ambient  temperature  below  -25°F.  The  normalized  power 
density  for  the  lapped  stack  is  about  23%  greater  than^that  for  the 
standaxd..  stack.  Power  density  in  Figure  2-61  again  Sows  about  a 
3  5%  decrease  for  the  standard  oven  stack  as  compared  with  the  standard 
freezer  stack,'  verifying  the  power  density  trend  indicated  by  subcooiing 
the  freezer  stacks.  As  an  additional  check  on  temperature  sensitivity, 
Stack  No.  2  was  transferred  to  the  oven  after  2085  hours  in  the  freezer. 
Power  density  decreased  rapidly  for  the  first  few  hours  of  operation  at 
150°F,  then  stabilized.  The  overall  reduction  in  power  density  was 
about  35%  after  3  40  hours  in  the  oven,  although  a  reduced  stress  swing 
resulted  in  only  a  10%  reduction  in  normalized  power  dens ity.  No 
electrical  breakdown  was  observed,  substantiating , the  theses  that  the 
breakdowns  which  occurred  during  initial  operation  of  Stack  No.  4  at 
1  50 °F  were  a  result  of  specific  stack  processing  problems  rather  than 
;  a  function  of  the  elevated  temperature.  A  significant  trend  toward 

j  power  density  degradation  is  shown  for  the  oven  stack  (Figure  2-61) 

•;  where  the  normalized  power  density  decreased  about  15%  in  800  hours, 

t  then  remained  essentially  constant. 

i  •“  The  scatter  of  efficiency  data  (Figures  2-59  and  2-60)  for  the  two 

\  proximitors  indicates  that,,  in  contrast  to  the  IRAD  stack,  considerable 

\  tippirig  of  the  mounting  plate  for  the  proximitors  or  reference  screws 

was  occurring.  Similar  relationships  between  the  two  proximitors  on 
l  each  stack  were  usually  maintained,  but  in  some  instances  the  tipping 

effect  appeared  to  reverse.  This  could  have  been  caused  by  minor 
position  changes  occurring  when  the  preload  was  reduced  for  data 
collection.!  The  general  efficiency  trends  (Figures  2-59  to  2 -6 1 )  range 
f  from  essentially  constant  to  a  modest  increase  in  efficiency  with  accu¬ 

mulated  operating  lifetime. 
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Section  3 
CONCLUSIONS 


As  a  result  of  the  testing,  analysis,  and  mechanical  design  effort 
during  the  contract  period,  the  following  conclusions  can  be  derived 
with  re'spect  to  the  STEPZ  concept. 


1.  Piezoelectric  generator  life  tests  demonstrated  a  power 
density  of  10  watts/in.  with  80'C  efficiency  at  about  7500- 
psi  stress  fluctuation,  at  which  level,  less  than  5"’,.  degrada  - 
tion  occurred  during  3000  hours  of  operation  at  60  Hz. 
Successful  operation  was  demonstrated  at  temperatures 
from  - 5 5 ° F  to  +150°F.  Significantly  better  performance 
was  achieved  with  optically  lapped  piezoelectric  discs  than 
with  standard  discs.  Extensive  processing  and  clean  room 
handling  of  discs  was  used  in  the  life  test  assemblies,  but 
the  degree  of  cleanliness  required  for  successful  operation 
was  not  firmly  established. 


2,  Development  of  a  sophisticated  dynamic  simulation  detailing 
the  STEPZ  system  with  more  realism  than  previous  computer 
programs,  used  successfully  by  DWDL  for  modeling  actual 
engines,  leads  to  a  high  degree  of  confidence  in  the  analytical 
projections.  ;  Projected  performance  is,  in  general,  encour¬ 
aging  although  the ;  highly  resonant  operation  makes  the  sys¬ 


tem 


.  *  ✓ 
sensitive  to  load  changes.  A  nearly  constant  net  load 


will  probably  be  required 


Extrapolations  of  successful  hardware  approaches  to  higher 
operating  pressures  and  bearing  loads  involve  reliable 
engineering  design  techniques.  This  provides  confidence  in 
the  soundness  of  the  mechanical  design  of  this  engine.  The 
mercury  coupler  is  straightforward  in  concept,  but  does 
utilize  unique  ‘'natures.  Recent  related  experience  with 
corrugated  metal  diaphragms  in  a  conceptually  similar  system 
has  revealed  operational  problems  which  raise  questions  with 
regard  to  system  safety  for  the  mercury  coupler. 


Ill 
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Appendix  A 

PH  YSICS  INTERNATIONAL  RE!§0RT 
PIEZOELECTRIC  GENERATOR  TESTS 


ABOUT  TIMS  APPENDIX 

Appendix  A  consists  of  two  reports  submitted  by  Physics 
International  Company  to  document  their  efforts  on  the  sub¬ 
coni  ,r  i  to  test  piezoelectric  generator  performance.  The 
,,.  j  .1  report  discusses  results  and  conclusions  of  the  tests 
wit1  ,tn  •  Meet  rical  circuit  analysis  of  piezoelectric  generator 
operation.  The  supplementary  report  provides  additional 
backg  round  ■  i nlormat  ion,  analysis,  and  a  tabulation  of  test  data. 
Certain  areas  of  analysis  and  conclusions  are  based  on  pre¬ 
liminary  and  incomplete  information,  but  arc  included  here 
without  comment.  More  extensive  testing  and  analytical 
development  has  since  been  performed  by  DWDL  and  is  reported 
in  Appendices  B  and  C. 
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SUMMARY  AND  CONCLUSIONS 


Two  piezoelectric  generators  (stacks  of  piezoceramic  discs) 
were  tested  by  driving  them  hydraulically  up  to  12,400  psi  com¬ 
pressive  stress  at  fluctuations  of  60  Hz.  The  first  generator 
consisted  of  5  cubic  inches  of  PZT  . 5-H  material,  a  ferroelec- 
trically  "soft"  ceramic;  the  second  consisted  of  0.8  cubic 
inches  of  LTZ-1,  a  ferroelectricaily  "hard"  ceramic. 

The  PZT  5-H  generator  produced  electrical  power  as  pre¬ 
dicted  up  to  stresses  of  5000  psi,  but  the  piezoelectric  proper¬ 
ties  were  degraded  when  the  stresses  exceeded  5000  psi.  The 
output  at  5000  psi  was  6  watts  per  cubic  inch  of  PZ  material. 

The  LTZ-1  generator  produced  approximately  9  watts  per 
cubic  inch  at  9700,  psi. 

Data  above  9700  psi  was  invalid  because  of  excessive  fric¬ 
tion  in  the  hydraulic  '•.on! . 

It  is  concluded  that  16  to  20  watts  of  electrical  power  per 
cubic  inch  of  PZ  material  can  be  extracted  with  greater  than  90- 
percent  efficiency.  This  estimate  is  made  by  a  small  extrapola¬ 
tion  from  9700  psi  to  12,000  to  14,000  psi. 


SECTION  1 

I  INTRODUCTION 

This  report  covers  the  experimental  work  performed  to 
determine  power-density  limitations  of  a  piezoelectric  power 
generator  driven  hydraulically  at  low , frequency  (60  Hz). 

The  tasks  to  be  performed  were : 


1.  Design  and  fabricate  a  hydraulic  system  capable  of 
driving  a  stack  of  100  PZ  discs,  each  1.25  inch  in 
diameter  by  0.040  inch  thick,  with  pressure  pulses 
of  60  pulses/min.  at  zero  to  14,000  psi  compres¬ 
sive  stress. 

2.  Build  a  stack  of  100  PZT  5-H  discs  and  install 
into  an  existing  housing. 

3.  Test  PZT  5-H  generator  up  to  14,000  psi  or  to  its 
limiting  stress. 

4.  If  PZT  5-H  is  not  satisfactory  up  to  the  required 
stress,  test  a  second  material  with  a  higher 
stress,  rating. 


1.1  CHOICE  OF  PZ  MATERIALS 

-  I  : 

Table  1  shows  piezoelectric  properties  for  several  PZ 
materials  from  various  manufacturers.  The  power  output  for  a  PZ 
generator  with  a  pure  reslstanr-:  load  matched  for  maximum  power 
transfer  is  given  by 

W  =  T6  d3 3  933  °2  V 
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where: 


W  =  electrical  power,  watts  ' 


=  frequency,  radians  per  second 


=  PZ  charge  coefficient,  10 


-12  coul/M 


N/M 


2 


g,-,  =  PZ  voltage  coefficient,  — 

■  :  **k  n/m 

2  I  2 

a  =  compressive  stress  N/M. 

A 


V  =  volume,  of  PZ  material,  M' 


This  expression  is  derived  in  the  appendix. 


In  Table  1,  the  product  of  d^  g33  is  used  as  a  figure  of 
merit  for  PZ  materials  with  respect  to  power  output.  Note  that 
PZT  5-H  and  LTZ  2-H  have  the  highest  output  coefficient.  How¬ 
ever,  stability  of  the  piezoelectric  properties  of  a  material  in 
the  operating  environment  is  air,:?  important.  Trie  Curie  tempera¬ 
ture  and  the  coercive  field  have  a  bearing  on  this  stability. 

PZ  properties  become  unstable  as  the  temperature  approaches  this 

Curie  ooint.  Th&»- coercive  field,  K  ,  is  the  electric  field 

c 

that  will  cause  the  material  to  reverse  its  polarization.  Thus 

it  is  more  difficult  to  change  the  polarization  of  a  material 

with  a  high  K  than  one  with  a  low  E  .  Vernitron,  Inc.,  recom- 
c  c 

mended  a  maximum  compressive  stress  (cycled)  parallel  to  the 
polar  axis  for  each  of  its  materials.  Note  that  the  rated 
stress  for  P^T  5-H  is  2500  psi ,  compared  to  12,000  psi  for 
PZT- 4 .  Note  also  that  the  internal  losses,  are  lower  for  PZT-4. 
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From  Table  1  it  can  be  concluded  that  of  the  commercially 
available  PZ  materials,  the  ferroelectrically  "soft"  materials 
such  as  PZT  5-H  and  LTZ  2-1!  give  the  highest  output  at  a  given 
value  of  stress  fluctuation;  and  the  ferroelectrically  "hard" 
materials  will  withstand  higher  compressive  stresses  without 
their  piezoelectric  properties  being  degraded.  The  hard 
materials  are  also  less  lossy  than  the  soft  materials.  Physics 
International  has  fabricated  experimental  piezoelectric 
materials  that  appear  to  have  high  output,  high  Curie  tempera¬ 
ture,  and  high-coercive  field.  These  materials  have  not  yet 
been  fully  characterized,  but  would  be  worth  testing  in  a 
generator. 

Power  transfer  is  greater  when  the  capacitive  impedance  is 
balanced  by  an  inductive  impedance  rather  than  a  pure  resistance 
load  as  shown  below  (referring  to  Figure  1): 


or ,  when 


W  =  -j  rt 


E?  R. 


(zc  -  Zjf+fcs  *  V' 


the  power,  W,  is  maximum  when  Z  is  minimum  or  when  Zc  =  Z^. 


1 . 2  EXPERIMENTAL  APPARATUS 


A  schematic  of  the  hydraulic  system  and  PZ  generator  is 
shown  in  Figure  2.  Basically,  the  PZ  generator  is  alternately 
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■V.U1IUVL. 


J  ~  - - -  ana  to  a  vent:  with 

ovciy  half  rotation  of  the  rotary  valve.  -  T!^System  wan 

designed  to  use  on-hand  hydraulic  components  wherever  possible. 
The  rotary  valve,  however,  was  fabricated  for  the  test.  A 
Vickers  axial  piston  pump  with  a  displacement  of  0.065  cubic 
inches  per  revolution  was  driven  at  3450  rpm  to  supply 
hydraulic  pressure.  The  pressure  in  the  generator  was  moni¬ 
tored  with  a  Kictlor  607L  pressure  transducer.  Hydraulic  fluid 
was  stored  in  an  accumulator  close  to  the  rotary  valve  to  pro¬ 
vide  fast  pressurization  of  the  generator.  The  schematic  of  the 
generator  equivalent  circuit  and  the  output  circuit  is  shown  in 
Figure  1.  The  inductor  is  used  to  match  the  impedance  of  the 
generator's  capacitive  impedance.  Electrical  power  output  was 
determined  by  measuring  voltage  across  the  load  resistors. 
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SECTION  2 
RESULTS 


2.1  PZT  5f II  MATERIAL 

The  data  plotted  in  Figure  3  show  the  hydraulic  power  as  a 
function  of  compressive-stress  fluctuation  on  the  PZT  5-H  stack. 
The  numbers  indicate  consecutive  runs.  Power  output  generally 
decreased  with  each  run. 

After  run  5,  t:ve  stack  was  ropolarizcd  with  2-kV  pulsed  dc 
at  2  pulses  per  second.  Run  6  showed  a  higher  output  than  run 
5,  but  the  output  was  still  lower  than  run  1.  The  power  reached 
30  watts  at  6600  psi  and  then  dropped  to  26  watts  as  the  stress 
fluctuation  was  raised  to  9300  psi  .  When  the  stress  fluctuation 
was  reduced  to  4400  psi,  the  power  was  13  watts,  compared  to 
20  watts  for  the  same  stress  fluctuation  at  the  start  of  the 
run. 

The  projected  curve  was  obtained  hy  extrapolation  from 
earlier  tests  at  stresses  up  to  3500  psi  on  a  similar  stack. 

The  power  output  exceeds  the  projected  curve  for  stress  fluctua¬ 
tions  up  to  a  stress  level  of  5000  psi.  For  stresses  above 
5000  psi,  power  dropped  below  the  projected  curve.  Data  taken 
after  the  stack  had  been  driven  above  5000  psi  showed  lower 
power  output.  Power  output  could  be  raised  by  repolarizing  the 
stack  with  dc  voltage. 
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2:2  LT  7,-1  MATERIAL 


A  second  generator  was  built  with  10  discs  of  LTZ-1  (1. 25- 
inches  in  diameter  by  0.040  inch  thick),  a  material  with  more 
resistance  to  depolarization  but  with  theoretically  lower  output 
at  a  given  stress.  A  0 . 125-inch-thick  disc  of  the  same  material 
was  placed  at  each  end  of  the  stack  to  reduce  shear  stresses 
produced  by  the  Poisson's  ratio  effect  of  the  steel  parts  at  the 
ends  of  the  .stack. 


Power  odtput  of  the  LTZ-1  generator  is  shown  as  a  function 
of  applied  stress  fluctuation  (Figure  4).  Data  points  are  shown 
for  seven  runs.  The  last  data  point  taken  was  7  watts  at  9700 
psi.  Two  curves  are  drawn  through  the  data  points  to  cover  the 
spicad  of  the  data.  The  curves  are  extrapolated  to  12,000  psi. 
The  data  is  discussed  in  Section  3. 


SECTION  3 
DISCUSSION 


3.1  PZT  5-H  MATERIAL 


Figure  3  indicates  that  PZT  5-H  and  probably  other  ferro- 
electrically  "soft"  materials  cannot  be  driven  to  the  stress 
considered  for  the  STEPZ  generator.  This  material  did  perform 
as  predicted' for  short  periods  up  to  5000  psi,  which  is  twice 
its  rated  stress.  Even  at  low-stress  levels  it  was  necessary  to 
cool  the  stack  with  a  cooling  coil  in  the  insulating  oil.  This 
material  may  be  suitable  for  applications  that  require  a  short- 
duty  cycle  or  one-shot  power  from  a  small  package. 


Notice  that  there  is  some  scatter  in  the  data  that  increased 
with  increasing  stress.  There  was  a  0.013-inch-diameter  clear¬ 
ance  between  the  piston  and  housing.  The  O-ring  seal  tended  to 
extrude  into  this  gap  oven  though  there  was  a  nylon  backup  ring. 
This  could  cause  enough  friction  between  the  seal  and  the  hous¬ 
ing  to  explain  the  scatter,,  since  all  of  the  hydraulic  force  on 
the  piston  was  not  transmitted  to  the  stack. 


3.2.  LTZ-1  MATERIAL 


The  LTZ-1  material  (similar  to  PZT-4)  showed  no  signs  of 
losing  its  piezoelectric  properties.  Several  runs  were  made, 
showing  that  the  power  output  repeated  after  the  material  had 
been  driven  to  high  stresses.  The  data  points  in  run  6  are 
numbered  to  show  the  order  in  which  they  were  taken. 
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The  following  evidence  supports  •'the  thtfory  that  the  devi  a¬ 
tion  of  the  data  from  the  exponential  power  curve  at  stresses 
above  10,000  psi  was  due  to  a  defect  in  the  test  apparatus 
r  a  titer  than  to  degradation  of  P".  properties. 


1.  Marks  on  the  piston  and  the  housing  indicated  that  the 
O-ring  had  extruded  past  the  nylon  backup  ring  into  the 

0 .013-inch-diamoter  gap  between  the  housing  and  piston. 

2.  During  run  5  the  output  power  shifted  from  the  upper 
exponential  curve,  6  watts  at  7900  psi,  down  to  3.1  watts 
at  8SC0  psi.  The  subsequent  four  data  points  formed  a 

1  distinct  icurve  significantly  lover  than  the  exponential 
curves.  However,  the  data  in  run  6,  taken  the  following 
day,  generally  fall  between  the  two  exponential  curves 
up  to  9700  psi.  The  data  above  9700  psi  arc  lower  and 
scattered.  The  last  data  point  taken  lies  on  the  lower 
exponential  curve,  7  watts  at  9700'  psi. 

It  is  concluded  that  the  lower  power  curve  determined  in 
run  6  is  the  result  of  the  O-ring  exeruding  into  the  gap 
between  the  piston  and  housing,  causing  the  actual  load  on  the 
stack  to  be  significantly  lower  than  that  indicated  by  the  hy¬ 
draulic  pressure  on  the  piston.  The  scatter  in  run  5  is  due  to 
varied  degrees  of  high  friction  caused  by  the  extruding  O-ring. 


Assuming  that,  the  scatter  in  the  data  is  due  to  seal  fric¬ 
tion,  the  output  power  at  12,000  psi,  the  maximum  recommended 
stress  for  this  material,  would  be  13  watts  for  0.8  cubic  inches 
of  material.  Thirty-eight  percent  of  this  material  was  in  the 
two  0 . 125- inch-thick  end  discs.  It  is  probable  that  the  output 
would  be  greater  if  all  the  material  were  in  0 .040-inch-thick 
discs.  It  is  also  probable  that  the  material  could  be  driven 
at  higher  stress  than  12,000  psi.  Therefore,  it  is  estimated 
that  80  to  100  watts  of  power  can  be  extracted  from  5  cubic 
inches  of  1/13-1  material. 
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3.3  .  GENERATOR  EFFICIENCY 

The  efficiency  of  the  PZT  5-H  generator  had  been  determined 
in  an  earlier  experiment  to  be  70  percent.  The  power  loss  is 
considered  to  be  divided  into  three  mechanisms: 

1.  Electrical  and  mechanical  losses  in  the  material. 

2.  Losses  due  to  friction  between  the  PZ  stack  and  the 
stee]  end  pieces  from  radial  expansion  (Poisson's  ratio 
effect)  . 

3.  Seal  friction. 


The  PZ  material  loss  is  about  5  percent  for  PZT  5-H  and 
about  1  percent  for  LTZ-1.  It  is  estimated  that  seal  friction 
is  the  greatest  of  the  other  losses.  Therefore,  it  is  esti¬ 
mated  that  generator  efficiency  greater  than  90  percent  can  be 
achieved  by  replacing  the  O-ring  .with  a  metal  diaphragm. 


3.4  TEST  APPARATUS 


The  rotary  valve  was  designed  and  fabricated  to  drive  'the 
PZ  generator  for  testing  purposes.  It  became  apparent  during 
testing  that  the  system  had  several  limitations.  The  pressure 
applied  to  the  stack  was  pulsatile  rather  than  sinusoidal.  At 
high  pressures  the  system  required  a  higher  flow  rate  than  the 
pump  would  produce.  Therefore,  a  large  air/hydraulic  accumu¬ 
lator  was  necessary  so  that  hydraulic  energy  could  be  stored 
for  short  runs.  Also,  the  rotary  valve  became  hard  to  start 
at  hydraulic  pressures  above  2800  psi. 


r;.  - 
• 
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3.5  recommen  dat  lows 

It  is  recommended  that  small  stacks  of  Physics  International's 
own  materials  be  tested  to  assess  their  applicability  in  a  PZ 
generator.  It  is  also  recommended  that  a  longer  stack  of  LTZ-1 
material  ibe  tested  -in  such  a  way  that  the  .limitations  of  the  test 
apparatus  do  not  effect  the  results. 

A  piezoelectric  driver  is  quite  satisfactory  at  low  stresses, 
but  at  1 4,000  psi  the  driver  stack  would  have  to  he  3.5  times  as 
long  and  twice  the  diameter  of  the  driven  stack. 

It  is  also  recommended  that  the  hydraulic  driving  system  be 
modified  to  overcome,  the  difficulties  encountered  in  this  program. 
Those  modifications  would  include: 


1.  Using  with  a  low-friction  seal  or  petal  diaphragm,  a 
larger  piston  so  that  the  required  hydraulic  pressure  and 
seal  friction  are  lower. 


2.  Modifying  the  ports  in  the  rotary  valve  so  that  the 
pressure  pulses  would  more  closely  represent  a  sine  wave. 
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Electrical  oov/or  output,  watts 


0  1  2  3  4  5  6  7  8  9  10  11  12 

Stress  fluctuation,  1000  osi 


iquro  4  Output  power  of  LTZ1  generator,  volume  of  material 
0.8  in'. 
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PZ  GENERATOR  EQUIVALENT  CIRCUIT 
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An  equivalent  circuit  for  the  generator  is  shown  in 
Figure  A-l.  If  the  value  of  the  load  resistor  on  the  output 
is  small  compared  to  Rp ,  the  circuit  can  be  simplified  to  that 
shown  in  Figure  A-2.  Assuming  the  generator  output  is  sinusoidal, 
the  circuit  equation  in  the  frequency  domain  is: 


The  output  voltage,  V0,  can  be  measured  for  two  different 
values  of  load  resistor,  R.  .  The  load  resistor  must  be  small 

ij 

compared  to  Rp .  The  equation  can  Lo  written  as  follows: 


The  internal  series  resistors,  ,  is  now  given  by: 
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o 


rr  tho  generator  is  driven  at  the  frequency  w  ,  the 
of  .^c  can  he  inserted  into  the  equation,  'j^rerefore 

R.  Rr  \V'  -  V  \ 

R  .  = - i - £ _ 2: _ IS  _  _ 1_ _ TAw  u),  t 


VoJbuJu. 


S  "  Vr>  Rt  oj  C 

L2  °2  L1  1  s 


If  the  generator  ,\s  excited  at  frequency,  c>2 ,  can  again 
be  expressed  in  «*w*^equations?''" in  terms  of,  Z^,.  Since  the  value 
of  Rs  is  independent  of  frequency,  the  two  equations  can  be 
solved  for  R„  and  C„ . 


For  maximum  power  transfer  from  the  generator,  the  load 
resistance^ must  be  equal  to  the  internal  resistance,  R^ ,  and 
the  capacitive  reactance  must  be  balanced  by  the  addition  of  an 
inductive  reactance  of  the  same-  magnitude. 


Power  from  a  Pz  doner.- y- or  Will-  a  Resistive  Load.  n  if. i i 
*  n  "£he  maximum  power  from  a  VI  generator  with  resis¬ 

tive  load  is  given  by: 


^Favq  g3  3 
4  RL 


L  *  2v  f  C, 


Gulton  Industries  Catalog. 
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therefore , 
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rf  ,  2 
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Or 


W  = 


~8  d33  y 33  0  V, 


pz 
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where : 


W  -  power  delivered  to  load  (watts) (for  a  single  disc) 

—  -  - 

f  =  ■  frequency  (hertz) 

R  =  load  resistance  (ohms) 

L 

t  -  thickness  of  the  disc  (meters) 

2 

A  =  area  of  the  disc  (meters  ) 
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SUPPLEMENTARY  REPORT  TO  PIFR-355 


Report  PIPIO  355,  issued  January  1972,  presented  data  on  the 
design,  fabrication,  and  testing  of  a  piezoelectric  generator  intended 
to  be  driven  by  a  Stirling  engine  developed  by  McDonnell- Douglas 
A s trona utics  Company , 

The  following  supplementary  information  is  submitted  relating 
to  tests  on  PZT  511  and  LT7-  2H  generator  materials. 

1.  A  description  of  the  effort  to  build  and  debug  a 
hydraulic  driving  system  for  the  generator. 

2.  A  discussion  of  temperature  measurement  of  tile 
PZ.  generator. 

3.  A  cal.  ulntion  of  theoretical  power  generated  for 
•  •on i pa r i son  with  actual  power  measured. 

■1.  A  i  onsecul  i  ve  tabulation  oi  the  lest  data. 

5.  A  calculation  of  circuit  constant:-,  n.e.  ,  series 
inductance  and  load  resistance)  required  t<> 
match  the  impedance  of  the  PZ.  I  511  generator. 

<  The  <  ircuil  constants  for  the  1,  !  /.  1  generator 
were  determined  experimentally,  set:  tabulation 
data. i 

i).  A  list  of  expenditures  a<  cording  to  task  and  date. 


PZ  GENERATOR  f  DRIVING  SYSTEM 


In  order  to  provide  data  to  establish  the  feasibility  of  a  high 
power /density  generator  it  was  necessary  to  perform  tests  upon  a 
stack  of  piezoelectric  discs  to  determine  its  performance  at  high 
driving  pressures  and  to  validate  the  extrapolation  that  has  been  made 
from  limited  available  data. 

Since  the  proposed  generator  concept  employed  hydraulic  coupl¬ 
ing  between  the  prime  mover  and  the  generator,  a  system  was  designed 
to  provide  an  approximately  sinusoidal  hydraulic  output  for  application 
to  an  existing  PZ  housing. 

A  rotary  valve  was  designed  to  connect  a  hydraulic  pressure 
source  with  the  generator  twice  for  each  revolution  of  the  driveshaft. 
The  design  work  was  performed  with  P.  I.  funds  before  the  start  of  the 
contract. 

Fabrication  of  the  rotary  valve  started  about  July  1,  1971.  The 
rotary  valve,  the  generator,  and  the  hydraulic  system  were  assembled 
during  the  first  3  weeks  of  August  1971.  The  PZT  5H  stack  and  most 
of  the  components  for  the  hydraulic  driving  system  were  available  in- 
house.  Testing  of  the  system  began  on  August  18,  1971. 

The  rotary  valve  was  a  first  experimental  model  and  the  subject 
of  some  development  itself.  The  pump  output  required  was  determined 
by  calculating  the  volumetric  displacement  necessary  to  compress  the 
stack  plus  the  volumetric  compliance  of  the  hydraulic  fluid.  However, 
it  was  found  that  hydraulic  compliance  exceeded  estimates  and  leakages 
occurred  in  the  rotary  valve.  Consequently  the  output  of  the  pump  was 
inadequate  to  raise  pressure  to  the  desired  value.  The  driving  system 
was  modified  several  times  in  order  to  increase  the  driving  pressure 


no 


attainable.  First  the  pressure  was  controlled  by  a  1/4-in.  direct 
act  inn  relief  valve  (Manatrol  1 /4-RA400-4),  which  was  found  to  be  of 
inadequate  capacity  and  sensitivity  to  control  the  pressure  accurately, 
i!  was  replaced  with  a  pilot  operated  relief  valve  /Parker 

Hannifin  MRFN-06-P1AA)  which  was  found  to  be  adequate.  However 
the  output  from  the  pump  (Vickers  axial  piston)  was  limited  when  driven 
at  1800  rpm  and  consequently  drive  speed  was  increased  to  1000  rpm. 
Foaming  of  the  hydraulic  fluid  and  excessive  noise  indicated  that  the 
pump  inlet  was  restricted  causing  cavitation.  Therefore  tint  inlet  tube 
was  increased  to  l /2-in.  diameter  and  the  reservoir  was  elevated  to 
about  2- ft.  above  the  pump. 

The  needle  valve  and  a0.062-in.  diameter  passage  between  the 
rotary  valve  and  generator  restricted  the  flow  and  caused  heating 
locally.  The  valve  was  replaced  with  a  ball  valve  and  the  passage  was 
drilled  out  to  0.  125-in.  diameter.  The  passage  was  in  a  brass  block 
ton  hand)  employed  to  hold  the  Kistler  pressure  transducer.  Evidently 
the  block  was  weakened  by  the  larger  drilling  or  was  damaged  during 
installation.  Consequently  it  broke  during  one  of  the  first  runs  and 
was  replaced  with  a  steel  part. 

The  driving  system  was  limited  in  the  pressure  delivery  capabi¬ 
lity  to  approximately  1700  psi  (7500  psi  on  stack)  and  most  of  the 
hydraulic  power  was  dissipated  as  heat. 

The  Ins;  three  runs  on  the  LTZ  1  generator  were  performed  with 
a  eas /hydraulic  accumulator  m  the  system  to  boost  the  pressure  capabi¬ 
lity  of  the  driving  system  lor  short  durations. 

It  was  found  that  friction  in  the  rotary  valve  made  the  valve  hard 
to  start  at  pressures  above  2000  psi  despite  having  balanced  pressure 

in  >  r  t  s . 
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During  the  test  program  it  was  found  that  much  more  effort  than 
expected  was  required  to  develop  the  driving  system.  Consequently 
this  diversion  reduced  the  amount  of  testing  and  analysis  that  could  he 
performed  within  the  limited  budget. 


if  more  testing  of  PZ  stacks  with  a  hydraulic  driver  is  planned, 
more  effort  will  he  required  to  refine  the  hydraulic  driving  system. 


TEMPERATUR E MEASUREMENT 


In  the  initial  experiment  and  modification  phase  of  testing  the 
PZT  5H  generator  it  was  found  that  the  insulating  oil  temperature 
rose  to  06° F  after  a  few  minutes  running  time.  (The  temperature  was 
measured  with  a  mercury-glass  thermometer.  )  A  cooling  coil,  2-in. 
in  diameter  with  26  turns  of  1.8-in.  O.  D.  copper  tubing,  was  inserted 
in  the  oil.  With  water  running  through  the  coil  the  temperature  of  the 
oil  increased  from  75°F  to  76°F  with  7  minutes  running  time  on  the 
generator. 


In  six  days  of  testing  and  modifying  the  apparatus  after  the  cool¬ 
ing  coil  was  installed  it  was  found  that  the  temperature  variation  of 

o  o 

the  oil  did  not  exceed  7  F  (71-78  FI  or  approximately  the  same  as  the 
room  temperature  variation.  The  cooling  coil  was  led  through  the 
same  hole  that  was  used  for  the  thermometer.  Consequently  the  fit 
around  the  thermometer  was  snug.  The  thermometer  was  removed 
after  the  first  power  run  to  avoid  breakage  due  to  vibration. 

It  was  noted  that  the  insulating  oil  temperature  does  not  repre¬ 
sent  the  bulk  temperature  of  the  PZ  stack.  However  (lie  effort  required 
to  place  thermocouples  in  special  PZ  dis<  s  and  incorporate  them  into 
l  |m:  slack  v.'.is  not  allowed  for  in  the  budget. 


POWER  CALCULA  T I  ON 


The 


>e  expression  for  power  from  a  piezoelectric  8e„era,or  as 
^rtved  in  the  appendix  of  the  report  is  -  t 


W  '  ’  r  d33  *33  ••’2  V 


pz 


where  W 
f 


d. 


S3 


«.n 


V 


r»z 


power  delivered  to  the  load  in  watts 
frequency  of  applied  stress  (hertz) 

piezoelectric  voltage  coefficient  (volts  meters/newton, 
piezoelectric  charge  coefficient  (coulombs/meter2) 
applied  stress  fluctuation  (newton/meter2) 
volume  of  piezoelectric  material  (meters3). 


It  is  assumed  that  the  generator  is  matched  for  maximum  power 

transfer  f-.t  the.  load  (see  page  18). 


The  volume  of  the  PZ T  5H  stack 


is 


v  ox  )  l8'  x  I(i*  80  CC  or  8  X  10-5  m3 


sing  values  from  the  Vernitron  catalog: 

W  .sot  x  Id' 12)  <  I 7  x  10“3)  a2  ,8  x  10'5) 


2.  2  s  10 


!  -5 


T  for  T  in  newton /r»* 


It..  2  >;  it)  l  (6805 1  .7  2  v.  ult.s 


■6  2 


vV  2.05  x  10  a  watts  (for  a  in  psi). 
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This  agrees  very  well  with  the  projected  curve  in  Figure  3. 

- Similarly  the  power  for  the  LTZ  1  genera toi^us: 

W  -  (284  x  io"12)  (29.  2  x  10“3)  (6395)2  aZ  (1.26  x  10‘5) 

-  fl  ? 

W  =  0.  117  x  10  a  watts  (or  in  psi) 

The  actual  power  from  the  L.TZ  1  generator  falls  below  this  theo¬ 
retical  relationship.  However,  this  generator  was  made  up  of  ten  1.25- 
in.  diameter  discs,  0.040-in.  thick  and  one  0.  125-in.  thick  disc  on  each 
end.  Ignoring  the  two  thick  end  discs  the  power  is 

W  =  0.074  x  10“6  cZ  watts. 

This  expression  falls  on  the  2nd  curve  from  the  top  in  Figure  4.  Since 
several  points  fall  between  the  two  curves  it  is  apparent  that  the  effective 
volume  of  PZ  material  lies  intermediate  the  full  stack  volume  and  the  full 
volume  less  the  thick  end  discs.  The  thick  discs  were  included  to  reduce 
the  end  effects  tin  the  short  10  disc  stack  (i.e.  ,  shear  stresses  on  the 
ends  generated  by  the  steel  end  pieces  and  the  Poisson  effect).  If  a 
similar  stack  were  tested  again  it  is  recommended  that  the  end  discs  be 
shorted  and  left  out  of  the  generator  circuit. 

It  must  also  be  considered  that  the  series  inductance  was  not  truly 
optimized.  However  it  is  estimated  that  the  increase  in  power  with 
optimum  indxictance  would  be  less  than  2nk  neglecting  re-optimization  of 
the  load  resistor  to  suit  the  larger  inductor  (see  page  11). 
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Electrical  power  output,  watts 


PZT  5H  GENERATOR 


Power  Huns 


Load  resistcnee:  426Q 

Series  inductance:  1.4  Henr.es 

Ap  -  hydraulic  pressure  fluctuation,  psi 

7  -  stress  fluctuation  =  4.4  ;\p 

V,  -  peak  to  peak  voltage  across  resistor 

P  =  power,  watts 
Ambient  temperature  72  -  73°F 

/  P  VL  P 

September  2,  1971 
Hurt  No.  1 


450 

1980 

130 

5.  C 

600 

2640 

170 

8.  5 

750 

3300 

200 

11.8 

900 

3960 

250 

18.  5 

1050 

4600 

300 

26.  5 

1140 

5000 

320 

30.  0  - 

temp.  78° F 

September  3,  1971 
Hun  No.  2 


450 

2000 

120 

4.25 

600 

2650 

150 

6.6 

750 

3  300 

190 

1 0.  6 

1  150 

4950 

300 

26.  5 

975 

4300 

275 

22.  2 

l-*.6 


Aa 

VL 

P 

|V  •  September  7,  1971 

|  Run  No.  3 

j  750 

j; 

!  1150 

1 

3300 

180 

9.  5 

4950 

235 

16.2 

1  1500 

j 

6600 

310 

28.4 

|  1875 

8300 

340 

34.  0 

[ 

!  September  8,  1971 

|  Run  No.  4 

(  750 

3300 

220 

14.  3 

1150 

4950 

280 

23.2 

1350 

6000 

300 

26.  5 

1500 

6600 

305 

27.  5 

1650 

7250 

300 

26.5 

i  September  9,  1971 

Run  No.  5 

600 

2640 

125 

7.  3 

1000 

4400 

180 

9.  5 

1500 

6600 

230 

15.4 

1900 

8400 

260 

20.  0 

Repoled  stack  to  50  V/mil 
Run  No.  6 

.DC  pulse. 

250  ms 

pulse  at  2  Hz. 

600 

2640 

150 

6.  5 

1000 

4400 

260 

20.0 

1500 

6600 

320 

30.  0 

2100 

9300 

300 

26.  5 

1000 

4400 

210 

13.0 
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LTZ  1  GENERATOR 


fioloction  of  I.o.td  Roe  is  tor 

Series  inductance  =  18  H 
Stress  fluctuation  =  5500  psi 


October  7,  V) , 


7.0 

350 

2.2 

8.  5 

400 

2,35  4-max 

0.  5 

420 

2.  32 

10.0 

430 

2.  31 

11.  5 

460 

2.30 

Selection  of  Series  Inductance 

Rl  -  8.5  K n 

Stress  fluctuation  -  5500  psi 


Power  Runs 


Load  resistance  R  =  8.5  K Q 

JLi 

Series  Inductance  L  =21  Henries 


Ap  =  hydraulic  pressure  fluctuation,  psi 
to  =  stresB  fluctuation  =  4.4  Ap 

--  peak  to  peak  voltage  across  resistor 
P  =  Power,  watts 

Estimated  running  time:  2  min,  per  data  point. 

October  7,  1971 
Run  No.  1 


ip 

Aa 

ZL 

P 

1250 

5500 

415 

2.  52 

1560 

6850 

500 

3.68 

1600 

7000 

520 

3  96 

935 

4100 

330 

1.60 

625 

2750 

250 

0.92 

312 

1375 

120 

0.21 

75(1 

3300 

280 

1.  15 

500 

2200 

200 

0.  59 

500 

2200 

200 

0.  59 

37  1 

LoO 

160 

0.25 

625 

2750 

250 

0.  92 

750 

3300 

290 

1. 24 

Ji£ 

6  a 

>1 

P 

October  3,  1971 

Run  No.  2 

623 

2750 

240 

0.  85 

935 

4100 

350 

1.S0 

October  8,  1971 

Run  No.  3 

1000 

4400 

380 

2.  1 

500 

2200 

220 

0.  71 

1300 

5700 

490 

3.  5 

October  11,  1971 

Run  No.  4 

600 

2650 

220 

0.71 

800 

3500 

290 

1.25 

1000 

4400 

360 

1.90 

1500 

6600 

540 

4.  30 

1700 

7500 

600 

O 

CO 

9 

U3 

October  13,  1971 

Rim  No.  4 A 

1400 

6200 

330 

Odd  pressure  wave 
with  2  peaks 

1400 

6  200 

470 

3.25  Normal  pressure  wave 

1500 

6600 

510 

3.85  Pressure  falling  with 
time,  could  not  read 

1700 

7500 

530 

4.15  scope  accurately. 

Ap  Aa 


October  14,  1971 
Run  No.  5 


1100 

4850 

1200 

5300 

1300 

5700 

1500 

6600 

1800 

7900 

1800 

7900 

2000 

2800 

1800 

7900 

1000 

4400 

1500 

6600 

2400 

10600 

2600 

11400 

1500 

6600 

P 

350 

1.80 

400 

2.  35 

440 

2.85 

510 

3.85 

630 

5.85 

640 

6.  00 

460 

3.  10 

Erratic  behavior. 

620 

5.65 

H  It 

440 

2.85 

H  1  • 

360 

1.9 

H  U 

360 

1.9 

1?  <» 

560 

4.6 

1!  1* 

580 

4.95 

1 1  •? 

350 

1.8 

?»  U 

*These  stress  figures  are  in  doubt  because  it  is  suspected  that 
the  nO"-ring  extruded  and  jammed  the  piston. 


Ap  A  a  VL  P 


October  15,  1971 
Run  No.  6 
10:00  A.M. 


1500 

6600 

530 

4.  15 

1400 

6200 

500 

3.  7 

1900 

8400 

640 

6.  0 

270C 

12000'' 

650 

6.  2 

Erratic  b 

ehavior 

2700 

12000 

670 

. 

6. 6 

9  ! 

it 

1900 

8400 

600 

5.3 

tl 

?  t 

2100 

9250 

>  *  660 

6.4 

1! 

It 

2700 

12000 

660 

6.  4 

I« 

tl 

2600 

11500 

620 

5.7 

1? 

H 

2600 

11500  j 

670 

6.6 

15 

It 

October  15,  1971 
Run  No.  6 
2:30  F.M. 


2300 

10100 

680 

6.8 

Erratic  behavior 

2200 

9700 

>  *  670 

6.6 

!J  M 

2  400 

10600 

) 

!  690 

7.  0 

it  it 

2200 

9700 

690 

7.0 

*These  stress  figures  are  in  doubt  because  it  is  suspected  that 
the  "0"-ring  extruded  and  jammed  the  piston. 
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IMPEDANCE  MATCHING 


A  knowledge  of  the  behavior  of  piezoelectric  materials  was  used 
to  postulate  the  equivalent  circuit  (shown  in  Figure  A  1 )  of  a  PZ  gen¬ 
erator.  The  voltage  source  shown  in  the  equivalent  circuit  represents 
an  ideal  voltage  generator.  Since  a  PZ  stack  is  capacitive,  a  capacitor, 

,  has  been  included  in  the  equivalent  circuit.  The  resistor,  R  .  , 
represents  the  internal  losses  of  the  PZ  generator  while  Rp  is  the  shunt 
resistance  present  in  all  PZ  stacks.  This  shunt  resistance  causes  an 
initial  voltage  on  a  stack  to  decay  and  is  included  in  the  circuit.  The 
series  capacitance  Cg  is  not  equal  to  the  capacitance  of  the  stack,  but 
is  calculated  as  shown  below. 

If  tiie  value  of  the  load  resistor  on  the  output  is  small  compared  to 
Rp,  the  circuit  can  be  simplified  to  that  shown  in  Figure  A2.  !f  it  is 
also  assumed  that  the  impedance  due  to  and  are  small  compared 
to  R  ,,,  then  the  generated  voltage,  V,..,  can  be  assumed  equal  to  the  open 

F  vj 

circuit  voltage  measured  at  the  output  of  the  generator.  Assuming  the 
generator  output  is  sinusoidal,  the  circuit  equation  in  the  frequency 
domain  is 

V'  -  l  «R.  4  Rc  f  Z„)  ICq.  (1) 

VI  La  »->  1/ 

or 

voc  *  57  ,rl  '  Rs  *  zc>  E<>-  121 

Li 

[f  the  generated  voltage,  V„,  is  independent  of  frequency  then,  by 

1.  A  circuit  similar  to  this  is  shown  on  page  17  in  Piezoelectric 
Technology  Data  for  Designers,  by  the  Piezoelectric  Division 
of  Clevite  Corp. 


assumption  VQC  is  also.  Hence  the  output  at  two  frequencies  can  be 
equated: 


Therefore  R 


or 


Substituting: 


0, 


R, 


O, 


fRL  +  RS  f  zc> 


R, 


f  RS  '  ZC» 


Eq.  (3) 


is 


R. 


O 


V  4  - 

°2  RL 


2  , 


V 


O, 


(zc) 


O,  R 


L 


(zc) 


V„  - 
°1  °2 
R, 


Eq.  (4) 


o. 

c 


o 


V<VV°1  '  RL,*ZC  ~  RU*l  CJ 


R, 


V 


o, 


o.. 


Eq.  (5) 


RS 


K  -  v°.) 


'V 


O- 


o 


R, 


iV~ 


■1 


1  1  JL 

C 


V  -  V 
°1  °2 


Eq.  (61 


RS 


?3  o  _  5  ?)  100  t  f _I 

•  *  V  197  377  /  C 


(r>.  7-3.91 


Eq.  (7) 


r  — - _  ioo 

77* 


iq.  ( 8 i 


I  <34 


Using  Eq.  (8)  and  the  experimental  values  for  the  PZT  5H  gen 
orator  shown  in  Table  3A  for  Eq.  (2)  yields: 


60.  i 


5.7 

100 


100  t 


2.61  x  10 
C 


100  t 


>77  C 


Eq.  (9) 


or 


C 


5.  7 

"07 


(5.26)  10"6 


5.  01  n  ( 


Eq.  ()()) 


Therefore: 


~3 

R  -  .?.«..61x..LQ  _  100  =  421  n  Eq.  (11) 

15  5.01  x  10"J 

The  capacitive  impedance  is: 

Z  =  377  (5.  01  xTF6  5320  Eq*  (12) 


For  maximum  power  transfer  from  the  generator,  the  load  resis¬ 
tance  must  be  equal  to  the  internal  resistance,  Rg,  and  the  capacitive 
reactance  must  be  balanced  by  the  addition  of  an  inductive  reactance  of 
the  same  magnitude.  Therefore: 


L 


60  Hr. 


532 

377 


1.4)  h 


Eq.  (13) 


It  is  necessary  that  the  voltage  measurements  used  to  calculate 
the  equivalent  circuit  be  accurate  to  prevent  errors.  For  instance. 


if 


Vq,  was  measured  5°/<-  high  (6.0  volts  instead  of  5.7)  and  Vq2  was  mea 

sured  s  !,  low  (3.  7  volts  instead  of  3.9)  the  series  capacitance,  C  , 

S 

increases  I  0'!' . 


R 


S 


(3.  7  -  6.0)  100  t 
_ _ 19  < 

67cT~.  7 


6.0  1 
377  C 


2. 87  x  10 
C 


-  loon 


Therefore: 


C  5.  52  x  10"6  f 


Although  the  induction  calculation  above  produced  a  value  that 
worked  very  well  in  tests,  later  calculations  on  different  stacks  pro¬ 
duced  values  that  were  not  realistic.  Another  method  was  used  to 
determine  the  circuit  values  which  worked  fairly  well  in  actual  tests. 
However  furhter  work  is  necessary  to  characterize  PZ  generators.  A 
vital  part  of  this  is  to  determine  whai  relationship  the  stack  capaci¬ 
tance  has  to  the  equivalent  series  capacitance. 
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Appendix  B 

INDEPENDENT  DEVELOPMENT  OF  PIEZOELECTRIC 
GENERATOR  TECHNOLOGY  USlN^t' CYCLE 
AVERAGED  FORMALISM 

This  appendix  discusses  the  analytical  closed  solution  formalism 
associated  with  continuous  and  pulsed  operation  of  the  piezoelectric 
generator.  For  both  modes  of  operation,  two  cases  are  considered: 
(1)  the  purely  resistive  load,  and  (2)  the  load  with  inductance  added  to 
compensate  for  the  capacitance  of  the  piezoelectric  stack.  For  con¬ 
tinuous  operation,  an  analytical  correlation  with  the  P/.T-SH  data 
presented  in  Appendix  A  is  also  made. 


B,  1  CASE  I  -  CONTINUOUS- POWER  OPERATION  WITH  RESISTIVE 
LOAD 

This  case  is  shown  in  Figure  B- 1  where  the  following  conditions  apply. 
Mechanical  Input:  Sinusoidal  Compression 
Load:  Resistive:  no  matching  inductance 


2238 


Figure  B-1.  Mechanical  and  Electrical  Equivalent  Circuits 
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'I  in-  shunt  resistor,  R  ^  for  crystals  of  useful  sixes  is  on  the  order 
of  10  ' ohms;  C  j  is  the  capacitance  of  the  leads  and  load.  For  this 
case,  Ch  is  very  small,  e.g.,  a  few  picofarads;  therefore,  X(;l  is 
also  very  large  (  1 0  1  ^  ohms).  The  electrical  equivalent  circuit  can 
then"  be  simplified  as  in  Figure  B-2. 

With  sinusoidal  compression  as  the  mechanical  input,  the  instantaneous 
v oil  age  is 


v  -  V  sin  wT 
m 

and  at  steady  state,  the  instantaneous  current  is 


i  I  .  iti  ( —jT* o) 
m . 

Where  the  peak  current  1  ,  and  peak  voltage,  V  are 

related  by 


V  V 


with 

R.  Load  resistance 

R  Internal  stack  series  resistance 

Stack  capacitive  reactance 
/.  :  Circuit  impedance 

I  * 


2?S9 


*C  '  u»c0 


Figure  B-2.  Simplified  Electrical  Equivalent  Circuit 
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tnd  the-  phase  angle,  p,  is  given  by 


O  -  tan 


X 

-1  c 


hl+rs 


Straightforward  mathematical  manipulation  leads  to  the  following 
expression  for  power  W,  delivered  to  the  load 


V  I  Rf 

nr  ■  m  rn  _ L 

i  /. 


By  definition  and  t rom  typical  piezoceramic  property  data. 


tan  6  «1 


The  maximum  power  transferred,  therefor.*,  occurs  very  nearly 
at  Xc,  For  this  case 

V  2  V  - 
nr  -  m  .  rms 


where 


rms 

As  was  discussed  in  Section  2.  1.3 
V  (F  /A\  gt 

rms  ^  rms  }  * 


The  re l‘o  re 


B.2  CASK  il  -  CONTINUOUS- POWER  OPERATION  WITH  MATCHING 
INDUCT  ant:  K 

Mo  t.- ha  Eli  cal  Input:  Sinusoidal  Compression 
Load:  Resistive  with  Matching  Inductance 

r  lie  equivalent  electric  circuit  for  the  ease  is  shown  in  Figure  R-i. 


When  Nj,  X(-,  this  circuit  reduces  to  the  form  of  Figure  R-4, 
where  X[,  ujL  inductive  reactance. 


Figure  B-3.  Simplified  Electrical  Equivalent  Circuit  with  Inductance 
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2242 


Figure  B-5.  Mechanical  and  Electrical  Equivalent  Circuits 


Current  a,  function  of  time  after  the  switch  is  closed  is 


_V 

R 


^  t-T|RsMiL)C0J 


This  leads  to  an  average  eurre 


em  over  a  half  cycle  of 


1 


IavH  ' T7T,-  f 


l/Zf 


_i_  r,/2f_Voc 
l/Zf  J  R  ‘R 


l/2f  /  idT"  1777  /  -7~-  <-t/(r»’Rl.c4 


S  L 

integration  and  evaluation  at  the  limits  leads  to 
^avy  '  ^ •  ^oc^o  |  1-  t  ~l/|(2r  C0  H  R  )|| 
The  average  power.  W,  delivered  to  the  load 


(IT 


can  now  be  determined 


W~  Iavs"RL  "-  Uf  VocCo)2rl  Ji_  2e -1/2Ft4  t  -1/FtJ 


whe 


re 


7  '  ,RL+l!s»Co 


iH  L  III  & 


R|.  to  pnrfucc  ie  meximSm^wel  i  ’1*°  "  v.lo 

following  arEumem:  When  is  v^rt  sms  ,  '  -"oostrated  by  the 

*'  r>  small,  the  quantity  outside 
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the  brackets  is  very  small;  when  R^  is  very  large,  the  quantity  inside 
the  brackets  is  very  small. 

_ _  _  —  • 

B.  4  CASE  IV  -  PULSED-POWER  OPERATION  WITH  INDUCTANCE 
This  case  is  shown  in  Figure  B-6  with  the  following  conditions. 

Mechanical  Input:  Periodic  Compression,  frequency  w 
Load:  Resistive  plus  Inductance 

This  case  is  a  little  more  difficult  mathematically.  The  solution 
lor  instantaneous  current  depends  on  the  relative  values  r.f 


(Rt  +R L.  and  C 
L>  S 

Case  IV -A: 


fJL 

\ZL 


1  - 


^or 

UJ  * 


Normal  Damping 
-RT/2L  <  6T  _t  '&T 
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Figure  B-6. 


Mechanical  and  Electrical  Equivalent  Circuits 


— r-  •  r  ,  -  r 

L  s 


(R/2JA2  -  l/LC 


Case  IV- B: 


ffi 


2  I 


’  Critical  Damping 


i'  r  c-kt/zl 


Case  IV-C: 


R  y-  1 

Tf  I  ”[  Oscillatory 


4l -e  RT/^D  sin  nT 


whe  v 


2]  0.  5 


As  in  Case  III,  average  currents  can  be  calculated  inr  Case  IV- 
A.  IV.  anri  : 


Again,  average  power  delivered  t<*  the  load  is 


\V  I“avg  K  i  _ 


Inr  a  given  stack  >  <a  pac  itance  >,  an  R  and  L  van  be  found  t< 


t  t  ] 

produce  tlie  maximum  power.  Again.  R  would  normally  be 
increased  above  this  value-  for  the  Ijest  tradeoff  between  power 
dens  i  t  v  a  nd  «•  tf  i.  ie  nev  . 


15.  CORRELATION  OF  EXPERIMENTAL  DATA  WITH  THEORY 

Calculations  were  performed  utilizing  the  formalism  for  continuous 
operation  presented  in  Cases  I  and  II.  I  he  particular  example  given 
he  i  e  is  tor  a  PA  I  -all  stack  consisting  ot  100  discs  connected  electri¬ 
cally  in  parallel.  Each  disc  is  0.  040  in.  thick  x  1 . 2r>  in.  din.  Such 
a  stack  was  used  in  the  experiments  reported  in  Appendix  A. 

'•  igut'e  B-  /  presents  the  results  of  applying  this  formalism  to  the 
operation  conditions  ot  Run  No.  1,  B-2-7  1,  dose r  third  in  Appendix  A. 
Also  shown  are  the  experimental  data  points  fathered  on  that  run. 

The  calculations  ot  P/.T-aH  stack  performance  make  use  of 
materials  parameters  provided  in  Cieviie  Catalog  No.  695. 

Generally,  low  signal  parameters  were  used  with  some  corrections 
for  high  stress  and  high  field  made  to  the  dielectric  constant  and 
dissipation  factor.  Data  for  these  corrections  are  only  strictly  valid 
to  2000  psi,  and  these  corrections  should  he  considered  qualitative 
at  higher  stresses.  This  is  especially  true  because  a  linear  depen¬ 
dence  of  each  correction  was  assumed  for  extrapolation  which  is  not 
necessarily  correct.  Series  resistances  *  internal,  lead,  and  inductor) 
are  *1  is  regarded  in  all  powe  r -calculations  using  this  formalism.  Even 
though  there  are  these  limitations  on  the  annlvsis,  some  correlations 
and  dependencies  are  worth  examination.  In  Figure  B-7.  calculated 
power  exhibits  the  same  general  dependence  on  applied  stress  as  the 
experimental  data.  These  calculated  values  lend  to  be  higher  than 
the  data.  When  a  low  field  relative  dielectric  constant  of  5200  is 
used  instead  of  5400,  this  difference  disappears.  Such  a  difference 
in  dielectric  constant  could  result  from  the  temperature  of  the 
material  or  time  after  poling. 

Figure  15-8  shows  the  calculated  dependence  of  the  P/.T-5H  stack 
output  power  mi  load  resistance  at  a  given  stress  fluctuation 
i  2  b  4  0  psi).  With  respect  to  power,  the  calculated  optimum  circuit 
resistance  of  about  700  ohms  with  the  given  inductance  ( X  -Xj  -  7O0S2  ) 
and  the  experimental  optimum  load  resistance  of  427  ohms  imply  as 
much  as  27;;  ohms  ot  series  resistance  elsewhere  in  the  circuit,  if 
the  calculations  are  correct.  Such  resistance  would  probably  bo 
primarily  in  the  stack  with  a  small  amount  in  the  inductor  and  leads. 
Increasing  stack  temperature  leads  to  an  almost  exponential  rise  in 
tan  h  and  increases  the  of  fee  live  series  resistance.  Such  a  rise  with¬ 
out  cooling  could  lead  to  rapid  stack  deterioration:  any  temperature 
increase  causes  more  self  heating  and  more  temperature  rise. 

This  same  analytical  in'  i.ialism  also  provided  good  agreement  for 
power  vs  stress  and  power  vs  load  for  ciata  gathered  at  DWDL. 
However,  with  both  Physics  International  and  DWDL  data,  the  depen¬ 
dence  on  inductance  was  in  error.  While  the  power  versus  inductance 
dependence  was  qualitatively  correct,  the  optimum  inductance  calcu¬ 
lated  was  ove  r  A ()''.>  higher  titan  measured.  By  reformulating  the 
equations  on  an  instantaneous  rather  than  time -averaged  basis,  and 
accounting  lor  the  charge  on  the  stack  throughout  a  given  cycle,  this 
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OUTPUT  POWER  <W> 


Figure  B-7.  PZT-5H  Stack  Output  Power  Dependence  Upon  Stress 


2240 


lb,R 


discrepancy  was  eliminated.  This  latter  formulation  is  described  in 
more  detail  in  Appendix  C  and  applied  to  Physics  International  and 
DW'DL  data. 


A  sail. pie  set  of  calculations  tor  obtaining  the  curve  in  Fipare  B-7  is 
”iven  in  Section  B.  6.  Symbols  used  are  defined  in  Table  H  -  1 . 

B.  6  SAMPLE  CALCULATION 
P/.T-SIl  Stack: 


100  discs,  0.040  in.  thick  x  1  .25  in.  dia  {connected  in  parallel) 
D  =•  1.2 a  in.  -  0.03  174  in 
A  D"  1.227  in."  0.  0007<*2  m“ 

t  •  0.  040  in.  0.  00  10  16  m 

1’Z  volume  /disc  ■  At  0.  049  1  in.  0.  H04  x  10  m 

P/.  stack  volume  -  100  At  4.91  si-.’  0.804  x  10  m^ 

t- 1  it  (8.1 -lax  10*  1 t  1400  fl  -  0.  04  fA<r)/2000]  5 

3  3  o  r  L  J 

Assume  A  o  -  1 9 8 0  psi 
t  |  0.  28r>  x  ! 0”  ' 

C  f  ^  A/t  0.224  x  10  *  fa  rati /disc 
o  3  3 

Stack  capacitance  {lOOh  0.224  x  10*'}  0.22a  x  10  farads 

x  — 1 — -----  {/  (2  iroOit  0.  224  x  I9'5)  ILLS  Si 

c  2  trtC. 

I  J 

Assume  R,  42  4  <2.  L  L  3  henries 


S' 

L 

!  2tt  t;  0  It  ! .  it 

:  X  - 

X.  1  1.8812 
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/  'J 

i . 
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.  ie.i 
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■xeept  the  load  is  dis  re->a  rded. 

0.  *26 

•  -r,i  rd  i  n.t,  ail  i  irnii!  resistance  except  the  load) 
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Table  B- 1 

SYMBOLS  USED  IN  SAMPLE  CALCULATION 


Symbol 


t 

D 

f 

E 

F 

Vol 
( 


}  1 


R 


L 

T 

'll 

F 

P 

F 

PP 

F :  \ 

r  ms 

Fj 

?.v  g 

/. 

X 

c 

C 

o 

I 

P 

L 

xl  ; 

tan  6 


Y 

i 

Y 

t 

Y 


rms 


(X* 


\- 


pp 

i. 


it 


Definition 


I 


Area 

Thickness 
Diameter 
Freqne  ncy 

Electric  field,  iv/mi 

I'n  rce 

Volume 

Piezoelectric  coefficient  t volts/applied  stress) 
Load  resistance 

Relative  dielectric  constant >  (free) 

Peak  force 
Peak- to- peak  force 
Rms  force 
Average  force 
Impedance 

Capacitive  reactance 

Capacitance 

Current 

Powe  r 

Inductance 

Inductive  reactance 

Dissipation  factor  effective/ 

(X  .  teffecrive  5,1 

c  series’ 

Rms  voltage 
Open  circuit  voltage 
Pcak-to-peak  voltage 
Voltage  across  the  load 
!  lea  t  gene  rated 

Efficiency  _ 
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£o<MXtftil>v  cK/>ra^>v'u«wcici«jiii 


(F/'^MAX  3980  psi  (198G)(6895)  0.  1365  x  lO^n/m2 

I F/A ?avg  990  psi  -  0.683  x  10  n(.m1' 

(F/A), 


(F/A), 


'RMS 


vTT 


700  psi  0.483  x  10^  n/m2 


V RMS  ~  ^  {0.483  x  loS  (0.  0!97)(0.  001106)  ^  96.  6 

V  =  VD.  .  ^  96.6  v 

oc  RMS 

V.  V  cosO  =  (96.  6)  (.  527)  ~  50.  9  v 
1.  oc 

V  -  >y/T~  V.  -  14  3  v 

PP  L 

V  nr  , 

I  - — r~  ■■  rrfwT  ;  0.  1195  amp 


/. 


808 


V  I  (50.  9H0.  1 195)  -  607  w 


Tan  6  =  ( Tan  6  ) 


A  ( ta  n  6  ) 


( R 


low  signal  A  E 

0.  02  ~  (0.04  -  0.02)  (0.483  x  10?)  (0.0197)  -  0.0834 
_  _____  ■  {0 .3)  — 

series'  .  -  X  tan  6  -  ( 1 178  ft)  (.  0834)  -  98  Q 

eff  c.  i 


H  -  I2(R  .  )  ,r  (0.  1  195)**  (98)  ^  1.  40  w 

series  eff  • 


9  - 


6.07 


P  -  'H  6.  07  -  1.40 


5)2 

-  0.81  -  81% 


This  value  of  efficiency  disregards  all  mechanical  losses,  lead,  and 
inductor  losses. 


-° 

IN xegration  formalism  bMERrc 


ICAL 


fcm»E“,  piez'oelectrdc  ™*r*?P  nUmeH.cal  iteration 
developed  at  DWDL.  This  <v  r  e  a  or  anaIysis  which  was 
integrated  with  the  dynamic  rnrL'S"°  per™Us  sV««th  studies  u.-l 
better  correlation  «  ™wZI  ^  *?  P"vM„ 

cribed  in  Appendix  R  Ushitr  th*<>  ***?.• than  the  techniques  dos 

ga  the  red  b/S ics i„te  r 1 l£n5  i  *•  LTZ-  1  data 

analyzed.  Data  gathered  at  DWDI*!*  r,ep°rle^  in  Appendix  A  is 
n  discussion  o,  tL 

c.  I  COMPUTER  PROGRAM 

compn.er  ihstnntaheously  Sonitoh  S  Co  t“^t ‘ 
gram  is  «iven  in  Figure  C-l  The  W*  V  •  COmPutcr  Pr°- 

The  basic  charge  balance  on  the  capacitor  is; 

Q  "  charge  generated  .  charge  lost 

-  (capacitance)  (voltage  generated; 

-  (average  current)  (time  increment) 


page  blank 


where  the  dielectric  constant  used  is  field  dependent: 


The  average  current,  I,  is  determined  by  an  iterative  process  of  cal¬ 
culating  voltage  on  the  stack  and  average  impedances  in  the  circuit  at 
each  time  instant  in  the  cycle.  AH  current  values  for  a  cycle  are  then 
used  to  calculate  an  average  current.  Work  delivered  to  the  load  is 
the  average  Irms)  current  squared  times  the  load  resistance.  While 
all  calculations  and  experiments  to  date  Have  Involved  sinusoidally 
applied  force,  the  computer'  program  is  capable  of  utilizing  whatever 
time  dependence  is  input.  For  forces  requiring  more  rapid  changes 
within  a  time  frame  than  the  60  Hz  sine  wave  utilized,  the  only  modi¬ 
fication  necessary  is  reducing  the  computer  time  step  increment 
which  results  in  additional  computer  running  time  for  each  case. 

The  work  input,  calculation  also  utilizes  the  charge  balance  to  deter¬ 
mine  the  instantaneous  compliance  of  the  slack. 


Mechanical  work  input  ~  Sum  of  all  F  (DX)  over  a  cycle 


Figure  C-1.  Piezoelectric  Stack  Equivalent  Circuit 
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Symbol 


Do.  fin  it  ion 


The  net  electrical  efficiency  of  the  generator  is  then  given  by 


work  to  load 


I  Rload 


-  mechanical  work  input  J-  F  { DX) 

Alternately,  the  efficiency  of  the  stack  alone  can  be  determined  by 
including  the  work  dissipated  in  the  inductor,  although  this  definition 
was  not  used  for  anv  of  the  efficiency  numbers  resorted  here. 


('V, 


d 


/  F  i  DXi 


r..l  CORRELATION  WITH  PHYSICS  INTERNATIONAL  LT/.-i  DATA 

The  formulation  developed  in  Section  C.  I  was.  utilized  in  analyzing 
LT/.-l  stack  performance.  Because  LT/.  -  I  mate  rial  is  essentially 
identical  with  P/.T-4.  the  materials  parameters  for  the  latter 
i  Reference  C-P  were  used  in  the  computer  program. 

These  parameters  include: 

(  MOO 


PC  0.  17 

p  ,  p>0. 000  volt /meter 

'  MA  X 

-•  !).  !)F>  I  volt  meter/ newton 


S< •  t  hat  t 
and  tank 


/  17  SVabt\ 

)0  ‘  <>  p»0,  ooo  t  / 

/.)  04  -  O.  t>t>4\  ,,  ,  p 

*  \  s 1  * n , non  /  A 


The  LT/.-I 
} ;  2. in.  <i ia 
«•  illations  we:* 
t  antes,  a  rut  serier 
of  fiO  II/.. 


stack  consisted  of  10  discs  each  0.040  in.  thick  x 
with  all  discs  electrically  connected  in  parallel.  Cal- 
•uii  tor  a  varietv  of  stress  fluctuations,  load  re  sis  - 

inductance.  All  calculations  were  for  a  frequency 


Figure  C.-l  shows  the  results  of  these  power  versus  stress  calcula¬ 
tions  superimposed  on  the  Physics  International  figure  previously  pre 
Rented  Appendix  A.  Particularly  good  agreement  occurs  tor  stress 
fluctuations  at  or  below  f>000  psi.  At  higher  stress  levels,  calcu¬ 
lated  values  arc  mostlv  higher  than  measured  perlormance.  This 
may  result  from  »  and  Tan  b  having  a  nonlinear  field  dependence. 

The  linear  dependence  of  these  two  quantities  was  assumed  based  on 
two  points  of  data  provided  in  Reference  C-l:  the  low  Hold  values 
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Figure  C-2.  DWDL  Analytical  Correlation  with  Physics  International  LTZ-1  Generator  Performance 


i  n 


of  tano  and  f ,  (ho  electric  field  for  tan6  0.  04  ( Ew  A  v  •  ton  onn  v/mi 

and  the  corresponding  change  in  e  for  that  field  (PC  ^0.  17). 

Figure  G- 3  shows  the  dependence  of  output  power  on  series  load 
inductance  for  the  experimental  conditions  where  inductance  was 
varied.  the  agreement  between  experiment  and^Beory  appears 
close.  Experimentally,  output  power  continued  to  increase  grad- 
ually  as  inductance  was  added.  The  largest  inductor  available  was 
Z1  henries.  Calculations  indicate  that  increasing  the  inductance 
further  would  have  resulted  in  a  maximum  output  power  for  an 
inductance  of  over  50  henries  (where  X  effective  X  ).  This 
maximum  of  over  10  watts  at  30  henries  is  over  4  times  the  mea¬ 
sured  value  at  21  henries. 


Figure  C-4  presents  the  dependence  of  output  power  on  load  resis¬ 
tance  tor  both  the  18  henry  experimental  case  and  55 -henry  calcu¬ 
lated  case.  For  the  calculated  case,  lowering  the  load  resistance 
increases  output  power  even  further.  With  all  circuit  reactance 
minimized,  a  load  resistor  equal  to  the  effective  internal  resis¬ 
tance  of  the  stack  produces  maximum  obtainable  output  power  for 
a  given  stress  fluctuation. 

Disadvantages  o!  operating  at  this  increased  power  densitv  include 


density  is  more  important  than  efficiency. 


C.  3  PIEZOELECTRIC  MEASUREMENTS 


C.  3.  1  Test  Equiprr-c nt 

Tests  were  carried  out  by  DWDL  to  establish  practical  operation 
parameters  for  electrical  power  generation  and  to  verify  analytical 
predictions  of  portormance.  Ail  tests  were  performed  using  the 
elect ruhyd  raulic  materials  loading  system  at  Battelle  Pacific  North¬ 
west  Laboratories  in  Richland,  Washington.  This  unit.  Model  312.21 
(  12.  5/25  KIP'.,  was  manufactured  by  MTS  Systems.  Inc.  Loads  were 
applied  sinusoidally  at  60,  90,  dr  120  Hz  with  a  minimum  of  1000  psi 
compression  and  up  to  1  3,  000  psi  maximum  compression.  At  each 
compressive  loading,  the  electrical  load  on  the  stack  was  changed  by 
varying  load  resistance  and  inductance. 

Equipment  utilized  in  these  piezoelectric  tests  experiments  is  shown 
in  Figure  C-5;  however,  none  of  the  MTS  System  control  equipment 
is  shown.  Included  in  the  electrohyd raulic  system  are  control, 
loading,  and  hydraulic  power  subsystems.  The  control  loop  utilized 
a  load  cell  in  series  with  the  test  specimen,'  together  with  the  neces¬ 
sary  electronic  signal  conditioning,  comparator,  and  readout  modules. 
A  pump  maintained  hydraulic  fluid  pressure  in  a  servo  valve  controlled 


by  the  servo  controller.  In  all  tests,  a  sinusoidally  varying  force 
was  maintained  on  the  loading  piston  by  the  servo  valve.  The  piezo¬ 
electric  stack  under  test  was  mounted  oh  top  of  the  piston,  just 
below  the  load  cell,  all  within  the  load  frame  shown. 

i  .'M?' 

~  Mechanical  work  input  to  a  piezoelectric  stack  under  test  was  deter¬ 
mined  as  the  area  within  the  force  versus  displacement  loop  dis¬ 
played  on  a  Tektronix  Model  503  X-Y  oscilloscope.  The  force  signal 
was  obtained  from  the  MTS  system  load  cell.  The  signal  of  1  volt/ 
1000  lb  provided  to  the  oscilloscope  was  displayed  on  the  Y  axis. 
Displacement  was  measured  using  a  Bently  Nevada  Model  302L30 
proximity  measuring  system  consisting  of  a  proximitor  (probe 
excitation  source)  and  probe  (sensing  element).  Displacement  is 
measured  by  the  non-contacting,  eddy  current  probe  with  an  output 
signal  of  the  proximitor  supplied  to  the  X  axis  of  the  oscilloscope. 
V/ith  the  Model  316L  probe  used  in  displacement  measurements,  the 
gain  was  0.20  volt/mil  displacement  (±10%)  according  to  manufac¬ 
turer  specifications.  This  unit  was  capable  of  measurements  over 
0  to  120  Hz  without  significant  phase  lag  or  amplitude  attenuation. 
Displacement  variation  under  static  load  conditions  could  be  veri¬ 
fied  with  a  Bendix  Indi-Ac  A.  T.  Gaging  System  Model  AT- 1  with 
a  Type  T-220  gage  head.  However,  this  unit  with  a  head  operating 
on  the  LVDT  (linear  voltage  differential  transformer)  principle  could 
not  produce  accurate  results  at  operating  frequencies. 

Figure  C-6  shows  a  piezoelectric  stack  on  the  MTS  System  with  the 
Bendix  probe  on  the  left  and  the  Bently  Nevada  probe  on  the  right. 

The  electrical  leads  from  the  piezoelectric  stack  are  also  shown. 
Electrical  power  output  was  measured  as  P  =  V£/RL,  wfaere  VL  was 
the  voltage  across  the  load  resistor  R^.  This  voltage  was  detected 
with  a  Ballantine  Model  355  digital  voltmeter  (±  0.  1%  accuracy  with 
sine  waves;  this  meter  does  not  indicate  true  rms  voltage  for  non- 
sinusoidal  imput).  In  later  tests,  all  load  resistors  were  Dale  wire- 
wourri  1%  power  resistors  each  with  25-w  heat  dissipation  capability. 
Inductors  were  Triad  filter  chokes  Model  C-  /X  (90  ma,  10  h, 

270  ohms).  Overnight  recording  of  stack  output  voltage  was  made  on 
a  Mosely  Model  7100B  strip  chart  recorder. 

|  C.  3.  2  Test  Results 

Initial  testing  was  done  on  a  16-  disc  unbonded  stack  (each  disc 
0.  10  in.  thick  x  0.  5  in.  dia)  consisting  of  unbonded  HDT-31 
(equivalent  to  PZT-4)  material  supplied  by  Gulton  Industries,  Inc. 
Both  stress  level  and  load  resistance  were  varied.  A  stress  fluc¬ 
tuation  of  8000  psi  (1000  to  9000  psi)  produced  4.  25  watts 
( 13.  5  w/in.  3)  with  a  load  resistance  of  175  K  ft.  The  load  resis¬ 
tance  for  maximum  power  decreased  from  180  K  at  5100  psi  to 
170  K  at  8200  psi  stress  fluctuation,  suggesting  a  stress  or  field 
dependence  of  the  piezoelectric  properties.  Table  C-2  summarizes 
test  results  for  all  tests  performed  by  DWDL 
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Many  tests  wore  then  conducted  on  a  20-  disc  unbonded  stack  of  the 
same  material  with  each  disc  0.016  in.  thick  x  0.  5  in.  dia.  This 
stack  had  a  much  lower  capacitive  reactance.  By  varying  the  load 
resistance  and  measuring  output  power  at  various  stress  levels,  a 
dependence  of  stack  capacitive  reactance  on  stresgywas  experi¬ 
mentally  determined.  This  dependence  is  shown  in  Figure  C-7. 

Later  analysis,  particularly  the  time  dependent  formalism, 
related  this  dependence  to  the  electric  field  dependence  of  the 
dielectric  constant  as  previously  stated.  This  stack  was  tested 
at  0  H/.  and  stress  fluctuationsUo  8000  psi  at  various  loads. 

These  data  are  plotted  in  Figure  C-8.  The  stack  was  then  held  at 
8000  psi  and  RL  =  20  K  for  over  one  hour.  At  these  conditions, 
output  power  dropped  7",o  during  the  hour.  This  decrease  in  per¬ 
formance  during  the  initial  few  hours  of  testing  occurred  on  all 
subsequent  tests  as  well.  For  stress  fluctuations  to  and  including 
8000  psi  (1000  to  0000  psi)  me  decrease  is  reversible;  that  is, 
the  stack  returns  to  its  undegraded  condition  when  allowed  to  relax. 
For  a  stress  fluctuation  ol  12,000  psi  (1000  to  13,000  psi),  an 
additional,  permanent  degradation  occurred.  This  degradation 
was  observed  when  the  stack  was  stressful  overnight  at  12,000  psi 
and  60  Hz.  Output  power  across  the  20-K  resistive  load  decreased 
from  an  initial  value  of  1.  14  watts  (18.  2  v/in.  ’)  to  0.  79  watts 
(12.  8  w/in.3).  The  resistive  load  was  not  detectably  warmer  to  the 
touch,  nor  had  the  equipment  changed  in  temperature.  These  results 
are  also  shown  in  Figure  C-8. 

Work  input  measurements  attempted  on  these  tests  were  not  valid 
because  the  Bendix  displacement  measuring  unit  was  not  able  to  per¬ 
form  satisfactorily  at  60  Hz.  Both  signal  attenuation  and  phase  lag 
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Figure  C-7.  Effective  Capacitive  Reactance  of  Piezoelectric  Stack  as  a  Function  of  Applied  Stress 
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Figure  C-8.  Test  Results  for  Twenty-Disc  Piezoelectric  Stack 


■  were  introduced;  subsequently  the  Bently  Nevada  proximitor  and 
another  type  of  LVDT  (linear  voltage  differential  transformer)  trans¬ 
ducer  were  ordered.  The  latter  did  not  arrive  in  time  to  be  used  on 
tests  described  in  this  report. 

When  the  inductors  were  available,  the  same  stack  was  again  tested 
as  a  function  of  stress,  load  resistance  and  series  inductance.  The 
strong  dependence  of  output  power  on  series  inductance  is  shown  in 
Figure  C-9.  While  the  analytical  prediction  of  optimum  series  induc¬ 
tance  (X^  =  effective  Xc)  was  70  henries  with  the  time -ave raged 
formulation  described  in  Appendix  13,  the  prediction  was  50  h  with 
the  time -dependent  formalism  described  in  this  appendix.  Figure  C-9 
shows  the  experimental  optimum  value  was  45  h. 

This  stack  was  then  tested  overnight  at  a  stress  fluctuation  of  6000  psi 
at  60  H/.  with  a  2 0 - K  load  and  40-h  series  inductance.  The  results 
of  this  test  are  also  plotted  on  Figure  C-8.  Coincidentally ,  the 
time  -  0  operating  point  for  this  test  with  40-h  inductance  on  the  stack, 
which  had  degraded  during  the  test  at  12,000  psi  stress  fluctuation, 
corresponded  with  (he  original  short  test  at  6000  psi  stress  fluctuation 
and  no  inductance  prior  to  stack  degradation. 
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Figure  C-9.  Dependence  of  Output  Power  Upon  Inductance 


With  the  stack  at  an  optimum  inductance,  the  load  resistance  was  then 
reduced  until  the  maximum  power  density  of  tt.  0  w/in.  ^  was  achieved. 
This  is  shown  in  Figure  C-10. 

The  next  stacks  tested  were  10  disc  (1.  23  in.  dia  x  0.040  in.  thick 
disc)  stacks  of  P/.T-4  equivalent  piezoelectric  ceraniic  suppliec,  by 
Linden  Laboratories,  State  College,  Pennsylvania.  The  discs  in 
these  stacks  were  bonded  together  with  conducting  epoxy.  During 
the  first  few  minutes  of  testing  on  each  of  these  stacks,  arcing 
occurred  in  the  epoxy  on  the  sides  of  the  stacks.  This  occurred 
at  stress  fluctuations  of  6000  psi  as  series  inductance  was  added 
to  the  load  (thereby  increasing  voltage  across  the  stack).  Prior  to 
breakdown,  output  power  was  only  25%  of  the  predicted  value  for 
one  stack  and  50%  for  the  second  stack.  As  load  resistance  was 
varied  with  no  inductance,  maximum  power  occurred  at  the  predicted 
load  resistor  value.  The  probable  cause  of  the  low  power  output 
was  a  relatively  low  shunt  resistance  path  provided  by  the  epoxy. 

After  continued  arcing  on  both  stacks  at  stress  fluctuations  of 
3000  psi,  all  testing  of  these  bonded  1.25-in.  dia  stacks  was 
tc'rminated. 

Testing  resumed  on  the  20-  disc  (0.  4  in.  dia  x  0.  016  in.  thick  disc) 
stack  of  HDT-31  material  previously  tested;  Figure  C-l]  shows 
this  stack.  For  this  series  of  runs,  the  Bently  Nevada  probe  was 
utilized  to  make  accurate  displacement  measurements.  From 
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these  measurements,  an  efficiency  can  be  calculated.  Figure  C-12 
is  a  typical  force  versus  displacement  photograph  from  the  XY 
oscilloscope.  The  operating  conditions  for  this  case  were: 

f  =  60  Hz 


The  photograph  was  cut  up  to  obtain  the  area  with  the  force  displac- 
ment  loop.  The  area  was  weighed;  a  section  of  the  grid  display 
2  cm  x  4  cm  on  the  same  scope  picture  xerox  copy  was  then  weighed 
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Figure  C-10.  Power  Density  Dependence  Upon  Load  Resistance 
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Figure  C-12.  Typical  Force-Displacement  Experimental  Curve  for  Twenty-Disc  Piezoelectric  Stack 


From  these  two  weights,  the  true  nren  of  the  loop  was  established. 
Input  power  was  then  dele  r mined: 


P.  =  (3.  87  cnt4*) 


in 


0.  2- V 
cm  ) 


/>  ooo 

lb\ 

l  v 

) 

n- 


448  newtons 


lb 

]  meter 


* 


/0.02v\  /I  mil\  /9.001  in.\  /]  meter  \1 
V  cm  /  \  0.  2  vy  \  mil  7  \3  9.  37  inj]  5 


x  (60  sec  *) 


-  0.326  watts 


Efficiency  n 


0.  394 
0.  323 


100 


Current  uncertainty  in  these  measurements  includes  (])  ±10"o 
uncertainty  in  the  0.2  v/mil  ealibrn'tion  factor  for  the  Bently  Nevada 
displacement  transducer,  and  (2)  ±2%  in  determining  the  area.  The 
latter  was  estimated  from  three  attempts  to  determine  the  area  of 
the  same' curve  where  one  measurement  utilised  the  inside  of  the 
scope  trace  line,  the  second  used  the*  middle  of  the  scope  trace  line 
and  the  third  used  the  outside  of  the  line.  Other  uncertainties  which 
appear  to  be  small  compared  to  these  include  oscilloscope  calibration 
of  X  and  Y  axes,  the;  load  cell  calibration  factor  ( 1  000  lb/ volt) ,  and 
the  imperfect  sine  wave  affecting  voltage  measurement  as  indicated 
on  the  Ballentine  voltmeter.  Incorporating  these  uncertainties, 
efficiency  becomes 


73  ±12";, 
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In  the  near  future,  this  uncertainty  should  be  reduced  as  the  actual 
calibration  of  the  displacement  transducer  is  established. 

I'  igure  C-13  presents  the  efficiency  and -output  power  of  the  same 
20-disc  stack  as  a  function  of  inductance.  The  load  resistance  of 
30  K  gave  maximum  power  with  no  inductance  in  the  circuit.  Stress 
fluctuation  was  again  6000  psi  and  frequency  was  60  Hz.  For  the 
same  operating  conditions,  load  resistance  was  then  varied  while 
holding  load  inductance  near  optimum  (X^  =  effective  X  ).  These 
results  are  given  in  Figure  C-14. 

Testing  then  began  on  14-disc  unbonded  stack  (0.  50  in.  dia  x 
0.  10  in.  thick  disc)  of  HDT-3  1  material  at  60  Hz..  This  particular 
configuration  resulted  in  a  stack  capacitive  reactance  for  which 
neither  optimum  load  resistors  nor  series  inductance  was  available. 

At  60  Hz  and  6000  psi  stress  fluctuation,  the  stack  produced  the 
following  results: 

1.  4  2  watts  ( 5.  2  w/in.  at  L  -  0,  R.  =  163  K;  0  =  57% 

1.  8  watts  (6.  5  w/in.  .")  at  L  =  100,  R  -  163  K;  0  =  57% 

After  running  16  hours  at  these  conditions,  stack  output  had  stabi¬ 
lized  at  2.  0  watts  (7.  2  w/in.  with  the  stress  fluctuation  drifting 
and  eventually  stabilizing  at  6300  psi.  Returning  to  the  original 
stress  fluctuation  of  6000  psi  reduced  the  output  power  to  1.6  watts 
(5.7  w/in.  ■*).  This  decrease  from  the  initial  6.  5  w/in.  3  is 
qualitatively  and  quantitatively  consistent  with  the  results  shown 
in  Figure  C-8  for  the  small  stack. 

A  second  unbonded  s.tack  of  12  identical  discs  was  then  placed  on  top 
of  the  14  discs  to  produce  a  26 -disc  unbonded  stack  with  an  overall 
height  of  2.67  in.  As  in  all  unbonded  stack  tests,  the  interdl.se 
electrode  leads  were  0.002  in.  ec.ppor  sheet.  Plexiglas  sleeves 
aided  in  maintaining  alignment  during  set  up;  these  were  not  removed 
during  testing.  Testing  was  carried  out  at  120  Hz  and  6000  psi 
stress  fluctuation  for  various  load  conditions.  With  Xj  ~ 
effective  X  ,  the  load  resistance  was  reduced  to  maximize  output 
power.  These  results  are  shown  in  Figure  C  -  1  5. 

(3.3.3  Thermal ;  Effects 

These  we  re  the  first  tests  to  establish  apparent  thermal  effects.  For 
those  tests  producing  more  than  10W,  output  power  drifted  gradually 
upward  over  a  few  minutes.  For  example,  at  R.  10  K  at  t  -  0, 

V  -  375  v;  al  t  -  1  min,  -"-385  v;  at  t  --  2  min,  Vj  397  v.  The 
stack  may  have  been  heating’ as  a  result  of  the  increase  in  power,  or 
the  load  resistor  might  have  been  healing. 
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Figure  C-15.  Power  Density  Dependence  on  Load  Resistance  for  26-Disc  Stack 


At  a  power  level  of  7  v.  (14  w/in.  ),  the  stack  was  tested  lor  an  hour 
with  a  1'fo  increase  in  power  in  contrast  to  previous  observations  at 
lower  power  densities  (about  a  7%  decrease). 

Temperatures  were  measured  only  on  the  1.25-in.  dia  stacks.  All 
stacks  were  positioned  between  the  piston  and  load  cell  of  the  MTS 
system.  While  the  load  cell  remained  essentially  at  ambient  tem¬ 
perature,  piston  temperatures  increased  because  the  hydraulic  oil 
is  maintained  at  12  5SF  by  a  water-cooled  heat  exchanger  in  the  sys¬ 
tem.  The  top  of  the  piezoelectric  stack  was  72°F  and  the  bottom 
was  98'' F .  After  stress  fluctuation  at  6000  psi  and  60  Hz,  with  1  watt 
output  for  several  minutes,  the  top  of  the  stack  was  still  at  72 °F, 
but  the  bottom  had  warmed  ai  98.  5°F.  This  agreed  with  analytical 
predictions  based  on  calculated  electrical  losses  for  those  operating 
conditions. 

The  Plexiglas  sleeves  of  the  26 -disc  stack  were  warm  to  the  touch 
during  these  high  power  density  tests,  but  never  became  objection¬ 
ably  hot.  These  tests  included  maximizing  power  with  a  resultant 
efficiency  decrease  to  42".o.  This  is  a  worst  heating  condition  which 
would  not  normally  he  used  in  practice.  Because  of  the  small 
diameter  discs,  relatively  large  copper  leacis  which  could  also  act 
as  heat  conductors,  and  the  Plexiglas  sleeves,  no  thermocouples 
were  attached  to  the  stacks  during  these  tests. 

Suppliers  literature  (References  C-l  and  C-2)  provided  materials 
property  data  as  a  function  of  temperature.  The  relative  dielectric 
constant  increases  as  the  temperature  rises.  This  change  appears 
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to  dominate  in  performance  calculations;  therefore,  power  is 
expected  to  increase  with  temperature.  However,  tne  maximum 
stress  allowable  drops  fr  om  12,000  psi  at  2 a°C  to  6000  psi  at 
lOO^C.  If  6000  psi  had  been  chosen  for  loriff- term  operation  at 
25°C,  then  a  much  lower  value  should  be  chosen  for  100°C  oper¬ 
ation.  This  interaction  between  stack  operating  conditions  and  ther¬ 
mal  design  considerations  remains  critical  to  successful  system 
operation.  Furthermore,  long-term  testing  with  accurate  temper¬ 
ature  measurements  under  varied  conditions  will  be  required  to 
unequivocally  determine  allowable  operating  conditions.  Tempera¬ 
ture  effects  under  realistic  operating  conditions  have  been  negli¬ 
gible  to  date  as  predicted  bv  a  simplistic  thermal  analysis  assum¬ 
ing  all  heat  is  generated  at  the  center  of  the  stack.  Should  a 
temperature  problem  with  full-size  stacks  arise  in  the  future,  stack 
cooling  is  a  technically  feasible  alternative,  as  previously  discussed. 
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Appendix  D 

DETERMINATION  OF  HEAT  TRANSFER  COEFFICIENTS 
AND  REGENERATOR  PERFORMANCE  IN  THE  STEPZ  ENGINE 

The  heat  exchangers  and  regenerator  in  the  STEPZ  engine  are  quite 
different  from  those  typically  used  in  Stirling  engines.  A  partial  cross 
section  of  the  annular  flow  passage  between  the  cylinder  and  displacer 
walls  is  shown  in  Figure  D-l.  Different  regions  of  the  annulus  serve 
as  heater,  cooler,  and  regenerator  in  this  approach.  There  is  a 
lower  surf ace-to- volume  ratio  and  less  heat  transfer  area  than  for 
conventional  Stirling  engines.  This  imposes  a  power  density  limita¬ 
tion  in  exchange  for  decreased  gas  dead  volume  and  less  mechanical 
complexity. 

Heat  exchanger  performance  is  readily  analyzed  for  the  STEPZ  engine 
geometry.  An  excellent  treatment  of  all  essential  aspects  of  STEPZ 
heat  exchanger  analysis  is  given  in  Reference  D-l,  especially  chapters 
2  and  6.  Heat  transfer  coefficients  are  determined  on  the  basis  of 
laminar  flow  with  fully  developed  velocity  and  temperature  profiles. 
These  are,  in  general,  valid  assumptions  which  in  any  case  give  con¬ 
servative  results  if  they  are  not  precisely  true  during  transient  con¬ 
ditions.  Specifically,  for  the  reference  design  case,  the  Reynold's 
number  varies  from  zero  to  5600,  with  a  mean  value  of  1600.  Turbu¬ 
lent  flow  may  possibly,  but  hot  necessarily,  be  established  for  momen¬ 
tary  Reynold's  numbers  above  2000.  If  turbulent  flow  at  a  high  Rey¬ 
nold's  number  does  occur,  the  heat  transfer  coefficient  will  momen¬ 
tarily  be  increased  by  up  to  a  factor  ot  three.  Similarly,  near  the  ends 
of  the  displacer  where  velocity  and  temperature  profiles  are  not  fully 
established,  the  heat  transfer  coefficient  is  higher  by  a  factor  of  two 
to  three.  This  length  is  typically  10  to  100  flow  widths  or  0.05  to  0.5 
in.  in  the  reference  engine.  A  third  area  where  a  conservative  heat 
transfer  assumption  is  made  is  in  the  end  regions.  Stagnant  gas  con¬ 
duction  is  assumed,  although  Substantial  mixing  certainly  occurs. 

This  approach  gives  confidence  that  the  heat  exchanger  model  will  not 
indicate  unrealistically  Vngh  performance. 


Mathematical  Relationships: 

The  Nusselt  number  Nu  is  defined  as 
Nu 

where 


jL2b 

K 


2hG 

K 


(D-l) 


heat  transfer  coefficient 


Preceding  page  blank 
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~  hydraulic  diamotor  -•  2G  for  a  narrow  annulus 
“  radial  clearance  between  displacer  and  cylinder 
K  thermal  conductivity  of  gas 


In  practice,  it  is  the  heat  transfer  coefficient  that  is  required  because 
the -Nus selt  number  is  constant  for  a  given  region  under  the  assumed 
operating  conditions.  The  equation  of  interest  is,  therefore, 


h 


Nu  K 
2G 


(0-2) 


Thermal  conductivity  of  th working  fluid  is  a  strong  function  of  temp¬ 
erature  and  must  be  separately  determined  for  each  region  or  node. 
The  Nusselt  number  is  dependent  on  geometrical  and  thermal  charac¬ 
teristics  ol  the  heat  exchangers.  The  a  nnular  gap,  which  is  very  small 
compared  with  the  diameter,  is  essentially  identical  with  flow  between 
infinitely  wide  parallel  plates  separated  by  a  distance  G.  In  the  heater 
and  cooler  regions,  boundary  conditions  are  essentially  isothermal  on 
the  outside  of  the  gap  and  adiabatic  on  the  inside.  From  Figure  6-1  of 
Reference  D-l,  this  leads  to  a  Nusselt  number  of  4.86  for  healer  and 
cooler.  In  the  regenerator  region,  a  large  amount  of  heat  transfer  to 
and  from  the  gas  occurs  regenerat ively  on  each  cycle,  although  the 
spatial  temperature  profile  remains  essentially  constant.  This  con¬ 
dition  represents  a  constant  heat  transfer  rate  per  unit  of  length  on 
both  inner  and  outer  surfaces,  for  which  the  Nusselt  number  is  8.24. 


This  approach  yields  conservative  beat  transfer  coefficients  for  the 
STEP'/,  engine.  These  coefficients  are  then  used  to  determine  engine 
there*'  ’  performance  as  outlined  in  the  paragraphs  which  follow.  Meat 
tr;  s  rate,  Q,  in  the  heater  or  cooler  is  given  by 


Q  -  h  A  AT  (D-3) 

where 

A  -  heater  or  cooler  area  <>n  the  outside  of  the  annulus 

AT  temperature  difference  between  the  heater  or 

cooler  and -adjacent  gas 
y 

Meat  transfer  in  the  hot  end  is  given  by  the  sum  of  all  nonpa rasitie  heat 
requirements  for  tile  engine,  while  in  the  cooler  it  is  this  value  loss 
t h o  mechanical  work  output.  Equation  (D-3)  can  then  be  solved  for  the 
gas  temperature  in  the  heater  and  cooler  regions.  These  temperatures 
then  affect  the  work  output,  Carnot  heal  input,  and  reheat  losses,  such 
that  heat  transfer  requirements 'are  revised  and  new  gas  temperatures 
determined  on  an  iterative  basis  until  a  convergence  criterion  is  sat¬ 
isfied. 
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Regenerator  performance  is  determined  using  an  effectiveness  approach 
as  m  Reference  D-l,  Chapter  2  and  mere  directly  in  Reference  D-2. 

The  regenerator  analysis  is  again  conservative,  in  that,  gas  specific 
heat  at  constant  pressure,  cp,  is  used  to  calculate  reheat  loss.  The 
effective  specific  heat  is  some  intermediate  value  between  c  and  cv 
(specific  heat  at  constant  volume).  This  true  sp>eci:E^c  heat  varies  in  an 
indeterminate  manner  over  a  cycle.  Using  the  larger  c  at  all  times 
yields  a  larger  reheat  loss  which  provides  for  some  analytical  perfor¬ 
mance  margin. 


Tile  STEP/',  regenerator  is  a  periodic  counterflow  type  with  alternating 
flow  of  the  same  gas  in  opposite  directions.  For  these  conditions,  the 
following  relationships  are  valid. 


( h  r  A  r )  c 

and 

(tinr  cp>c 

wher  e 


mr  - 


(hr  Aj,)^  -  hr  Ar 


(mr  cp)h  =  rhr  c. 


(D-4) 

(D-5) 


heat  transfer  coefficient  in  regenerator  region 
total  surface  area  in  regenerator  region 
mass  flow  rate  through  regenerator 


Cp  =  working  gas  specific  heat  at  constant  pressure 

and  the  subscripts  c  and  h  refer  respectively  to  the  gas  flowing  from  the 
cold  region  and  the  hot  region.  Relationships  needed  to  evaluate  regen¬ 
erator  effectiveness  are  found  in  Reference  D-2  as  follows. 


NTU  - 

h  r  Ar 

t 

mr  Cp 

NTUU  = 

NTUC  1  4 

where  the  last  equality  results 

f  -  NTUq 

.1  +  NTUC 


(hr  Ar)c 
(hr  Ar)h 

from  applying 


=  \  NTU 

Z  c 

equation  (D-4). 


(D-6) 

(D-7) 

F  urther, 
(D-8) 


where 


NTUC  -  number  of  trcinsfer  units  for  cold  flow 
NTU(J-  overall  number  of  transfer  units 
f.  -  regenerator  effectiveness 
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(D-9) 


The  eti  ectiveness  may  also  bo  expressed  as 


who  re 


AT, 


AT, 


AT 


AT 


overall  temperature  difference  between  hot  region 
gas  and  cold  region  gas 


AT, 


temperature  difference  between  hot  region  gas  and 
gas  entering  the  heater  from  the  regenerator 


ATj^  is  the  temperature  rise  which  must  be  accomplished  in  the  heater 
as  a  result  of  regenerator  ineffectiveness,  which  gives  rise  to  the  main 
reheat  loss.  The  two  expressions  for .  regene  rator  effectiveness  may  be 
equated  and  solved  for  ATu. 

4  rp 

AT,  =  _ o _  (D- 10) 

1  +  NTU 

o 


Using  Equations  (D-6)  and  (D-7)  in  (D-10)  reduces  AT^  to  an  expression 
in  more  basic  quantities, 


2  mr  cD  ATq 


2  rhr  cp  +  hr  Ar 


(D-  1 1) 


The  reheat  power  which  must  be  added  to  achieve  this  temperature 
rise  is 


P 


rh  =• 


•n  r  cp  ATh 


(D-  12) 


where  the  factor  of  1  / 2  is  included  because  the  gas  is  only  flowing  to¬ 
ward  the  hot  region  half  the  time  as  opposed  to  the  continuous  flow 
case  in  Reference  D-2.  Combining  Equations  ( D - 1 1 )  and  (D-12)  gives 
the  final  result  for  main  reheat  loss. 


(mr  cp)2  AT0 
2  mr  cp  +  hr  Ar 


(D  -13) 


The  basic  heat  exchanger  and  regenerator  analysis  summarized  here  is 
common  to  both  the  isothermal  and  nonisothe rmal  computer  programs. 
The  method  of  implementing  the  analysis  differs,  and  the  isothermal 
program  typically  indicates  significantly  lower  reheat  losses  and  temp¬ 
erature  drops.  This  is  shown  in  Table  2-2.  The  reason  for  the  reheat 
loss  discrepancy  is  that  the  isothermal  program  calculates  it  from 
time  averaged  flow,  temperature  drop,  and  heat  transfer  coefficient  in 
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integration  of  EquatLm^L^TVve^LTy'c ^"’“The^"  US<,S  a.num<'ri<:al 

at  times  wh7„  ^  ^  f’”r e™aX™rnbl  *  rt^ 

on  [y-ovplanat  ion  for  the  cioseoo  rreia,  ion  of 

>17,  flTY"  Tabl°  f°r  th°  lW°  m"'la's  «  that  most  heaf  tranTL 
occurs  under  favorable  conditions. 

In  summary,  the  mechanically  simple  but  unconventional  heat  exchanger 
I  a.PPro«h  »•-<•  »"  the  STEF/,  engine  can  be  rea.mj  " 

e  ce  D  1  Wb  y  experimentally  verified  analysis  given  in  Refer- 

e  ers  n  ’  .  Iv  M/  T'*"™  <'xis,‘  rr-la,iva  to  appropriate  pa, -am- 

ond  it  ions,  the  most  conservative  approach  is  taken.  While 

imit??;rr°n  °'  th|  r‘'S““S  *°  th£'  iSathormal  P'-OBf»m  is  necessarily 
imitec.  in  accuracy  by  assumptions  in  the  isothermal  model,  the  non- 

isothorma1  simulation  should  offer  a  high  degree  of  confidence  in  the 
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t  *  Appendix  E 

SAMPLE  DESIGN  OF  A  LINEAR  ELECTROMAGNETIC 
DISPLACER  ACTUATOR 

_ _ -  'ir 

A  dynamic  coil  linear  electromagnetic  displacer  actuator  has  been 
designed  at  DVVDL  for  a  66  w(m)  STEPZ  engine  module  for  space  ap¬ 
plications.  The  configuration  of  the  actuator  is  shown  in  Figure  E-l. 

A  weight  analysis  indicated  that  the  magnetic  circuit  containing 
Alnico-9  was  lighter  than  one  with  samarium-cobalt,  as  a  result  of 
flux  concentration  in  the  center  pole.  Advantageous  use  of  SmCoj  is 
anticipated  for  higher  windage  engine  modules.  The  actuator  for  the 
66  w(m)  engine  module  produces  1  w( m)  at  60  Hz  with  an  efficiency 
of  27%.  The  actuator  weighs  0.  5  lb  and  requires  1.  9  amperes  (rms) 
excitation  for  rated  output.  The  impedance  of  the  coil  is  2.  0f2with  a 
lagging  phase  angle  of  59°. 

E.  1  DESIGN  PHILOSOPHY 

For  space  applications,  reliability  of  the  STEP/,  displacer  drive  must 
be  of  the  same  order  as  that  of  the  displacer  suspension  and  other 
engine  components.  The  intrinsically  high  reliability  of  flexural 
mounting  may  be  more  than  an  order  of  magnitude  superior  to  that  of 
couplings  involving  rotating  motors  and  mechanical  linkages.  The 
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Figure  E-1.  Displacer  Actuator  for  STEPZ  Engine 
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oscillatory  displacer  motion  suggests  the  use  of  a  linear  motion  elec¬ 
tromechanical  actuator,  such  as  a  solenoid  or  a  dynamic  coil  (loud¬ 
speaker)  motor  system.  Both  devices  permit  the' suspension  of  the 
oscillating  component  directly  on  the  displacer  and  require  no  seals 
or  mechanical  Icadthroughs  to  function. 


Oscillating  linear  motion  is  produced  by  the  interaction  of  alternating 
currents  and/or  magnetic  fields.  The  most -efficient  use  of  magnetic 
material  in  a  solenoid  is  achieved  with  direct-current  excitation  of  the 
field  coil.  With  ac  excitation  (Reference  E-l),  the  magnetic  circuit 
is  used  only  one  half  as  effectively  as  that  of  a  dc  system  to  produce 
the  same  average  force  on  the  armature.  As  the  armature  is  a  com¬ 
ponent  of  the  magnetic  circuit,  its  mass  reflects  this  inefficiency  and 
contributes  a  several-fold  increase  in  the  weight  of  the  displacer. 

The  force  acting  on  a  solenoid  armature  is  a  monodi rectional  force  of 
attraction.  Use  of  a  solenoid,  therefore,  also  necessitates  incorpora¬ 
tion  of  a  spring  or  double  solenoid  to  provide  the  restoring  force  for 
oscillatory  motion. 

The  dynamic  coil  linear  motor  depends  on  the  interaction  of  an  alter¬ 
nating  electric  field  and  a  permanent  magnetic  field  for  its  action. 

The  instantaneous  force  (F)  on  a  coil  system  carrying  a  sinusoidally 
varying  current  is  given  by 

F  --  (B  x  I  )  2sinwt  (E-l) 

nr 


where 


B  -  flux  density  in  permanent  magnet  gap 

1  ~  maximum  coll  current 
m 

2  -  length  of  coil  winding 

u>  -•  angular  frequency  of  coil  current 

Oscillatory  motion  of  the  coll  results  from  the  sinusoidal  force- 
function.  This  response  is  significantly  superior  to  that  of  the  sole¬ 
noid.  The  dynamic  coil  also  follows  complex  excitation  waveforms 
which  permit  norisinusoidal  excitation  if  desired  to  effect  sophisticated 
coni  :  ol  functions.  For  these  reasons,  therefore,  a  dynamic  coil  drive 
is  considered  the  prime  choice  to  activate  the  displacer. 


E.  2  DESIGN  ASSUMPTIONS  AND  CONSTRAINTS 

The  design  of  the  displacer  actuate’  was  based  on  the  following 
'assumptions. 

1.  Actuator  power  supplies  the  displacer  windage  loss  in  an 
otherwise  dynamically  resonant  mechanical  system. 

2.  Coil  mass  added  to  the  displacer  has  negligible  influence 
on  the  displacer  resonant  response. 
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'■  ^agneUc  flux  fringing  the  motor  gap  i-s  the  same  order  as 
the  flux  between  pole  faces.  1  as 

4.  The  internal  field  energy  of  the  magnet  provides  the  nan 

?sVc7l  a"d1frLnSC  field;  magnetization  of  the  pole  pieces 
is  neglected.  1 

5.  The  engine  working  fluid  absorbs  the  Joule  heating  produced 

by  the  coil  current.  1 

Assumed  characteristics  of  the  displacer  were: 


1.  Windage  loss: 

2.  Stroke: 

3.  Diameter: 


1  w(m)  at  60  Hz 

1.  27  mm  (0.  05  inch) 

2.  54  cm  ( 1  inch) 


The  mechanical  work  done  by  the  displacer  in  one  half  cycle  is  the 
product  of  average  force  (F^)  and  stroke  (s).  From  Equation  E-l, 


B  I  I?  rT‘ 

m  s  / 


f  . 

/  sin 

hi 


(E-2) 


0.  637  B  I  l 
m 


(E-3) 


where  0  is  the  angle  in  radians  from  the  instant  of  z 


ero  current. 


At  60  Hz,  there  are  120  half  cycles  per  second;  therefore,  mechanical 
power  of  the  displacer  (P^)  is  given  by 


120  F  s  Joums/sec 

3.  V 

76.  38  I  Bfs  wjml 
m  ' 


(E-4) 

(E-5) 


With  s  =  1  27  mm  and  P  =  1  w(m)  Equation  E-5  leads  to  an  expression 

t  I  V  1  »■%  1-  Vr  i^v  f  T3  1 


ofT  (  in  term's  of  B 


10.  31 


(E-6) 


Joule  heating  (P^)  in  the  dynamic  coil  is  expressed  by 


m  P  f 


(E-7) 


a 
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f 

f 


f 

i 


where  P  is  the  resistivity  of  coil  material  and  a  is  the  cross  sectional 
a  roa. 


With  reference  to  Equation  E-6,  and  for  a  given  coil  geometry, 


P 

e 


K 


(E-H) 


In  Figure-  E-2  Joule  heating,  Pc,  is  plotted  versus  gap  flux  density,  B, 
in  the  range  0.  1  <  B  <  1  webers  per  square  meter  where  permanent 
magnet  materials  of  interest  achieve  their  maximum  energy  product. 

The  flux  density  in  the  working  gap  was  selected  in  the  range  0.  3  to 
0.  5  Wb/m^  corresponding  approximately  to  the  range  in  which  Pc  is 
between  3  and  1 0°!  of  its  value  where  B  --  0.  1  Wb/m^.  This  does  not 
provide  an  optimum  configuration  and  was  pursued  on  a  quasi-intuitive 
basis  to  produce  a  compact  magnetic  circuit  with  acceptable  weight 
and  performance. 

A  major  region  of  uncertainty  occurs  in  estimating  flux  distribution  in 
the  working  gap.  Computer-aided  mapping  of  this  region  is  necessary 
to  refine  this  effort,  the  fringing  flux  (i.  e.  .  flux  not  linking  the 
dynamic  coil)  was  considered  equal  to  the  useful  flux.  An  aspect  ratio 
of  5  was  selected  for  tint  working  gap  to  maintain  the  general  validity 
of  this  assumption. 

Two  configurations  of  the  magnetic  circuit  are  shown  in  Figure  E-3 
consistent  with  a  design  which  suspends  the  dynamic  coil  on  an  exten¬ 
sion  of  the  displacer.  The  placement  of  the  permanent  magnet  shown 
in  A.  Figure  E-3  is  the  more  desirable.  However,  with  gap  dimen¬ 
sions  shown  in  Figure  F.-l  and  flux  density  in  the  gap  approaching 
0,  •  W  |,  /  m“ .  the  geometry  demands  flux  densities  in  the  center  pole 
ere.  ter  than  I  Wb/m*".  In  Figure  1  of  Reference  E-2,  0.  75  Wb/ni^ 
is  the  maximum  flux  density  at  which  an  acceptable  energy  product  is 
aehi<  veil  with  available  materials.  For  Alnico  9,  this  is  in  the  region 
of  its  1  P»I  l|mux  point.  However,  tor  SniCoj,  |  Bli  {max  occurs  at. 
approximately  (>.  4-  \Vi;/m',  The  magnetic  ci  rcuit  design  is  forced, 
therefore,  to  that  of  the  "outboard"  magnet  (B.  Figure  E-.3)  to  permit 
sizing  the  magnet  independent  ol  displacer  diameter.  This  geometri¬ 
cal  constraint  will  be  superimposed  on  both  optimum  and  sub-optimum 
designs  in  all  engines  with  a  displace  r  diamete  r  less  than  approxi  - 
mately  1 .  -1  in.  j 

Detailed  calculations  were  performed  of  the  Alnico  9  "outboard" 
magnet /coil  configuration,  and  a  similar  design  was  completed  for  a 
SmCoc;  magnet  ad  rcuit  to  provide  the  same  output.  Table  1M  sum¬ 
marizes  a  weight  comparison  of  the  two  systems. 


t 

1  : 
!.  1 
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Table  E-l 


WEIGHT  SUMMARY  OF  DISP FACER  ACTUATOR  COMPONENTS 
_ .  W*** 


Component 

Weight  (g) 

SmCoc  at  0,  45 
b 

Wb/m*'  Alnico  at  0.  75  Wb/m^ 

Magnet 

34.  0 

52.  5 

Inner  Pole  (1)  (2) 

110.  0 

100.  0 

Outer  Pole  (2) 

ioy.  o 

63.  5 

Coil 

9.  3 

9.  3 

Displacer  Extension 

1.  5 

1.  5 

Total 

262.  8 

227.  0 

(1)  With  central  conical  void 

(2)  With  ellipsoidal  transitions  from  magnet  element 


The  lighter  magnet  is  achieved  with  the  Alnico  9  despite  Us  lower 
energy  density  as  a  result  of  |BH|max  oc  curring  at  0.  75  Wb/m2.  The 
superior  energy  product  of  SmCoj  would  be  reflected  more  favorably 
in  a  design  for  a  larger  displacer.  Table  E-2  summarizes  character¬ 
istics  of  the  Alnico  9  actuator. 


E.  3  CONCLUSIONS 

If  the  assumption  of  flux  distribution  is  realistic,  this  design  exercise 
provides  encouragement  for  the  application  of  electromechanical  linear 
device  actuation  and  control  of  the  STEP/,  engine.  The  particular 
engine  investigated  imposed  an  unscalable  constraint  by  requiring  the 
center  pole  to  fit  inside  the  extension  of  a  small  displacer.  The  re¬ 
sulting  high  flux  density  excludes  the  favored  configuration  of  the 
magnetic  circuit.  Alnico  9  provides  the  lightest  magnetic  circuit  of 
any  available  material  working  at  its  BHniax  point.  The  superior 
energy  product  of  .SmCo5  would  benefit  magnetic  circuit  designs  for 
larger  displacers.  Its  availability  is  established  and  it  is  considered 
both  performance  and  cost  Competitive  with  Alnico  materials. 

Flux  mapping  of  the  working  gap  is  necessary  as  a  basis  for  optimiz¬ 
ing  the  configuration  of  the  linear  actuator.  In  particular,  flux  con¬ 
centration  should  be  traded  against  increases  in  fringing  flux  for  lower 
aspect  ratio  gaps. 


Table  E-2 


ALNrco  «  displacer  actuator  characteristics 


Characteristic 

Weight  (g) 

Mechanical  Output  (w  } 

nr 

Specific  Mechanical  Power  (w  /lb) 

Coil  Dissipation  {w  ) 

0 

Efficiency  (%} 

Coil  Impedance  (ft) 

Coil  Phase  Angle  ( °) 


Value 

22?  (0.  5  lb) 


*3.  75 
27 

2.  0  (1.  0  +  j  1.7) 

59 
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